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CHAPTER 1 INTRODUCTION 


This report documents the research accomplished from May 
16 , 1976 through July 31, 1977 under the funding of NASA Grant 
NSG-2145 , The research conducted was concerned with the design 
optimization of short haul and connnuter airplanes as proposed in 
References 1 and 2. The intent of the research was to look at 
the problem of commuter airplane configuration design for the 
minimization of Direct Operating Costs (DOC) . A more detailed 
explanation of the purpose and objectives of this project is 
provided in the following section. The status of the research 
and finances of the project are discussed in Section 1.2. 

1.1 Purpose and Objectives 

It was proposed to look at the problem of commuter airplane 
configuration design from the following view points: 

’ Assume that a specific cabin volume is needed to meet 
some mission criterion. This will be called a utility 
constraint. The investigation will focus on commuter 
type airplanes with a crew of two and a passenger load 
up to approximately thirty. 

• Assume that specific stability and control requirements 
must be met. These will be called stability and control 
constraints. 

• Assume that specific mission profiles must be flown. 

These will be called mission constraints. 

• Assume that specific performance (for example, field length, 
minimum speed, single-engine climb, etc.) requirements must be met. 
These will be called performance constraints. Attention 

will be focussed on field lengths in the 2500 ft. to 
4500 ft. range. 

The problem then was to design an airplane configuration which 
would have the lowest DOC. Figure 1.1 illustrates the problem . 

The approach taken to the solution of this problem was along the 
following lines: 


1.1 


Given: 


1 ) 


o o o o o 


Cabin -utility Constraint 


2 ) 



Mission Constraint 


3) Stability and Control Constraints 

4) Performance Constraints 



The Configuration Which Minimizes D.O. C. 


Figure 1.1 Illustration of the Design Problem 
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a) Investigate methods for determining the fuselage shape 

and size which would yield minimum fuselage drag and ] 

weight. The fuselage was to be treated as three separate 1 

sections (the nose, the cabin, and the tail cone) , and ij 

methods for the design of each section were to be considered. ! 

Chapter 2 discusses the configuration design approaches } 

considered for the fuselage. " 

b) Determine how the methods of item a) might be integrated ( 

into the NASA-Ames General Aviation Sjnithesis Program 

(GASP) . The GASP was chosen to be the most effective i 

way to apply configuration design methods to the overall 
preliminary airplane design process. Chapter 3 discusses 

the ways in which the methods of item a) were to be | 

integrated into the GASP. 

c) Determine the critical stability and control constraints 

which needed to be integrated simultaneously into the GASP. i 

Chapter 5 discusses the stability and control constraints 
that were considered and explains the ways in which they 
were to be integrated into the GASP. 

d) Perform a short study of wing sizing methods to consider 
trade-offs in wing loading for effects on performance, 
ride qualities, etc. One method in particular that was to 
be considered for wing sizing was the method proposed in 
Reference 3. Also to be considered was a rational approach 
to wing placement, bearing in mind weight and balance, 
stability and control, passenger acceptability, and feasibility 
of associated structural arrangements. These studies are 
discussed in Chapter 3. 

e) Consider the problem of estimating the dr'ag of airplane 
components submerged in a propeller slipstream. This 
problem is discussed in Chapter 4. 

As has been stated, the GASP was to be the medium through which " 

all of the configuration design methodologies were to be applied. 

Unfortunately particular difficulty was encountered in getting 
a version of the GASP to run on the University of Kansas Honeywell 
66/60 computer. For this reason it has not been possible to com- 
plete the integration of the design methodologies into GASP. A 
complete discussion of the transliteration process is presented in 
Chapter 5. 
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1.2 Research and Financial Status 


1.2.1 Research Status 

At the time of this report much of the work outlined in the 
objectives has been completed. The fact that the GASP is still 
inoperational has provided the greatest set-back. Revisions to the 
GASP to include the newly developed methods could not be made with- 
out first being assured that the GASP could be operated properly in 
its existing form. Also^ more research needs to be conducted in the 
area of stability and control methods. Table 1.1 summarizes the 
status of each research area. Also, Table 1.1 refers to the -chap, tees 
where each facet of the research will be discussed. 

1.2.2 Financial Status 

Tables 1.2 and 1.3 present the budgets for the two phases of the 
project as proposed in References 1 and 2. Together, they represent 
a total allocation of $50,590 through NASA Grant NSG-2145. Although 
a complete breakdown of expenses through July 31, 1977 was not 
available at the time of this report, it was possible to arrive at 
an estimate as to the remaining balance as of that date. It is 
estimated that as of July 31, 1977 the balance of funds remaining in 
NASA Grant NSG-2145 will be approximately $100. 
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TABLE 1*1 DESIGN OPTIMIZATION GRANT STATUS 
31 JULY 77 



ITEMS AMD/OR OBJECTIVES 

ITEMS 

ITEMS 

ITEMS 


PROPOSED 

ACCOMPLISHED 

PENDING 

TO BE PROPOSED 

PROPOSAL 1 

PROPOSAL 1 1 



PROPOSAL 1 1 i 


I) Procedure to find fuselage/ 
empennage shape for min* 
drag 

a)Roskam/ri Mman Method 


b)Cabtn Arrangement and 
Baggage Compartment 
Studies 


c)Wetted Area Studies 


2)Weight Estimation 
Methods 

a) Torenbeek * s Method 


OContInued Cabin arra. ;, 

& Baggage Compartment Studies 


' 1)KU-FRL 902; Program 
written and tested. 
Results questionable 
Chapter 2, this report 

2)KU“FRL 901; Prelimin- 
ary correlations 
Chapter 1, this report 


3)KU-FRL 901; Prelimin- 
ary correlations 
Chapter this report 


If)KU-FRL 902; included 
in Boskam/Fi 1 Iman 
Method Program 
Chapter 2* this report 


l)Test Roskam/FI 1 Iman 
Method with Wing 
Balance included in GASP 


2) Compar i sons with 
McDonnell Douglas and 
Boeing Vertol Methods 
Chapter 2, this report 

3) Correlat ion of Actual 
Wetted Areas with 
FliSE-computed Areas 
Chapter this report 


4}Checking use of GASP 
Weight Estimation 
Me thods 
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3)StabnitY Constraints 

2)Ana lysis oF GASP designs 

5)KU”FRL 902; Empennage 

Sizing by C £ C 

fn n„ 

a 0 
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est, to enable dynamic 
stabi 1 ity analysis, 
and to locate wing 

OChapier 6, thir report 

^i) 5 i i pst ream Drag 
Analysl s 


6)KU-FRL 902 
Chapter h 



S) Include Methods Into 
GASP 

3)Foselage Config. Studies 

7)KU-FRL 9D2; Flowcharts 

6)Updat1ng Cessna version 
of GASP to NASA version 
Chapter 5. this report 



a)Nose section model 

8)KU-FRL 313-3; FUSE 

7) Use of FUSE wi th NC5U 
BODY Program to Analyze 
Nose 



b)Cabrn sizing 

9) Same as Above; Design 
Mode of FUSE 
Chapter 2, this report 




TABLE 1,1 DESIGN OPTIMIZATION GRANT STATUS (CONT’D) 


ITEMS AND/OR OBJECTIVES 

ITEMS 

ITEMS 

ITEMS 

PROPOSED 

ACCOMPLISHED 

PENDING 

TO BE PROPOSED 

PROPOSAL 1 PROPOSAL 1 1 



PROPOSAL Ml 


^)Wing Config. 




a)Si£mg 

tO) Chapter 3t this report 

8)Loftin*s Method to 
be Considered 


h) Placement 

lOChapter 3, this report 

9)Considering combina- 
tion of GASP and 
Torenbeek‘s method 


5)Deve1op algorithms 
for GASP 

M 

12)Chapter 3. this report 

10)Finding that It Is 
possible to simply 
amend GASP methods 

2)Complete addition of 
developed algorithms 
Into GASP and check 
them out 




Table 1.2 Proposed Budget 
May 16, 1976 - December 15, 1976 


Salaries & Wages 

NASA 

KU 

TOTAL 

Principal Investigator (Rosk£\m) 
20% for 4 mos. academic* 

67% for 3 mos. summer 

1,271 

5,807 

1,271 

0 

2,542 

5,807 

Research Assistants (3 1/2) 
50% for 7 months 

9,800 

0 

9,800 

Secretary 
2 man-months 

1,000 

0 

1,000 

Total Salaries & Wages 

17,878 

1,271 

19,149 

Fringe Benefits 




16% Faculty & Staff 
7% Students 

1,293 

686 

203 

0 

1,496 

686 

Computer 

2,000 

0 

2,000 

Supplies & Reproduction 

200 

0 

200 

Telephone 

200 

0 

200 

Travel {2 trips to Ames) 

750 

0 

750 

Total Direct Costs 
# 

23,007 

1,474 

24,481 

Indirect Costs - 53.6% of 
Total S & W 

9,583 

681 

10,264 

Total Proposed Costs 

32,590 

2,155 

34,745 

*Principal Investigator’s academic 
equally- Each paying only 10% for 

time cost shared by KU and NASA 
4 mos. academic time = 20%., 

Salary Schedule FY76 

FY77 



Roskam 2,889 /m.o . 
Research Assistant B90/rao. 
Secretary 500/mo. 

3,178/mo 
• NA 
NA 

- 
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Table PROPOSEO BUDGET 

January 1, 1977 - June 30. 1977 


Pi rect Costs 

Salaries & Wages 

Principal Investigator - Roskam 
25% for 1 mo, academic FY 77 
15% for 1 mo. summer FY 77 

Graduate Research Assistant 
15% for h.5 mos. academic 77 
100% for 1 mo. summer 77 
75% for .5 mo. summer 77 

Undergraduate Research Assistants 
75% for 5 mos. academic 77 {1 student) 
100% for 1.5 mos. summer 77 (1 student) 
100% for .5 mo. summer 77 (l student) 

Secretary 

2.33 man-months (student) 

Total Salaries and Wages 

Fringe Benefits 

17% staff 
7% students 

Other Direct Costs 
Computer 

Supplies and Reproduction 

Telephone 

Travel 

Total Direct Costs 

Indirect Costs @ 53-6% of Salaries 
6 Wages 

Total Proposed Costs 

Less unexpended balance in NASA 
Grant ^'SG 2145 

TOTAL PROPOSED COSTS 


NASA 

KU 

TOTAL 

$ 416 

$ 417 

$ 833 

2,500 

0 

2,500 

3,038 

0 

3,038 

900 

0 

900 

337 

0 

337 

2,625 

0 

2,625 

1,050 

0 

1,050 

350 

0 

350 

1,398 

0 

1,398 

12,614 

417 

13,031 

496 

71 

567 

679 

0 

679 

2,540 

0 

2,540 

400 

0 

400 

500 

0 

500 

1 ,000 

0 

1 ,000 

18,229 

488 

18.717 

6.761 

224 ■ 

6,385 

$ 24,990 

$ 712 

$ 25,702 

6,990 




$ 18,000 
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CHAPTER 2 FUSELAGE CONFIGURATION STUDIES 


As stated in Chapter 1 one of the objectives of the research 
program was to investigate methods for determining the fuselage 
shape and size which would yield minimum fuselage drag and weight. 

A two“phase approach was taken to meet this objective. The first 
phase consisted of evaluating existing methods for sizing the fuselage 
or, where the need arose, to develop new methods. Also included in 
this phase was an evaluation of existing methods for determining 
fuselage drag and weight. The second phase of this approach was to 
integrate the methods chosen as a result of the first phase into the 
GASP and to perform trade-offs in the fuselage configuration design 
to find trends from which optimum configurations might be determined. 

In this manner it was also hoped that it might become possible to 
locate critical configuration determining parameters, and eventual- 
ly to develop routines to optimize the fuselage directly, for 
a given set of constraints. Unfortunately, because of the difficulties 
encountered in putting the GASP into an operational status, the second 
phase of the approach has not been completed at the time of this report. 
For this reason, this chapter is concerned primarily with the first 
phase - the evaluation and development of fuselage design and analysis 
methods. 

To facilitate the fuselage design procedures, the fuselage was 
considered to be composed of three distinct sections - the nose cone, 
the cabin and the tail cone. This is illustrated in Figure 2.1. Design 
methods were considered for each section individually and then were 
integrated to provide a method to determine the overall fuselage 
shape and size. The methods considered for each section are discussed 
in the following sections of this chapter. Section 2.1 compares 
different methods for sizing the cabin given a specified number of 
passengers. Section 2.2 describes a method for sizing the nose cone 
and crew compartment. Section 2.3 discusses a method for determining 
the fuselage/empennage configuration to meet specified static stability 
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criteria. Section 2.4 presents the results of a study to determine 
baggage compartment requirements for commuter airplanes and the 
alternatives for satisfying those requirements. Finally a method for 
sizing the overall fuselage is discussed in Section 2.5. 



Figure 2.1 Definition of Fuselage Sections 


2.1 Cabin Arrangement 

The fuselage cabin design is primarily dependent upon the 
number and arrangement of passengers. However, possible use of the 
airplane for cargo payloads, the location of the wing carry-through 
structure, and location of main landing— gear storage can also 
enter into the cabin design. Although the main emphasis 
of this research was to determine and evaluate methods to size 
the cabin for utility (i.e. passengers) constraints, a conscious 
attempt was made to keep other factors in mind. 

Before looking at existing sizing methods, or developing any 
new ones, a survey was conducted to examine the dimensions of existing 
commuter airplane cabins. The results of this survey are presented 
in Table 2.1, In Table 2.1 the seating configuration is indicated 
in an X/Y format where the values of X and Y indicate the number of 
seats abreast on either side of a single cabin aisle. The sole 
exception to this is the Britten-Norman Trislander which has no aisle. 
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Table 2.1 Cabin Dimensions of Existing Coimnuter Aircraft 


No 

Aircraft 

Seating 

Config. 

Passenger 

Seats 

Cabin* 

Width 

(in.) 

Cabin* 

Height 

(in.) 

Cabin** 

Length 

(in.) 

1 

Beech 99 

1/1 

15 

55 

57 

276 

2 

Cessna 402 

1/1 

7 

56 

51 

174 

3 

Britten-Norman Islander 

1/1 

10 

43 

50 

156 

4 

Swearingen Metro 

1/1 

19 

62 

57 

306 

5 

Queen Air 

1/1 

7 

52 

57 

108 

6 

GAF Nomad 

1/1 

15 

51 

62 

331 

7 

Piper Navajo Chieftain 

1/1 

8 

51 

52 

- 

8 

Saunders ST -27 

1/1 

22 

54 

69 

336 

9 

Beriev Be-30 

1/1 

16 

60 

60 

264 

10 

Pilatus Twin Porter 

1/1 

8 

46 

50 

120 

11 

Mitsubishi MU-2L 

1/1 

11 

59 

51 

- 

12 

Twin Otter 

1/2 

20 

63 

59 

222 

13 

Short Skyliner SC-7 

1/2 

22 

78 

78 

223 

14 

Falcon 30 

1/2 

30 

96 

73 

445 

15 

Short SD3-30 

1/2 

30 

78 

78 

372 

16 

YAK- 40 

1/2 

33 

80 

62 

336 

17 

Nord 262 

1/2 

29 

84 

70 

300 

18 

Casa C.212 

1/2 

18 

82 

70 

192 

19 

Handley Page Jetstream 

1/2 

1 18 

73 

71 

288 

20 

Piper PA35 

1/2 

16 

76 

70 

204 

21 

Omnipol L-410 

1/2 

14 

69 

65 

216 

22 

Fokker-VFW- 614 

2/2 

44 

105 

78 

444 

23 

Antonov AN-24 

2/2 

50 

109 

75 

444 

24 

Breguet 941c 

2/2 

56 

102 

88 

432 

25 

deHavilland DHC-5 

2/2 

53 

105 

82 

480 

26 

deHavilland DHC-7 

2/2 

48 

103 

77 

444 


* Inside Dimensions 

** Approxiraate-varies with configuration 
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No Aircraft 


Seating Passenger Cabin* Cabin* Cabin** 

Config. Seats Width Height Length 

(in.) (in.) (in.) 


27 

Fokker F.27 

2/2 

52 

88 

79 

- 

28 

GAC-100 

2/2 

32 

98 

75 

372 

29 

Handley Page 
Herald 

2/2 

54 

109 

76 

_ 

30 

Canadair CL-215 

2/2 

30 

94 

75 

360 

31 

Hawker Siddeley 
HS. 748 

2/2 

52 

97 

75 

558 

32 

NAMC YS-11 

2/2 

60 

106 

78 

528 

33 

lAI Arava 

2/2 

18 

75 

68 

144 


* Inside Dimensions 

** Approximate-varies with configuration 


Wiere the data were available seat, aisle and seat pitch dimensions 
were also examined, in addition to the inside cabin width and height 
dimensions. Figures 2.2 and 2.3 define the cabin and cabin seating 
dimensions that will be referred to throughout this section. 




Figure 2.3 Definition of Seating Arrangement Dimensions 
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The cabin sizing problem is one of trading-off passenger comfort 
and seat-miles for drag and weight. Passenger comfort may be related 
to the spacing between seats. As spacing is increased, so is pas- 
senger comfort. However, as comfort, and for that matter the number 
of total passenger seats is increased so are the drag and the weight 
of the overall airplane, which tends to reduce cost-effectiveness. 
Although this section only deals with sizing the cabin for specified 
comfort levels and seating arrangements, the trade-offs will be 
considered to some extent in Chapter 4, 

Comfort level is difficult to define. In lieu of trying to 
survey passengers as to what might be considered as 'comfortable' or 
not, and then having to extract the size-related factors from noise, 
vibration, ride quality, etc., it was decided to take what data were 
available for cabin seating arrangements and attempt to define 
general comfort levels. A letter from the project technical monitor 
at NASA-Ames, Tom Galloway (Reference 4), defined three comfort 
levels as shown in Table 2.2. For the most part the data available 
supported these values. The one exception was in the case of aisle 
width. A number of aircraft were found to have aisle widths lower 
than 16 inches. For this reason the comfort levels that will be 
used in this report are as presented in Table 2.3. 


Dimension 

Single Seat 
Width 

Seat Pitch 

Aisle Width 


Table 2.2 Seating Arrangement 
Comfort Level 
(From Reference 4) 


Minimum 

18" 


Adequate 

20 " 


Maximum 

22 " 


28" 30" 32" 
16" 18" 20" 
64" 70" 


Standup Headroom 


76" 


Table 2.3 Revised Seating Arrangement 
Comfort Levels 


Dimension 

Minimum 

Adequate 

Maximum 

Seat Width 

18" 

20" 

22" 

Seat Pitch 

28" 

30" 

32" 

Aisle Width 

12" 

18" 

20" 

Standup Headroom 

64" 

70" 

76" 

2.1,1 De terminat ion 

of Cabin 

Cross-Section 



As was stated in Chapter 1, the GASP was to be the medium of 

application for all design methods that were developed under this 

project. The GASP presently assumes that the cabin may be represented 

by a right circular cylinder of constant cross-section. For the 

purpose of this project the constant cross-section assumption was 

maintained. However, the cross-section was not limited to circular 

shapes. The smaller commuter airplanes frequently have cabin 

cross-sectional shapes which are far from being circular. 

For this reason, circular, elliptical, and what will be referred to 

as rounded-rectangular shapes were considered. This report will 

deal with the circular and round-rectangular shapes only. The 

elliptically shaped cross-sections are inherently covered by the 

methods that will be discussed here. The geometry for the rounded- 

rectangular cross-sectional shapes is presented in Figure 2.4. Note 

that the rounded-rectangular cross-section becomes an ellipse (or 

circle) for r^ and r, =1. 

tc be 



Figure 2.4 Geometric Definition of 

Rounded- Rectangular Cross-sections 
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For both circular and rounded-rectangular cross-sections the 
values for the inside cabin height and width were expected to be 
largely dependent on the seating arrangement selected and the comfort 
level. To help in finding the best methods for determining these two 
dimensions the data of Table 2,1 were plotted to correlate total 
passenger seats and seating arrangement with inside cabin width and 
height. These plots are presented as Figures 2.5 and 2.6. 

From Figures 2.5 and 2.6 note first that for the range of airplanes 
considered there seemed to be more or less definite passenger capacities 
for which each seating configuration was used. Table 2.4 presents 
the apparent passenger capacity ranges for each seating configuration. 
Note also from Figure 2.5 that the cabin width ranges for each 
seating configuration are well defined. These are also presented 
in Table 2.4. The cabin heights for each seating configuration are 
not so well defined in Figure 2.6. This would seem to indicate that 
cabin height is not as dependent upon the seating arrangement as was 
originally expected. This will be discussed further at a later point 
in this section. 


Table 2.4 Passenger Capacity and Cabin 
Width Ranges for Seating Configurations 
of Existing Commuter Aircraft 


Seating 

Config. 

1/1 

1/2 

2/2 


Passenger Capacity 
Range 

0-22 
15 - 33 
30 - 60 


Cabin Width 
Range (in.) 

43 - 62 
62 - 96 
88 - 109 


It is interesting at this point to compare the cabin width 
ranges of Table 2.4 with the comfort levels of Table 2.3. Using 
values resulting from the sum of the seat widths and aisle width as 
as estimated cabin width Table 2.5 compares the values for minimum 
and maximum comfort levels with the observed cabin width ranges. 
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Inside Cabin Width 










Table 2.5 would seem to indicate that the comfort level criteria are 
reasonable. The minimum comfort level does indicate larger cabin 
widths for the 1/1 and 1/2 conf igiircticns , i.n a number of 

cases aisle width values lower than 12 inches were observed. Also 


it should 

be noted that the 43 

inch 

cabin width was 

the Britten- 

Norman Islander which has no :i 

Is Id ♦ 




Table 2.5 Comparisons of Existing 



Cabin Widths 

with 

Comfort Level 



Cabin Widths 


Seating 

Cabin Width 


Minimum 

Maximum 

Config. 

Range 


Comfort Level 

Comfort Level 


(Table 2.4) 


Cabin Width 

Cabin Width 


(in.) 


(in.) 

(in.) 

1/1 

43 - 62 


48 

64 

1/2 

62 - 96 


66 

86 

2/2 

88 - 109 


84 

108 


In an attempt to derive relationships for cabin width and cabin 
height as a function of passenger capacity and/or seating arrangement 
least-squares linear and logarithmic curve fits were applied to both 
Figures 2.5 and 2.6. The results of these curve fit applications are 
presented as Figures 2.7 through 2.9. Note that the cabin width 
data are only presented in the "thumbprint” forui. This has been done 
purely to add clarity to the figures. All of the curve fit calculations 
are based on the data of Table 2.1. Note in Figure 2.7 that the 1/1 
and 2/2 configuration lines are almost parallel and that the 1/2 con- 
figuration line seems to represent a transition between the two. A 
linear regression analysis was performed for each seating configuration 
for the cabin height case but the resulting relationship reinforced 
the conclusion made earlier that cabin height is not very dependent 
upon the seating arrangement. Both of the logarithmic curve fits 
(Figures 2.8 and 2.9) show that cabin width and cabin height may be 
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Figure 2.7 Linear Correlation of Inside Cabin Width with 
Number of Passengers 
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de Cabin Width (in.) 



Eigure 2.8 Logarithmic Correlation of Inside Cabin Width 
with Number of Passengers 
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Figure 2,9 Logarithmic Correlation of Inside Cabin Height 
with Number of Passengers 








relatsd to total passenger capacity. 

As both cabin width and cabin height may be related to total 
passenger capacity but only cabin width may be easily related to the 
number of seats abreast, then it was thought that perhaps cabin 
height could be related to the cabin's other primary dimension, 
length. It seems reasonable, at least from a comfort point of view, 
that the longer the cabin, the less passengers will like having to 
put up with a low ceiling aisle. Figure 2.10 relates cabin height 
with cabin length. Figure 2.10 does seem to indicate that there is 
some relationship but the correlation is not as good as was expected. 

No further research was pursued along this line. 

2 . 1 . 1 . 1 Existing Cross-Section Sizing Methods 

Only three existing cross-section sizing methods were available 
for evaluation - the Boeing Vertol Method (Reference 5) , the McDonnell 
Douglas method (Reference 6), and the GASP method. Each of these 
methods assumes a circular cross-section. 

The Boeing Vertol method is based on the assumption that a 
circular cabin cross-section can be designed for a particular seating 
arrangement that will clear certain body-seated control points. The 
method was developed in the following manner. Five control points 
were chosen to constrain the inner cabin wall about a seated passenger 
as shotm in Figure 2.11. Using these control points layouts for 
different configuration layouts were made. An empirical relationship was 
developed based on the number of seats abreast, seat width, number of 
aisles, and aisle width using the data derived from the layouts. 

This relationship is presented in Figure 2.12. The outside body 
radius was then computed as a constant percentage of the inner radius. 

A comparison of actual outside body radii and computed outside body 
radii is presented in Figure 2.13. 

Although reference will be made to the McDonnell Douglas method, 
no actual sizing method was available. Instead the results of a study 
of the operational requirements for medium density air transportation 
were available. These results are presented in References 6, 7 and 8. 








Figure 2.12 Empirical Relationship for Cabin Radius by 
the Boeing Vertol Method (Reproduced from 
Reference 5) 


BTAKR 
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Predicted Outside Body Radius (in.) 




Figure 2.14 is a reproduction from Reference 8 of the cross-section 
developed for the baseline feederliner of the McDonnell Douglas 
study. The cross-section was designed not only with passenger 
seating in mind, but also to accept palletized cargo. The only 
control points specified to constrain the shape for passenger seating 
were at the tangent to a seven inch circle drawn about the center of 
the head for a 95 percentile seated man and at the top-center of the 
cabin. A further control point was defined at the base of the outer 
seat leg, six inches from the seat center line. With the seat dimensions 
and aisle height (i.e. cabin inside height), these points could be 
located and the cross-section shape determined. 

Having specified three control points on the inner wall of the 
cabin, and defining these points relative to the cabin center point, 
the cabin cross-section is over-constrained. Only two constraints are 
required, but by having three it is possible to choose which constraints 
are critical in the case being considered. Equations 2.1, 2.2 and 2.3 
define the inside cabin radius from each control point in terms of the 
seating configuration and the cabin center point: 

A. Ceiling Constraint (Control Point A) ; 

r = h - z (2.1) 

c c 

B. Passenger Headroom Constraint (Control Point B) : 


r^ =-y/[l/2((Ng ' 

C. Floor Constraint (Control Point C) : 

=-\/tl/2((Ng - l)Wg + N^w^) + 6]2 + (2.3) 

where, in the above equations, 

N = number of seats abreast 
s 

w « single seat width (in.) 
s 

N, = number of aisles 
A 
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= aisle width (in,) 

z = vertical location of the cabin center point with respect ■ 
to the floor (in,) 

h^ = inside cabin height from the floor (in*) 

It is assumed that the only unknowns in equations 2.1j 2.2, and 
2.3 are r and z. The value of cabin height will be assumed to be 
a given comfort level value. To find the cabin radius any pair of 
constraints may be chosen. By then solving the resulting pair of 
equations for the cabin is sized for those particular constraints. 

The third equation may then be solved to determine if the cabin 
radius is adequate for that constraint. 

In most cases, constraints B and C will probably be critical. 

By solving equations 2.2 and 2.3 simultaneously for these constraints 
can be met* However, in this case the cabin height must also become 
a variable. After determining the cabin radius and the center point 
locations, the cabin height may be checked. If the cabin height is 
inadequate the cabin may be sized for new constraints. 

Equations 2,4, 2.5 and 2.6 provide the solutions for r^ for 
each constraint pair. 


Constraints A and B : 

+ (48 - - 49 

r = 

2(55 - h^) 


Constraints A and C t 



Constraints B and C: 


-B ± - 4AC 


(2.4) 


(2.5) 


2A 


( 2 . 6 ) 


where 



For all of the above equations; 



1 

2 


W3 - l)w^ 


Va 



Figure 2.14 McDonnell Douglas 30-Passenger 
Feederliner_ Cross-Section 
'(Reproduced From Reference 8) 
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Table 2.6 Cabin Widths* for the 


Minimum Comfort Level 


Seating 

Config, 

Boeing 

Vertol 

Method 

McDonnell 

Douglas 

Method 

GASP 

Method 

(Outside) 

1/1 

69.65 

67.04 

60.00 

1/2 

80.97 

79.26 

78.00 

2/2 

93.49 

93.53 

96.00 


Table 2.7 

Cabin Widths* for the 




Adequate Comfort Level 


Seating 

Config. 

Boeing 

Vertol 

Method 

McDonnell 

Douglas 

Method 

GASP 

Method 

(Outside) 

1/1 

75.79 

72.14 

70.00 

1/2 

89.19 

87.00 

90.00 

2/2 

104.06 

103.75 

110.00 


Table 2.8 

Cabin Widths* for the 




Maximum Comfort Level 


Seating 

Config. 

Boeing 

Vertol 

Method 

McDonnell 

Douglas 

Method 

GASP 

Method 

(Outside) 

1/1 

79.66 

74.90 

76.00 

1/2 

94.96 

91.88 

98.00 

2/2 

112.05 

110.77 

120.00 


* Cabin Width in inches 
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The GASP method only sizes for the outside dimension of the 
cabin. As with the other methods a circular cross-section is 
assumed with a diameter determined as in equation 2,7. The constant 
twelve inches are intended to account for clearances about the seats 
and for the structure. 

D = N • w + N. • w. + 12 (inches) (2.7) 

s s A A 

where; D - outside cabin diameter (inches) 

c 

o 

N ~ number of seats abreast 
s 

N, = number of aisles 
A 

w = seat width (inches) 
s 

w. “ aisle width (inches) 

A 

2. 1.1. 2 Evaluation of Cross-Section Sizing Methods 

To evaluate the application of the cabin sizing methods 
described above to the design of short haul and commuter airplanes, 
it was decided to compare cabin widths computed by these methods for 
each comfort level with the data of Table 2.1 and Figure 2.5., 

As the cross-section is assumed to be circular at this point the 
inside cabin width, w^, may be assumed to be equal to the inside 
cabin diameter. Tables 2.6, 2.7, and 2.8 present the results for 
the minimum, adequate, and maximum comfort levels respectively. 

Figures 2.15, 2.16, and 2.17 plot the results of each of these 
tables respectively on the seating arrangement thumbprints of 
Figure 2. 5, for each sizing method and each seating arrangement. 

From Figures 2.15 through 2.17 it is apparent that the methods 
generally size for cabins larger than those of existing aircraft. 

This is especially true for the 1/1 configurations. In this case the 
GASP method seems to do the best job when it is taken into consider- 
ation that the GASP sizings are for outside widths. The Boeing Vertol 
and McDonnell Douglas methods do fairly well for the 1/2 and 2/2 
seating configurations however. It is also noted that the McDonnell 
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Inside Cabin Width Wc (in.) 


1) Cabin Sizing Method 

Boeing (inside) 

Douglas (inside) 

GASP (Outside) 

2) Thumbprints Indicate 
Actual Aircraft from 
Figure 2. 5 


Passenger Seats PAX 


Figure 2.15 Comparison of Computed Minimum Comfort Level 
Cabin Widths with Existing Aircraft Data 




Inside Cabin Width Wc (in.) 


1) Cabin Sizing Method 

Boeing (inside) 

Douglas (inside) 

GASP (Outside) 

2) Thumbprints Indicate 
Actual Aircraft from 
Figure 2.5 


20 30 

Passenger Seating ~ PAX 


Figure 2,16 Comparison of Computed Adequate Comfort Level 
Cabin Widths with Existing Aircraft Data 





Passenger 


Figure 2.17 Comparison of Cor 
Cabin Widths witl 





Douglas and Boeing Vertol methods represent very similar results 
in all cases. 

To determine why the Boeing Vertol and McDonnell Douglas methods 
over-sized cabin width for the smaller seating configurations (i.e. 

2 and 3 seats abreast), the extreme case was considered - the 1/1, 
minimum comfort configuration. The cabin was sized using both 
methods in Table 2.6. Figure 2.18 presents a drawing showing the 
resulting cross-sections and seating arrangement. Note that for the 
Boeing Vertol cross-section control points B, C, and D are not 
needed and that points A arid S become, the critical constraints* This 
results in a large amount of wasted space between the outer armrest 
and the cabin wall* The same occurs for the McDonnell Douglas 
cross-section. To provide the seated headroom constraints for these 
smaller seating configurations the width must be oversized- 

Figure 2.19 shows the result of the GASP method which produced 
a cabin width which was comparable with existing aircraft. The 
seat used in both Figures 2*18 and 2*19 was a conceptual economy- 
class seat sized to the minimum comfort level dimensions of Table 
2-3. The seated headroom constraint shorn in Figure 2.19 is the same 
as that for the McDonnell Douglas method. The seven inch clearance 
radius is drawn about the location of the center of a 95 percentile 
seated man (48 inches above the floor at the seat center line) . As 
the GASP method only sizes for the outside cabin dimension, the 
Boeing Vertol method for estimating structure was used to determine the 
inside cabin radius (See Figure 2.13). Figure 2.19 indicates that 
with these constraints it would be extremely difficult to fit a 95 
percentile man in the cabin. Although seat dimensions for smaller 
seats were not readily available, it is assumed that with such a 
seat (i.e. lower seat height, greater nominal recline angle, etc.), 
such a cabin would become cramped but practical. Indeed, several 
aircraft such as the Gates-Learjets already do this. For a pressurized 
airplane where a circular cross-section means weight savings, this 
Chen would be a viable cabin sizing method. 

In the case of unpressurized commuters, a better solution is 
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2) Comfort Level: IVlinimum 

3) Control Points 

® McDonnell Douglas 
A Boeing Vertol 


Figure 2,18 Boeing .and Douglas C oin Cross-sections for 
Minimum Comfort Level 





available. The rounded-rectangular cross-section of Figure 2.4 
can result in a more comfortable cabin without giving up much in the 
way of weight or drag. In this case, however it may no longer be 
assumed that outside body width equals body height. Therefore each 
must be sized. The values of the round-off radii are for the ma‘?t 
part a matter of preference. 

Equation 2.8 may be used to determine the inside cabin width. 

In equation 2,8 it is assumed that the window-side armrest of the 
seat is placed directly up against the cabin wall. 


w 


c 


N w 
s s 


+ N.w. 
A A 


( 2 . 8 ) 


To determine the fuselage inner height (i,e, total outside 
body height minus structure) , first a control point is established 
as in control point C for the McDonnell Douglas method. It is assumed 
that to minimize fuselage inner height, the floor is to be placed 
as low as this control point will allow (as shown in Figure 2.20). 
Using the bottom round-off radius fraction, should be 

possible to determine the necessary fuselage height. Let the height 
below the floor be represented as h^. The inside fuselage height is 
then as represented in equation 2.9. 



Control Point 


Figure 2.20 Location of Floor Control Point 

for Rounded-Rectangular Cross-Sections 




page buane not 


h = h + h. 
c f 
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The round-off shape is elliptical. Therefore at the control 
point : 



where , 

’"cSc 
y ^ 



z 



- hf 


( 2 . 10 ) 


Substituting equation 2.9 and solving for h results in: 



( 2 . 11 ) 


Table 2.9, presents inside cabin widths and inner fuselage 
height for each comfort level and seating configuration. In Table 
2.9 round-off radii ofr, = -2, .5, .8, and 1.0 are used. From 

DC 

Table 2.9 it is apparent that for large round-off radii, r^^, the 
minimum inside body height, h^, required to meet comfort level 
constraints can become overly large. Therefore it would seem that 
if cabin wetted area is to be minimized, a smaller r^^ would be 
preferable. 

To check that adequate room for a passenger is indeed provided , 

the minimum comfort level 1/1 seating configuration was again used. 

Figure 2.21 provides the resulting cabin cross-sections for upper 

round-off radii values of r^ = .2, .5, .8, and 1.0. In all cases 

‘'c 

the lower round-off radius is r, =0,2. Note that this gives very 

be 
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Table 2^9 Cabin Dimensions for 
Rounded-Rectangular Cross-Section 


Comfort 

Level 

(Table 2,3) 

Seat 

Config. 

Cabin Width, w 
(in.) 


Inside Body Height, 


r, = .2 
be 

'be ■ -5 

r, = .8 
be 

'be = 

Minimum 

1/1 

48 

64.6 

73.3 

94.8 

128.0 

Minimum 

1/2 

66 

65.3 

75.2 

97.1 

128.0 

Minimum 

2/2 

84 

65.8 

76.5 

98.6 

128.0 

Adequate 

1/1 

58 

70.4 

79.4 

102.7 

140.0 

Adequate 

1/2 

78 

71.1 

81.5 

105.4 

140.0 

Adequate 

2/2 

98 

71.7 

82.9 

107.0 

140.0 

Maximum 

1/1 

64 

76.2 

85.2 

110.1 

152.0 

Maximum 

1/2 

86 

76.9 

87.6 

113.3 

152.0 

Maximum 

2/2 

108 

77.5 

89.2 

115.2 

152.0 
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reasonable results for < 0,8. Although other examples are 
not presented here for the sake of brevity, in general it was found 
that for r^^ <0-8 all seated headroom constraints could be met. 

It should be noted again that this method sizes the inside 
cross-section only for the minimum dimensions, given the comfort 
level constraints. It might be desirable, for instance, to increase 
the under- the floor height to allow for baggage, wing carry-through 
structure or landing gear retraction. These considerations are not 
inherently taken into account by this method# 


Methods have been discussed to size both circular and rounded- 
rectangular cross-sections. It was found that all of the circular 
cross-section methods provided results which were reasonable for 
configurations with three and four seats abreast, but that appeared 
to be either too small (GASP) or too large (Boeing and Douglas) 
for configurations with only two seats abreast. Of these methods 
only the GASP method did not directly account for adequate clearances 
for a seated passenger. Of the other two methods the McDonnell 
Douglas method seemed to provide the best results for the two-seats- 
abreast case. Also the McDonnell Douglas method is slightly easier 
to check and to alter should different seat sizes be chosen. For 
these reasons it was decided to adopt the McDonnell Douglas method 
to incorporate into GASP for circular sections., The rounded-rectangular 
method also appears to provide reasonable results, and in this case, 
particularly for the smaller seating configurations. This method has 
the distinct advantage of providing a way of sizing short haul 
aircraft for cargo payloads, also. The round-off radii may be input 
as desired, although it was found that for round-off radii greater 
than 0.8 a circular sizing method would be preferable. Therefore 
this method will also be used to provide added versatility to the 
cabin sizing process. 

Up until this point only the inner cabin dimensions have been 
dealt with. Only one method was found to account for structure. 
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upholstery, etc., in general. The Boeing method determined the inside 
cabin diameter to be 94% of the outside diameter. Although this is 
somewhat conservative for smaller unpressurized airplanes, this method 
will be used to determine the outside dimensions. 


2.1.2 Determination of Cabin Length 

At this point the cabin length will be determined only as that 
length required for passenger seating. It will be assumed that for 
short haul/coramuter airplanes of size and range being considered that 
cabin attendants and lavatory facilities will not be necessary. 
Baggage compartment requirements will be considered in Section 2.4. 

To determine the cabin length for passenger seating equation 
2.12 will be used: 


a 

“pax 


= N 


ROWS 


a 

s , 


where , 

N „„„ = number of rows of seats 
ROWS 

a = seat pitch in inches, 
s 

2.2 Nose Cone Configuration 


( 2 . 12 ) 


The primary consideration in designing the nose cone and wind- 
shield geometry is the cockpit arrangement and visibility. The cockpit 
constraints must be met if the crew are to be able to fly the airplane. 
Other considerations however include the drag increment resulting from 
the windshield, nose gear retraction, and some stability characteristics 
The main emphasis in this section will be placed upon the design of the 
nose for adequate cockpit volume and visibility. 

The GASP default* size for the cockpit has a flight deck 4.44 ft. 
long. Based on the data available in References 9 and 10 this 
might prove inadequate. Figure 2.22 presents the recommended 


*Note: A "default" quantity is a value stored in the computer 

routine which is used as the value of a particular 
variable when no value is input by the user. 
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Figure 2.22 Recommended Flight Deck Dimensions 
For Transport Aircraft 
(Reproduced from Reference 9) 






Using the dimensions of Figure 2 -22 it was possible to define 
a more or less standard shape for the crewbox. The dimensions for 
this crewbox are presented in Figure 2.25. Note from Figure 2.25 
that the width of the crewbox is defined as in equation 2.13. The 
minimum distance between crew seat center lines for the class of 
aircraft in question should be about thirty inches, according to 
Reference 9. Also it is assumed that, as a minimum^ the crew seats 
are of the same width as the passenger seats. 

w , ” w +30 (inches) (2.13) 

cb s 

The length of the nose shell, i.e. the long thin cone of 
Figure 2.23 is difficult to specify in any generalized form. It 
was found that the best method is to specify the ratio of the 
lengths of the two cones. In choosing a length ratio certain 
factors should be kept in mind. Some of these are: 

• Forward and downward visibility on approach 

• Nose gear location and retraction method 

• Possible installation of radar 

• Desirability of a nose baggage compartment 

• Longitudinal and directional static stability. 

Although attempts were made to develop a method to size the 
nose shell based on these factors no successful method was derived. 
Instead, the length ratio approach seems to provide adequate results. 
It was found that for commuter aircraft the length ratio was normally 
between 1.5 and 2.5. 

2.3 Tail Cone/Empennage Configuration 

In line with the objective of minimizing the fuselage/empennage 
drag and weight while maintaining stability constraints, the Roskam/ 
Fillman method was chosen to be used to size the tail cone and 
empennage. The Roskam/Fillman method as originally published in 
Reference 11 represented an approach to optimizing fuselage and 
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empennage size with respect to zero-lift drag while maintaining 
certain stability and control characteristics. The zero-lift 
drag of the fuselage is a direct function of the fuselage fineness 
ratio, If it is assumed that the size and shape of the aircraft 

nose and cabin are dictated by the utility constraints then the 
fineness ratio may be varied by altering the ‘tail cone’ length, 

If certain stability and control characteristics are specified 
as constants, then the optimization process becomes a trade-off 
between fuselage and empennage drag* Fuselage zero-lift drag may be 
approximated by equation 2.14 (from Reference 12)* 



= C [1 + 

B 


60 

(ag/D)3 


+ .0025 


(Swet) bod y 
^Wing 


(2.14) 


The zero-lift drag of either the vertical or horizontal tail 
may be estimated by equation 2,15 (from Reference 12). 


\ ■ '' " (2.15) 

° Vor H ^Wing 

As the fineness ratio, l-o/B, is increased by lengthening the 

. u 

tail cone the following occur; 

• Increasing ii^/D will increase the fuselage Reynolds number 

n 

decreasing C- 


* Increasing will decrease the value in the brackets of 

equation (2*14) 

* Increasing will increase ^ , 

D wet Body 

* Increasing 5.-/D will reduce the empennage wetted area 

D 

requirements, for constant stability levels. 


One simplifying assumption made in the method as originally 
published in Reference 11 was that the aircraft aerodynamic center 
location could be considered as a constant as the tail cone was 
lengthened. The overall intent was to trade off a reduction in 
empennage drag for an increase in fuselage drag in hopes of finding 


the minimum drag configuration. 

To use the Roskam/Fillman method to size the fuselage and 
empennage for minimum weight and drag a number of changes became 
necessary. Weight estimation of the fuselage and empennage was 
added as well as a c,g, location method. With these additions 
static stability constraints could be specified as a constant C 

a 

for the longitudinal case and a constant C for the directional 

case. As with the original method, it was assiimed that the airplane 
geometry foirward of the tail cone would remain unchanged. 

To test the method a computer program was written in FORTRAN IV 
to facilitate the methods application to both existing and conceptual 
aircraft. This section will describe the program methods and 
discuss the results it produced. 

2,3.1 Roskam/Fillman Program Description 

A flowchart of the program is presented in Figure 2.26. A 
program listing and an explanation of how the input data are prepared 
are presented in Appendix B of this report. Explanations of the 
methods used in each of the main subroutines are provided in the 
following paragraphs. Note that on Figure 2,26 Section numbers 
refer to the section in which a discussion of each method is pre- 
sented* 


2. 3.1.1 Computation of Wetted Areas 


The wetted areas for the fuselage and empennage are computed 
in the following manner in the subroutine SI^T. 

The fuselage wetted area is computed in three components: the 
cabin wetted area, the nose wetted area, and the tail cone wetted 
area. The cabin is assumed to be a right circular cylinder. The 
cabin wetted area is expressed as in equation 2.16* 


^Wet 


Cabin 


irDJl 

u 


(2.16) 
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Figure 2.26 Flowchart of Roskam/Fillman 
Method Program 
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The wetted areas of the nose cone and tail cone are computed 
using Torenbeek’s elliptical cone methods of Reference 9. This 
approach is documented in Appendix A of this reports By using Figures 
A-1 and A-3 of Appendix A to determine the shape parameters, 
for the nose and tail cones, it is possible to determine a wetted area 
correction factor, k^, from Figure A-2 of Appendix A. The wetted 
areas for the nose and tail may then be expressed as in equations 
2*17 and 2*18; 

NOSE 

^ ^NOSE NOSE 

TAIL 

(SWET)^^j^ ^ ^WaIL 

The wetted area of the fuselage is then expressed as: 

(She«3 . + (=WET)^„3j + 

The horizontal and vertical tail wetted areas were computed with 
the aid of a surface area correction factor, ^ from Coming 
(Reference 13) . This surface area correction factor is a function 
of the airfoil thickness ratio, t/c, as shown in Figure 2.27. 
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factor 
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Figure 2.27 


Determination of the Surface Area Correction 
Factor, (Reproduced from Reference 13) 


2.44 



Using Figure 2.27 the wetted areas of the horizontal and vertical 
tails may be determined as: 


or V 

H or V 

^'/here , 

K = 0.52C-) + 1.987 
P c 
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2. 3. 1.2 Calculation of Zero-Lift Fuselage/Empennage Drag 


The zero-lift drag of the fuselage and empennage are computed 
in the subroutine CDO using equations 2.14 and 2.15 as stated 
earlier. 



(2.15) 


The mean skin-friction coefficients, C^, are a function of 
Reynolds number, R^, and Mach as shown in Figure 2.28. The factors 
L and g are as defined in Reference 11. 

2. 3. 1.3 Calculation of Fuselage and Empennage Weight 

The structural weights of the fuselage and empennage are mainly 
dependent on their areas and the maximum speed of the aircraft. The 
following equations were taken from Reference 9. The equation for 
fuselage shell weight estimation is: 

"f - '(=g)'' 
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I’Jhere =* .021 for in LBS 
W = Fuselage Weight (lbs) 

r 

Vjj = Design dive speed (kts) 

S = Gross shell area (ft^) 

G 

= Maximum fuselage width (ft,) 

= Maximum fuselage height (ft.) 

= Distance from wing quarter chord to horizontal tail 
quarter chord (ft.) 

- 1.00 for unpressurized fuselage 

= 1*08 for pressurized fuselage 

= 1.00 for wing mounted engines 

- 1.04 for fuselage mounted engines 

= 1.00 for main gear attached to wing 

= 1.07 for main gear attached to fuselage 

= 1.00 for main gear bay in the fuselage 
= 0.9b for no main gear bay in the fuselage 

As can be seen from equation 2.21, the weight is dependent on the 
square root of the dive speed and the wetted area to the 1.2 power! 

The empennage weight is dependent on the load factor and the 
square of the wetted area for dive speeds less than or equal to 
250 kts. 


^EMP [^^ULT 

where = 0.04 for in KTS 


^EMP 

^EMP 


'ULT 


EMP 


= Ultimate airplane load factor 
= Empennage Planform area 


( 2 . 22 ) 
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For dive speeds greater than 250 kts. the empennage weight is a 

direct function of the planform area and the terms and F„ which 

ri V 

are defined below. 


= ‘Si • <2-23) 

”vT ■ S ■ \ 


where 

"h 


Dive speed in kts- 

Horizontal Tail Planform area 

Vertical Tail Planform area 

1-0 for fixed stabilizer 

1.1 for variable incidence stabilizer 

1.0 for fuselage mounted horizontal tails 

1.0 + .15 for T^Tails 


where b = Height of Horizontal Tall above Fin Root 
H 

= Vertical Tail span 



1000 -^cos 



1000 yjcos 


(2.25) 


(2.26) 


where A = sweep in degrees at the maximum airfoil thickness 


The values of the parameters dependent on F and F can be determined 

n V 

from Figure 2,29. These equations were programmed and are used to 
calculate the shell weights used in the Roskara/Fillman program in the 
subroutines FUSWGT and EMPWGT. 
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Figure 2.29 Empennage Weight Function 

(Reproduced from Reference 9) 

One of the input parameters for the Roskam/Fillman program is 
the baseline airplane gross weight. Using this gross weight and 
the baseline shell weights of the fuselage and empennage, it is 
possible to correct for change in gross weight due to the change in 
tail cone and empennage configuration. This is done in the following 
manner. Equation 2.21 shows that the fuselage weight is proportional 
to fuselage gross wetted area to the power 1.2. It is assumed that 
to approximate the shell weight of the tail cone equation 2.27 may 
be applied. 


W. 


TAIL 


(^WET) 

(®WET) 


TAIL 

BODY 




(2.27) 


This weight and the weight of the empennage are then subtracted 
from the baseline gross weight to provide an adjusted gross weight, 

. It is then assumed that C^gr^ ' remains constant and that on)/ 
the tail cone and empennage weights vary. 

2.3. 1.4 Location of the Center of Gravity 

The longitudinal c.g. location is required to size the empennage 
for stability constraints. The c.g. location subroutine, SBARCG, 
performs this function. The routine is not Intended to locate the 
c.g. accurately, but rather to estimate the shift in c.g. location 
due to incrementally lengthening the tail cone. The equation used 
to locate the c.g. was as expressed in equation 2.28, 
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Where, 


^^GROSS^^ 


+ £ 


eg 


(2.28) 


^GR = adjusted gross weight (lbs) (See Section 2. 3. 1.3) 




- baseline c.g. location relative to the nose (ft.) 


W. 


TAIL = tail cone shell weight (lbs) 


("eg)' 


W, 


tail = tail cone c.g. location relative to the nose (ft.) 

H = horizontal tail weight (lbs) 

= horizontal tail c.g. location relative to the nose (ft.) 

H 


W, 


V 


= vertical tail weight (lbs) 


(X j ~ vertical tail c.g, location relative to the nose (ft.) 

V ^/V 


(^A 

i 


gross! ~ total gross weight for i'th tail cone /empennage 
/i configuration (lbs) 


eg 


= c.g. correction found as the difference between 
estimated c.g. location for baseline and actual c.g. 
location for the baseline. 


The tail cone c.g. location X 


eg 


was located at the centroid 


TAIL 


of a body of revolution having a planform shape modelled by the equation: 




= 1 


(2.29) 


With a cone shape of this type, it may be shown that the centroid 
of the tail is located by equation 2.30. 


^AIL 



'l f 1 ^ 

! f - n* ^ 

/ 1 - 3m + 2m^N 


2 \^ra(n + Z) j 

(2n + 2)y 

(3n + 2) j 


[i f ^ Y 1 

T l-m ^ , 

/l - 3m + 2m^\ 



“^m(n +1) f ' 

^4m2(2n + 1) y ’ 

6m3(3n + 1) j 



(2.30) 


Therefore, the c.g. of the tail, relative to the nose is 
located as: 


X 


eg 


TAIL 






(2.31) 


The baseline c,g. locations for the horizontal and vertical tails 
were input parameters, although for preparation of the data, where more 
accurate data were lacking, it was assumed that the c.g, was at 60% 
of the mean geometric chord,. To correct for tail cone length 
changes equation 2.32 was used. 
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(2.32) 


where. 



= baseline c.g, location it-iative to the nose (ft,) 


= change in tail cone length relative to the baseline 

(ft.) 


As stated before, this routine was intended to provide only 
quick and dirty estimations for c.g. location based on Ihe limited 
input parameters that the program used. It was anticipated that if 
the method had been integrated into the GASP, that the GASP c.g. 
location method would provide better . results . 


2. 3.1. -5 . Empennage Sizing 

The horizontal and vertical tails were to be sized to meet 

specified static stability constraints. The horizontal tail was to 

be sized to provide a specific C value, and the vertical tall a 

a 


specific C value. In addition to this, it was decided to optimize 

the effectiveness of the horizontal and vertical tails for each tail 
cone length. These functions were performed for the horizontal and 
vertical tails in the subroutines STABAREA and VERTAREA respectively. 
The effectiveness of each surface was to be optimized by 
maximizing the product of the lift curve and moment arm for each 
surface. The variable was to be the sweep of the surface. For 
instance, by varying the sweep of the horizontal tail, both the 
horizontal tail moment arm, and the horizontal tail lift curve 
slope, , are affected- The tail may be said to be most effective 

where the product C« x ^i_ is a maximum. The value of JL, is 

L h n 

easily defined geometrically, and the value of C , from Reference 

“h 

14, may be expressed as' 
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A = aspect ratio 

3 = Prandtl-Glauert transformation, yjl - m 2 

K = ratio of the actual section lift curve slope to 2ir 
half-chord sweep 


The program assumes that the aspect ratios of the horizontal 

and vertical tails remain constant, allowing this procedure to be 

more or less independent of the sizing. 

The horizontal tail is to be sized for a constant C value. 

m 


From Reference 14: 



(2.34) 
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where all of the above values are total aircraft values. 


From Referer.ee 14; 



(2.35) 
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All of the variables on the right-hand side of equation 2.36 
are determined using the methods of Reference 14. In particular, 
however, the X is determined as: 




= X 
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Whereas the X of the wing, X , 
ac acy’ 

aerodynamic center shift due to the body 


is an input 




parameter, the 
computed by performing 


a Multhopp strip-integration as described in Reference 14. This is 
accomplished in the subroutine MULTOP. 

The total aircraft lift slope in equation 2,34 may be expressed as: 
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The program uses equations 2.36 and 2.38 to iterate to find a 

value for S., which will meet the C constraint. 

H m 

a 

The vertical tail size is determined in a similar manner. 

From Reference 14, the value for the C of the aircraft may be expressed 
as; 
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where, 

C = the body contribution to C (rad 

"Sb 


n„ 


= the wing contribution to C (rad 
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C = the vertical tail contribution to C (rad 
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It is conservative to assume that the wing contribution is 
negligible. Therefore, equation 2.39 becomes; 
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or, 
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= C - C 
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From Reference 14 the body contribution may be determined as; 

Jl 

\ ^ ^ (rad-1) (2.41j 

^B ^ 

where , 

Kjj = an empirical factor for body and body + wing effects 
determined from Figure 7.19 of Reference 14. 

= a Reynold’s number factor for the fuselage determined 
^ from Figure 7.20 of Reference 14. 

S = body side area (ft^) 

a 


= body length (ft) 
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The side body area is computed using the methods of Appendix A 
of this report. 

The vexcical tail contribution, C , may be expressed as shown 


in equation 2.42 (Reference 14) 
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For an airplane in the cruise configuration it is reasonable 
to assume that the angle of attack is small. Therefore: 
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From Reference 14: 
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(rad 1) (2.44) 
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V 

where , 

k = an empirical factor from Figure 7.3 of Reference 14. 

C = the vertical tail lift curve slope determined from 

Lt 

equation 2.33 using A = value of the effective 

aspect ratio, is determined by the methods of 

Reference 14. 

The value of the factor 1 4- r\ is assumed Le approximately 

CEp 


Substitution and -solving equation 2.40 for renders: 
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By iterating between the equation for C and equation 2.45, 

Li 

“v 

the value for S„ which meets the C constraint may be determined. 
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2.3.2 Results of the Roskam/Flllman Program 


The Roskam/Fillman method was tested using data for both actual 
and conceptual baseline aircraft. Although preliminary results for 
the Gntes-Learjet Model 24 seemed to verify the method, later results 
for both other existing aircraft, and for the conceptual aircraft in 
particular, were rather disappointing. For the sake of brevity 
only a generalized discussion of the results and the conclusions 
based thereon will be presented here. A more detailed discussion 
of the airplane configurations tested and the results is presented 
in Reference 15. 

When the tests were made, as was expected the vertical tail 
size required to maintain C decreased with increasing tail cone 

length. On the other hand, with many of the aircraft tested, as 

the tail cone length wa5? increased so did the horizontal tail size 

required to maintain C . This was not as expected. As a result 

a 

the increasing drag and weight of the horizontal tail and tail 
cone together over-rode the decreasing drag and weight of the 
vertical tail- 

After carefully evaluating the Roskam/Fillman algorithm the 
following conclusions as to the cause of the difficulties were 
arrived at. The method assumes that the wing may be kept at a 
constant location relative to the nose of the aircraft* Because 
of this and the fact that the aft portion of the fuselage has 
very little effect on the shift in aerodynamic center due to the 
body, the location of the wing + body aerodynamic center remains 
relatively constant* As the tail cone is lengthened the aircraft 
c.g. moves aft. It is probable that inherent inaccuracies in the 
method used to locate the c.g. pushed it even 
far the X' aft than should have been the case. It is also reasonable 
to assume that the lift slope of the total aircraft is not going 
to change drastically with changes in tail cone and horizontal 
tail size. Therefore by referring to equation 2.34 note that 
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as the c.g. is pulled aft by the tail cone, the aircraft aero- 
dynaaiic center is also pushed aft to maintain the static margin. 
Realizing this and the fact that the wing + body aerodynamic 
center is more or less fixed, note that the numerator of equation 
2.36 is going to become considerably larger. At the same time, 
the (X - X ) factor in the denominator may also be decreasing, 

compounding the effect. 

The program was written so that some of the computations 
made for each tail cone increment relied on some of the data from 
the previous tail cone increment* This was done to reduce the 
number of iterative cycles required. In particular, to locate 
the c.g. at a new tail cone length, the empennage sizes from 
the previous tail cone length were used. This was originally 
felt to be reasonable as long as the increments in tail cone 
length were small. Unfortunately, this also resulted in pushing 
the c.g. even farther aft in light of the oversized horizontal 
tail. Thus the problem was doubly compounded. 

From this evaluation of the algorithm employed the following 
conclusions and recommendation as to the usefulness of the method 
may be made: 

Before the method may be applied reliably to determine if 
an optimum tail cone length (from a weight and drag vie^^point) 
may be found: 

A) An accurate c.g. location routine is needed. Also 
if this is to be done, it would be wise to include 
a more accurate estimation of the component weights 
and how they might be altered by changing tail cone 
length. 

B) The re-balancing of the wing on the fuselage for each 
new tail cone cone configuration might be investigated 
to prevent the distance between the wing + body aero- 
dynamic center and the total aircraft aerodynamic center 
from becoming overly large. One disadvantage to the 
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re-balancing of the wing would be to reduce the effect 
of an increased horizontal tail moment arm, with the 
possible result that the overall objective of the 
method is cancelled. 

C) The algorithm needs to be revised such that the 

calculations for each new configuration are independent 
of all previous configuration data* This would reduce 
the probability of compounding errors. 

The most obvious means of accomplishing these revisions was 
to include the empennage sizing subroutines into the GASP and then 
to change the tail cone lengths by external manipulations. This 
would also allow checking the method both with and without re- 
balancing the wing. Work was being accomplished to this end 
before the writing of this report but, due to the difficulties 
encountered in trying to put the GASP into an operational status, 
could not be completed- 

The technical monitor of this project at NASA-Ames, Tom Galloway, 
performed some preliminary calculations using the GASP to test the 
method’s application. The empennage sizing method in this case was 
the V method currently used by the GASP. The results from Tom 
Galloway’s study (Reference 16) indicate that the increasing weight 
of the tail cone structure with increasing length will more than 
overcome the decreasing weight of the empennage. Drag, on the 
other hand, did decrease as expected. The intent of this procedure, 
it must be remembered, was to minimize DOC. However, DOG is much 
more dependent upon weight than drag- Consequently, Tom Galloway’s 
results indicate that a shorter tail cone is better from a DOC 
point of view. Nevertheless, we believe that by revising the GASP 
in the manner suggested above to implement the Roskam/Fillman 
method, studies in this area might prove beneficial. 

2. A Baggage Compartment Study 

It is necessary to allow enough room for passenger baggage* 

Since the shorthaul/commuter airplane is frequently used to 


2-58 


X =^+ h + w 



Figure 2.30 Baggage Dimensions 


Table 2.10 Baggage Allowance 


Large Suitcase x = 62 inches 

Small Suitcase x =55 inches 

Total Weight 70 pounds 

where x = i+ h + w 



transport passengers to the major airlines, it is desirable to 
allow for the same amount of baggage as the major airlines- specify. 

The problem of storage location was investigated in Reference 15, 
and three storage methods were studied, two of which were cabin 
oriented and the other, tail cone oriented- For the cabin-oriented 
storage methods the limitation was availability of required volume. 

For the tail-cone storage method the critical consideration was c-g- 
location- 

To establish the baggage allowances for shorthaul/commuter 
airplanes, several major carriers were questioned about their 
baggage allowances- It was discovered that most major carriers 
specify baggage allowance by the sum of the baggage dimensions 
and total weight rather than by volume. The major carriers allow 
two pieces of checked luggage plus one carry-on piece- These 
baggage allowances are summarized in Table 2-10- Figure 2-30 
defines the baggage dimensions- One major airline allows 5.0 
cu- ft- per passenger for its Trijet service. 

To attempt to determine a volume corresponding to the 
allowances of Table 2*10 a survey of baggage dimensions was 
performed using the catalog of a major retail chain. From the 
survey the dimensions listed in Table 2-11 were selected as 
representative. The original baggage compartment study documented 
in Reference 15 used the baggage volume indicated by Table 2.11. 

In that study the required baggage compartment volumes for 12, 21, 
and 30 passenger configurations were computed- It was assumed 
that no additional allowances were needed for carry-on baggage, 
as this baggage may be stored under the passengers’ seats. 

Table 2.11 Representative Baggage 

^ h w vol- 

3 

Large Suitcase 30 in 22 in 10 in 3.82 ft^ 

Small Suitcase 28 in 18 in 9 in 2.62 ft 

3 

Total Volume 6.44 ft 
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BAGGAGE = 4.566 ft 

The result, although slightly small, looks reasonable. 

This method will be used to size for the baggage compartment. 

2*5 Fuselage Shape Simulation Program, FUSE 

A computer program was written to simulate and design commuter 
aircraft fuselages. The purpose of this program was twO“fold. 
First, the program was intended to provide a method for modelling 
actual aircraft fuselages to provide node coordinates for use 
with finite element analysis procedures. This modelling method 
was to be simple to apply and versatile enough to be used with 
most commuter aircraft slopes. The second purpose of the program 
was to provide a means for applying and evaluation the design 
methods previously discussed in this chapter. In addition to 
these objectives for the program, an option was included to allow 
for interactive graphic display of the resulting design or 
simulation to ensure reasonable configurations* The final 
version of the program at the time of this report is intended 
to be used with a Tektronix 4014 or 4015 graphic display terminal 
using the PLOT 10 graphics software package (Reference 18) . 

The computer program, which will be referred to as FUSE, was 
written in FORTRAN IV for a Honeywell 66/60 timesharing system. 
Appendix C of this report contains a complete listing and user’s 
guide for the program. 

The program was originally documented in Reference 19. 

At that time the program was only used in batch operations 
for actual aircraft simulations and plotting of the simulated 
fuselages was performed on a Benson Lehner platter. Although 
the design mode and the interactive graphics have been added 
the modelling routines are essentially the same. Section 2.5*1 
will discuss the approach taken to model the fuselage. Section 
2.5*2 will provide the program description. Section 2,5.3 will 
discuss the program results. 
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The three methods of baggage storage considered in Reference 15 
are shown in Figures 2,30 through 2,33. Reference 17 provides a 
more detailed analysis of each method. The results of the analysis 
indicated that only the fuselage lobe baggage compartment could be 
considered feasible from a volume requirement standpoint. The 
carry-on baggage compartment would have required that the cabin be 
lengthened by approximately 50% and the tail cone compartment 
resulted in center of gravity difficulties. One very important 
factor was neglected in considering the fuselage lobe compartment, 
however* No allowance was made in the study for either the wing- 
carry through structure or the possibility of main landing gear 
storage. This would be especially critical for the smaller 
coimiiuter airplanes where the wing root chord is significant in 
comparison with cabin length. 

For these reasons another form of carry-on compartment 
was considered. Figure 2.31 shows the baggage compartment to 
be only on one side of the cabin. McDonnell Douglas found in 
Reference 6 that baggage compartments on either side of the aisle 
were more efficient from a design-to-cost point of view than the 
fuselage lobe method. It was decided to try and develop a method 
to define the length required for this type of baggage compartment. 
Also, at the suggestion of Tom Galloway, the project technical 
monitor at NASA-Ames, it was decided to limit the baggage to 
5 cu. ft. per passenger. 

It was assumed that the baggage compartment could be represented 
by a trapezoidal cross-section for either the circular or rounded- 
rectangular cabin sections. It was also assumed that an aisle of 
1.5 the normal aisle width would be maintained. The resulting 
effective baggage compartment cross-sections were to appear as 
shown in Figure 2.34. 

It will be assumed that the aisle-side of the compartment 
is essentially equal to the inside cabin height, The base 

of each side may be represented as in equation 2.46, if McDonnell 
Douglas control points are used: 
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Figure 2.31 Carry-on Baggage 
Compartment 



Figure 2.32 Fuselage Lobe 

Baggage Compartment 



Figure 2.33 Tail Cone Baggage 
Compartment 




a) Circular Cabin Section 



b) Rounded-Rectangular Cabin Section 


Figure 2.34 Definition of Effective Baggage Compartment 
Area 


W, 


be 


^ [' 


= 4 (N^ - 1) - .5N, w, + 12 


A A 


(2.46) 


By solving the cross-section equations for a rounded- 
rectangular section it is possible to find the wall-side compartment 
height for both circular and rounded rectangular sections as: 




= h - - 
c c 2 


1 - r. 




- w + 


+ r . w \ 2 
: tc c A 

^ / 

tc c ' 


(2.47) 


where , 

h = inside body height (in.) 

floor width (in.) = (N^ - 1) -i- 12 


Finally, the additional cabin length required to allow for 
5 cu. ft. of baggage per passenger may be expressed as: 


^^^u^ BAGGAGE 


720 PAX 


w. 


be 


2hc - 2 


V - /w--w “T r^ w \ 2~ 


where , 

PAX = number of passengers 


As a check of the method the baggage compartment length for 
a 30-passenger cabin of the type presented by McDonnell Douglas 
in Reference 6 was computed. The actual length was chosen as 
the average of the lengths of the compartments on either side 
of the aisle. In this case, that implied an actual baggage 
compartment length of 5.375 ft. 


Given : 
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2.5.1 Approach to Fuselage Shape Simulation 


For the majority of short haul/ commuter aircraft, the 
fuselage can be considered to be made up of three distinct 
sections: the nose, the cabin, and the tail cone, as shovm 

in Figure 2.35, This was the basis for the development of FUSE. 

The intent was to model each of these sections individually 
and then assemble them. The methods used to model these three 
fuselage sections will be explained in Sections 2. 4. 1.3 through 
2. 5. 1.5. As all of che cross-sections will be modelled in the 
same manner the cross-sections will be discussed in Section 2. 5. 1.2. 




Figure 2.35 Definition of Generalized Fuselage Sections 
2.5. 1.1 Coordinate Systems 

As FUSE was to be used to provide the nodes for finite- 
element analysis methods, it was necessary to define a book-keeping 
system by which the nodes might be located and used. As it was 
hoped to use FUSE in conjunction with the NCSU BODY program of 
Reference 20, the node numbering and coordinate system were chosen 
to be compatible with the BODY program. The coordinate system 
is such that x is positive aft, y is positive left, and z is 
positive up. Tbii aircraft is assumed to be symmetric about the 


XZ-plane* Therefore, nodes are only simulated for half of the 
fuselage on the positively side* 

The fuselage will be subdivided both in the lengthwise 
direction, and radially about some specified c;;ntral axis. 

The nodes will be numbered ( i,j ) in a coordinate system, 
where the value of i represents the number of the corresponding 
lengthwise subdivision and the value of j represents the number 
of the radial segment- For each node, (i,j) values for x, y, and 
z will be computed and stored* 

The NCSU BODY program is written so that the number of 
lengthwise and radial subdivisions are assigned to each fuselage 
section (nose, cabin, or tail) individually as input parameters. 

FUSE was written such that the number of lengthwise divisions 
are input for each section, but that the number of radial divisions 
is constant for all sections. This facilitates the definition 
of each panel. 

2*5*1. 2 Cross-Section Determination , CRSSEC 

The original program as documented in Reference 19 assumed 
that the airplane cross-sections could be modelled by elliptical 
or circular sections. Also the cross-sections for the nose, cabin 
and tail cone were determined independently. This sometimes 
resulted in a discontinuity at the juncture of two fuselage 
section::. To avoid this problem and to allow for round-rectangular 
cross-sections, the subroutine CRSSEC was substituted for a""' cross- 
section calculations. 

The approach taken to locating the cross-section nodes is 
similar to that used in Reference 19- The main difference is 
that the method has been generalised to account for the use of 
rounded rectangular cross-sections. 

Figure 2.36 will be used to help describe the node location 
procedure. It will be assumed that the cross-section center is 
vertically offset from the y-z origin by some value, z^. The 


nodes will be located by rotating a radial line through 180° 

(from -90° to +90° relative to the y-axis) from a center point 
at y = 0, z = z . The nodes will be located at the intersection 

CLi 

, 'll- 

of the radial and the cross-section shape every — degrees, 
where "m” is the number of radial segments desired. 



Figure 2.36 Wode Location Method 


The equation for the radial is as expressed in equation 2.49; 


z = by + z 


where, 


b = tan 


CL 

L- 1 

m 


(2.49) 


■-!] 


for j - 1, 2, 3 * , , , (ra + 1) 


m = number of radial segments 


Note from Figure 2.34 that for z-<z<z*y= y . Similarly » 

° sb st m 
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for b < 0 and y<y-,z=z -z, and for b > 0 and V < y 

z = z + z * The values for y z y , ^ and z , are easily 
o m St -^sb' ^ ^ 

determined by equations 2,50: 


^st = 

^m 

(1 - r^) 


(2.50a) 

z ^ = 

z 

(1 - c, ) 

+ z 

(2.50b) 

St 

m 

t 

o 

^sb = 

^m 

(1 - r^) 


(2.50c) 

sb 

^m 

<'b - 

+ z 

o 

(2.50d) 

where r 

^ and r, are 
t b 

the upper 

and lower round-off radii expressed 


as a eonstaflt fraction of y. and z . 

‘'m m 

This reduces the node location problem significantly. Only the 
rounded corners of the section must be dealt with. By solving the 
cross-section corner equations at either the top or bottom corner 
it can be shown that: 


y ^ 

where , 

A - 1 + 


- 3 V - 4AC 
2 A 


(2.51) 


m 


B = -2y’ + 2b ( 2 ^, - z’) 


C = 2y y’ + 
m 


and where, 


^CL” ^ 


for b > 0 
for b < 0 
for b > 0 
for b < 0 


-1 (yj' 


With a value for y at any radial position from equation 2.51, 
the value for z is easily determined from equation 2.49. 
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An added advantage o£ the subroutine CRSSEC is that the 
discontinuities that could occur at the juncture of two fuselage 
sections may be "faired-out.” This has been done at the nose- 
to-cabin juncture in the following manner. At any nose section 
radial the node coordinates of both the nose and cabin are 
determined as if the cabin shell hai been extended forward over 
the nose. The coordinates x?hich ar : the least distant from the 
center line of the fuselage are use< . This is not implemented 
for the tail ccne to allow for upswt pt tails. 

2. 5. 1.3 Nose Shape 


The aircraft nose will be modeled as the locus of two 
superimposed elliptical cones as shown in Figure 2.23. Each 
of the elliptical cones is constructed such that the top and 
side view planforras may be modeled by the generalized expression 
of equation 2.51. 

({)* •©■ • ‘ 

Figure 2.37 defines the geometric parameters used to describe 
the shape of the nose. Using the parameters shox-m in Figure 2,37 
the coordinates at any point on the nose will be determined by the 
following procedure. 

Each of the elliptical cones used for the nose will be such 
that the top and side view planforms may be described by the 
general equation 2.51. Appendix A of this report explains the 
use of elliptical cones as described in Reference 9. Briefly, 
the shape of the elliptical cone is determined as follows. 

Given a specific cone shape, a rectangle may be superimposed 
about half the planform as shown in Figure 2.38. By constructing 
the diagonals, OS and OT, it is possible to define the shape 
parameters and <j> 2 . In general, the shape parameter may 
be considered to be a taper shape parameter, while may be 




defined as the bluntness parameter for the cone. These shape 
parameters may be used to determine the exponents of equation 
2.51 as described in Appendix A. 

It will be assumed that Cone N1 may be described by its top 
and side view shape parameters and that Cone N2 may be described 
by only its side view shape parameters. The method for determining 
these shape parameters for FUSE is discussed in detail in Appendix C. 

After determining the proper shape parameters and the corre- 
sponding exponents for equation 2.51, equation 2.51 may be applied 
to determine the planform shapes for each cone. The planform 
equations for Cones N1 and N2 are expressed in equations 2.52 
through 2,55, as derived from equation 2.51. 



For Cone N1 


Side View: 


z 


ml’x 



(2.52) 
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Top View: 


ymil* - 


*N1 - == 


where: and are derived from Cone N1 side 

view shape parameters and 

and; m , and n - are derived from Cone N1 top 
yl yl 


For Cone N2 


view shape parameters, and ^^ 2 * 


Side View: 




for X > (fHX-“N2> 


for X < 


Top View: 


r( 


Si - * 


for X > (^arSa^ 


for X < Wsr^Nz) 


where; ^2 derived from Ccne N2 side 

view shape parameters. 


With the planform values, y and z at any x, for both cones, 

m 

the cross-section shapes may be superimposed according to the 
methods of Section -,5,1,2, As stated in Section 2, 5, 1,1, Che 
cross-section at each lengthwise segment will be divided into 
a specified number of radial segments. To attempt to place the 
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nodes as effectively as possible, for use in the NCSU BODY 
program, the following method will be used. 

Figure 2.39 pictures the orientation of cones N1 and N2. 
A line, CL, has been constructed from the tip of Cone N1 to 
the centroid at the base of Cone N2. 


z 



Figure 2.39 Geometric Definition of the Centerline 
for Radial Divisions 

The equation for CL may be expressed as in equation 2.56. 


'CL 


ft)- 


+ z 


NO 


(2.56) 


This line will be used as a centerline from which the radial 
divisions will be constructed. 

For X ^ nodes will only be determined for Cone Nl. 

For X > (£., 


- S. , the distance of the j'th node for each cone 

4.1 ^ IN ^ 

from CL will be used to decide which node will be retained. The 
distance will be determined as: 


- V 


y2 + (2 - 


(2.57) 
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The node resulting in the greater value o£ r will be used. 


2.5*1. 4 Cabin (Utility Section) Shape 

The cabin will be modelled by a cylinder of a constant: 
cross-section. No offsets will be used for the cabin. The 
length of the cabin will be represented as The cross- 

section at any point may be represented by the method of Section 

2 . 5 . 1 , 2 , 

2. 5. 1.5 Tail Cone Shape 


The tail cone shape will be modelled as an elliptical cone 
in much the same manner as Cone N1 of the nose. At present, 
the shape parameters used by FUSE for the tail cone are the 
same for both the top and side planform shapes. For this reason, 
an average of the top and the side view shape parameters is 
advised for input data to the program. This is discussed in 
Appendix C in more detail. 

With the shape parameters and subsequently the exponents 
for the equation, the planform shapes may be expressed as in 
equations 2.58 and 2,59. Figure 2.40 describes the geometric 
parameters used to describe the tail cone. 



(2.58) 


(2.59) 


As with Cone Nl of the nose, to allow for upswept tail cone, 
provision has been made in FUSE for a vertical offset for the tail 
cone of the form; 
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= q 


cl 


(x') + 


The cross-section will, again, be determined according to the 
method of Section 2.5*1. 2. 

2.5.2 Program Description 

Using the model developed in Section 2.51 a program was 
written in the Fortran IV computer language to compute the 
coordinates of nodes for finite-element analysis applications 
for coEmnuter aircraft fuselages. Figure 2.41 presents a simplified 
flowchart of the program. A complete listing for the program is 
presented in Appendix C of this report. 

The program calls one subroutine to aid in the computation 
of the node coordinates, CONSHP. CONSHP is a short iteration 
subroutine used to determine the planform equation exponents 
from the elliptical cone shape parameters, as stated in Section 
2.5.1. A description of the method used in CONSHP is presented 
as part of the explanation of the elliptical cone method in 
Appendix A. 

It should be noted that fuselage node coordinates are stored 
in a three-dimensional matrix, 3FUS ( I, J, K) . SFUS has been 
dimensioned as a 60 x 30 x 3 matrix. This allows a maximum of 
1,711 panels for each fuselage. Although this may appear to 
be an excessive number of panels, the dimensioning was chosen 
to allow maximum flexibility in choosing lengthwise and radial 
segment distributions for structural or aerodynamic applications. 
The I value is the number of the lengthwise segment. The J value 
is the number of the radial segment, numbered from bottom to top. 
The K value is the x, y, or z coordinate, with the K ^’alues of 
1, 2, and 3 corresponding to x, y, and z, respectively. 

Although at present, they are not output, the coordinates 
for the nose cones N1 and N2 are also stored in three-dimensional 






Figure 2.41 (continued) , 
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Figure 2.41 (continued) 


2.80 




arrays. These arrays, SN0S1(I, J, K) and SN0S2(I, J, K) , are 
set up in a similar manner to SFUS. The intent is to eventually 
evaluate the individual nose shells for crew compartment volume, 
visibility, landing gear stowage, etc. 

The program may be operated in either a design mode or a 
simulation mode. In the simulation mode, the program will read 
shape parameters for a particular aircraft configuration from 
a data file already stored on disc. The logical unit number 
for the READ statements in this case is "03." In the design 
mode the program will ask the user to input the necessary variables 
in a quest ion-and-answer process. 

One of the original intents for FUSE was to prepare data for 
use with the NCSU BODY program (or Reference 20) to enable 
parametric studies for fuselage drag optimization. For this 
reason the program will re-orient the fuselage coordinates and 
output them to a disc permfile in a format useful to the BODY 
program, if desired. The logical unit number for the WRITE 
statements in this case is "07." Several other output options 
are available as well. These are discussed in Appendix C. 

2.5.3 Program Results 

To test FUSE, several simulations of actual aircraft were made. 
Two of these simulations are presented here as Figures 2.42 and 2.43. 
Figure 2.42 represents a simulation of the Gates-Learjet Model 35/36. 
Figure 2.43 represents a simulation of the Fokker-VFIiJ F28 Mk. 4000. 

In both cases, the simulation seems to represent the fuselage 
shape with reasonable accuracy. Discrepancies, however, are apparent. 
Because of the straight-line-segment method used, for instance, the 
nose and tail become pointed. The accuracy by which the windshield/ 
body intersection is modelled depends upon the number of panels 
specified for the nose section of the fuselage. Note, ai.so, that 
the "flat" portion of the F28 tail cone has been averaged into the 
overall rounded shape. The wismatch of tail cone lengths is a 
result of a miscalculation for input data. 
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Figure 2.42 Gates Learjet 35/36 Simulation -- 580 Panels 
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It is believed that when the analysis of the fuselages 
is completed using BODY, the results will be close. Also, at 
that time, to check the sensitivity of the method slight modi- 
fications in the shape and geometry will be made to document the 
effects of simulation discrepancies on aerodynamic predictions 
of BODY. The most considerable problem to be encountered was 
the choosing of the fuselage divisions to insure accuracy 
according to the constraints of Reference 20. 

Although several conceptual designs have also been made, 
hard-copy plots were not available. Preliminary indications 
are that FUSE-generated designs are reasonable. 
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CHAPTER 3 WING CONFIGURATION STUDIES 


The project objectives as outlined in Chapter 1 called for 
performing a short study of wing si 2 ing methods and considering 
a rational approach to wing placement. Although it has not been 
possible to completely meet the objectives in this area, this 
chapter will discuss briefly one alternative method to the GASP 
method for wing sizing, and a method that was considered for wing 
placement. Section 3.1 xd.ll discuss the wing sizing method, and 
Section 3.2 will describe the xdng placement method. 

3 . 1 Wing Sizing 

It was proposed to consider the methods of Laurencf^ Loftin 
of NASA-Langley (References 3, 21, 22, and 23) as an alternative 
method to sizing the aircraft wing for optimization studies. Time 
permitted only a quick look at the methods involved, and no formal 
evaluation was possible. 

The methods described by References 3, 21, 22, and 23 may be 
applied to either jet-propelled or propeller-driven aircraft. The 
methods used are baaed on data collected from several aircraft in 
each category. 

The propeller-driven aircraft sizing methods were based on 
the characteristics of over one hundred forty different aircraft 
in a gross weight range from approximately 1000 lbs. to over 
100,000 lbs. The maximum speed range of these aircraft covered a 
range from about 100 MPH to over 500 MPH. The methods are Intended 
to size the aircraft to one or more of the following performance 
objectives (Reference 21); 

A. Airport Performance 

1. Stalling Speed 

2. Landing Field Length 


3, Take-off Field Length 

4. Climb Performance 
B. Cruise Performance 

1. Maximum or cruising speed usually at a specified I 

I 

altitude and power setting I 

1 

2. Range, again at a specified altitude and power ! 

i 

setting 1 

{ 

3. Payload j 

i 

By specifying these performance objectives, and by using the I 

aircraft data and the analysis methods proposed, it is possible j 

I 

to rapidly estimate the following aircraft characteristics; | 

Gross Weight i 

Empty Weight 
Fuel Weight 

Wing Area and Wing Loading 
Power and Power Loading 

Performance Characteristics at Values of Altitude and 
Power other than those Specified 

Figure 3.1 provides a simplified flowchart of the method as 
applied to propeller-driven aircraft. 

The performance and size estimation method for jet-propelled 
aircraft was accomplished in a similar manner. The methods were 
based on the characteristics of approximately 35 aircraft with 
gross weight ranging from about 10,500 to 800,000 lbs. For the 
jet-propelled aircraft the following performance objectives were 
considered (Reference 22) : 

A. Airport Performance 

1. FAR (Federal Air Regulations) landing field 
length including missed approach requirement 

2. FAR take-off field length including second segment 
climb requirement 
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Figure 3,1 Aircraft Sizing Flow Diagram for Propeller-Driven 
Aircraft (Reproduced from Reference 21) 






B, 


Cruise Performance 


1, Cruise speed and altitude 

2 , Range 

3 , Payload 

In this case by using the performance objectives, the aircraft 
data, and the proposed analysis methods, the following aircraft 
characteristics may be rapidly estimated: 

Gross Weight 
Fuel Weight 

Wing Area and Wing Loading 

Thrust and Thrust Loading 

Altitude for Cruise, if not specified 

Figure 3,2 provides a simplified flowchart of the methods 
as applied to jet-propelled aircraft. 

Although a formal evaluation of the methods was not made, 
a preliminary study of the method indicated that very good 
results could be expected. The intent of Loftin’ s method x^as 
to allow one person to rapidly estimate the performance and 
size of an aircraft without having to rely on expensive computer 
aids. As a result the method is very much oriented towards 
graphical approximations. The GASP, on the other hand, accomplishes 
the same basic objectives through computer applications. At 
this time, it is not felt necessary to pursue the evaluations 
of Loftin’ s approach for use in the design optimization of 
commuter airplanes* 

3.2 Wing Placement 


When the proposal for the first continuation of this 
research project was submitted to NASA-Ames in December 1976 
(Reference 2) , a misunderstanding of the wing location method 
used by the GASP inclined KU~FRL personnel to believe that an 
improved method was available* This method, provided in 
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Figure 3.2 Aircraft Sizing Flow Diagram for Jet- 
Propelled Aircraft 
(Reproduced from Reference 22) 
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Chapter 8 of Ueference 9, however, was in fact essentially 
what the GASP already used. For this reason v;ork to Incorporate 
a wing location routine based on Reference 9 was discontinued. 

One other aspect of the GASP wing location routine did raise 
some questions. In locating the wing to provide a specified 
static margin value, the GASP assumes chat the aerodynamic 
center of the wing (alone) is at 25% of the MAC. This is 
satisfactory for essentially unswept wings ^ 10 ^) at low 

Mach numbers (M < 0.3). For business jet configurations this 
is not necessarily the case* The methods of Reference 14 were 
to be used to locate the aerodynamic center, given a specific 
wing geometry at cruise Mach number, but time did not permit 
the completion of the associated routines. Also, it appeared 
that, in the GASP, it was assumed that the effect of the body 
is to shift the aerodynamic center aft. Subsequently it was 
assumed in the GASP that any static margin allowed would be 
conservative. On the other hand, the effect of the body is, 
in fact, to shift the aerodynamic center forward, tending to 
destabilize the aircraft. For this reason it was intended 
to use the Multhopp strip integration routine mentioned in 
Section 2*3 to compute the body effect and then to iterate 
to find the wing location that would result in the proper 
static margin. Figure 3.3 presents a flow diagram to illustrate 
this approach. 

As the GASP is still inoperational, these methods have not 
yet been implemented. 
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CHAPTER 4 l^TTED AREA, DRAG, AND WEIGHT STUDIES 


Tox^jards the objective of determining a minimum weight and drag 
configuration to aid in minimizing DOC, studies were made to determine 
whether a minimum drag cuaf i,;: ;::‘atioa could be defined for the fuselage. 
Weight and drag studies were both to have been conducted using the GASP, 
but due to the inoperational status of the program this has not been 
possible. It was also hoped that additional drag data could have been 
compiled to compare with the GASP results by using the NCSU BODY 
program of Reference 20, Higher priorities had to be placed on the 
GASP program, and therefore the BODY program transliteration process 
was put aside. 

Wetted area studies were conducted using FUSE for both actual 
and conceptual aircraft* As wetted area may be related to both drag 
and weight, some implication of the wetted area study may be used 
to predict what might be expected with drag and weight. The x^7etted 
area study will be discussed in Section 4.1. Section 4.2 redocuments 
the slipstream drag study presented in Reference 15- 

4.1 Wetted Area Studies 

The zero-lift drag of an airplane is proportional to its wetted 
area as shorn in equation 2.4. The weight of the fuselage structure 
may be said to be proportional to its wetted area to the power 1-2, 
as shown in equation 2.21. Therefore by studying the wetted areas 
of various airplane configurations, it should become possible to 
derive some conclusions as to the trends in weight and drag to be 
expected. 

Early in the research program, an attempt was made to correlate 
both total fuselage wetted area and fuselage section (nose, cabin, 
and tail) wetted area with their characteristic length. Data were 
compiled using wetted areas either acquired directly from the 
airframe manufacturer, or acquired by estimations from available 
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drawings. In an attempt to derive mathematical relationships, 
linear, exponential, and logarithmic least-squares curve-fits 
were applied to these data. These relationships are presented 
as Figures 4,1 through 4.12. The results of these correlations 
were not very promising. 

With the completion of FUSE (Section 2,5) another approach 
to determining wetted areas became available. Several aircraft 
were modelled using FUSE and the resultant wetted areas were 
compared with their actual wetted areas. These results are pre- 
sented as Figure 4.13. With the exception of the Fokker F28, 

FUSE seems to provide good results. 

Having determined that the wetted areas computed by FUSE 
were reasonable, a quick study was performed to determine the 
effect of cabin seating arrangement and total passenger capacity 
on fuselage wetted area. Conceptual aircraft were designed by 
FUSE to consider 10, 20, and 30 passengers in 2, 3, and 4 seats 
abreast configurations. Also these designs were made for circular 
cross-sections and for round-rectangular sections with a round-off 
radii value of 0.5. In all cases the default values for a 
piston aircraft were used for the shape parameters and ratios. 
Adequate comfort level values were assumed, and baggage compart- 
ments allowing for 5 cu. ft. of baggage per passenger were also 
included. The significant parameters for each configuration are 
presented in Table 4.1. Figure 4.14 presents a plot of fuselage 
wetted area as a function of passengers for the circular fuselage 
configurations. Similarly, Figure 4.15 presents the results for 
the rounded-rectangular fuselage configurations. 

From Figure 4.14 note that for all of the passenger capacities 
considered the number of seats abreast seems to be the controlling 
factor. For the twenty to thirty passenger range, however, the 


In figures 4.1 through 4.12, gross wetted area is the wetted 
area of the entire fuselage; net wetted area is gross wetted area 
minus the area of wing-body and tail-body intersections. 

_ 

See footnote page 
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Table 4,1 Conceptual Aircraft Design Parameters 
for the Wetted Area Study 


Note ; 1) Adequate Comfort Level is Assumed 


2) Baggage Compartment ^ 5 cu, ft, /passenger 


Conceptual 
Conf ig. 

Total 

Passengers 

Seats 

Abreast 

Round-off 

Radii 

Outside 

Width 

(ft) 

Outside 

Height 

(ft) 

Nose 

Length 

(ft) 

Cabin 

Length 

(ft) 

Tail 

Length 

(ft) 

Wetted 

Area 

(ft2) 

A 

10 

2 

1.0 

6. AO 

6.40 

9.19 

17.57 

22.38 

878.17 

B 

20 

2 

1.0 

6.40 

6.40 

9.19 

35,13 

22.38 

1237.46 

C 

30 

2 

1.0 

6.40 

6.40 

9.19 

52.70 

22.38 

1494.43 

D 

10 

3 

1.0 

7.71 

7.71 

10.87 

12.43 

27.00 

1017.94 

E 

20 

3 

1.0 

7.71 

7.71 

10.87 

22.36 

27.00 

1264.58 

F 

30 

3 

1.0 

7.71 

7.71 

10.87 

32.29 

27.00 

1526.47 

^ G 

10 

A 

1.0 

9.20 

9.20 

12.09 

8.98 

34.11 

1331.11 

H 

20 

A 

1.0 

9.20 

9.20 

12,09 

15.47 

32.19 

1403.99 

I 

30 

A 

1.0 

9.20 

9.20 

12.09 

24.45 

32.19 

1683.11 

J 

10 

2 

0.5 

5.25 

6. 74 

10.55 

17.68 

20.98 

929.07 

. K 

20 

2 

0.5 

5.25 

6.74 

10.55 

35.35 

20.98 

1328.55 

L 

30 

2 

0.5 

5.25 

6.74 

10.55 

53.03 

20.98 

1686.99 

M 

10 

3 

0.5 

6.92 

6.88 

11.54 

12.78 

24.15 

955.47 

N 

20 

3 

0.5 

6.92 

6.88 

11.54 

23.06 

24.15 

1201.46 

0 

30 

3 

0.5 

6.92 

6.88 

11.54 

33.33 

24.15 

1470.86 

P 

10 

A 

0.5 

8.58 

6.99 

11.78 

9.40 

27.26 

1017.50 

Q 

20 

4 

0.5 

8.58 

6.99 

11.78 

16.30 

27.26 

1198.44 

R 

30 

A 

0.5 

8.58 

6.99 

11.78 

25.70 

27.26 

1460.11 


Net Wetted Fuselage Area (ft. 






Figure 4.2 Exponential Correlation Between Fuselage 
Length and Fuselage Wetted Area 


Net Wetted Fuselage Area (ft. 


lO 


1) Data from Ref. 17 

2) Logarithmic Curve Fit; 
Swof = “1894 + 665. 9/n(£D) 


0 40 

Fuselage Length, (ft.) 


Figure 4*3 Logarithmic Correlation between Fuselage Length 
and Fuselage Net Wetted Area 






^Figure 4*4 Linear Correlation Between Nose Length 
and Nose Wetted Area 
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Figure 4*5 Exponential Correlation Between Nose 
Length and Nose Wetted Area 


Net Wetted Nose Area (ft. ) 





'Figure 4.7 


Linear Correlation Between Cabin Length and 
Cabin Wetted Area 


Net Wetted Cabin Area (ft. ) 






Net Wetted Cabin Area (ft. ) 



Length of Cabin Section, (ft) 


Figure 4,9 Logarithmic Correlation Bet\i?een Cabin 
Length and Cabin Wetted Area 




Net Wetted Tail Cone Area ~ (ft?) 


NOTE: 

1) Data from Ret. 17 

2) Linear Curve Fit; 
Swet - 9. 33 Ig- 3. 92 

r^-0.90 


Wose Cabin 





Length of Tail Cone, ^ (ft) 


Figure 4,10 Linear Correlation Between Tail Cone 
Length and Tail Cone Wetted Area 




Figure 4*11 Exponential Correlation Between Tail 
Cone Length and Tail Cone Wetted Area 



Figure A. 12 Logarithmic Correlation Between Tail Cone 
Length and Tail Cone Wetted Area 


Predicted Wetted Area, Sv^,et ^Fuse), ft 


Fokker F-28 


Fokker F-27 


1600 


Lockheed Jetstar 

Ae;?o Commander 690A 
J3€tes-Lear]et 35/36 I _ 


Beech Baron 58 


0 800 1200 1600 200 
Actual Fuselage Wetted Area, S^et (Act), ft.^ 


figure 4.13 Comparison of Actual Fuselage Wetted Areas 
with Wetted Areas Predicted by FdSE 







Fuselage Wetted Area, 


1200 


1) Round-off Radii =1.0 

2) Wetted Area Predicted by Fuse 

3) Baggage - 5 Ft^/Pass. 

4) Adequate Comfort Level 

5) Seats Abreast; I I 


Number of Passengers PAX 


Figure 4.14 Gross Wetted Area as a Function of Passengers 
for Circular Fuselages 



Fuselage Wetted Area, 


2000 


1800 


1600 


J 1400 


1) Round-off Radii = 1.0' 

2) Wetted Area Predicted by Fuse 

3) Baggage 5 ft. /Pass. 

4) Adequate Comfort Level 

5) Seats Abreast: I i 


Number of Passengers -- PAX 


Figure 4.15 Gross Wetted Area as a Function of Passengers 
for Rounded-Rectangular Fuselages 


wetted areas for the 2 and 3 seats abreast configurations are very 
close. Considering drag, in this case, the 2 seats abreast config- 
uration is still b^st, A longer fuselage length results in higher 
Reynolds numbers and therefore lower skin friction coefficients. 

From a weight standpoint, as weight may be considered as propor- 
tional to wetted area to the power of 1,2, the 2 seats abreast 
configuration also prevails. 

Figure 4.15 produced some rather unexpected results. Note 
here, that the wetted areas for the 2 seats abreast configurations 
rapidly become very much greater than the other seating arrangements. 
Checking Table 4.1, it is apparent that the nose cone lengths for 
this configuration appear overly large. This might indicate that 
the standard crewbox is too large for this configuration. Looking 
at the 3 and 4 seats abreast configurations, note that 4 seats 
abreast become more efficient, from a wetted area standpoint above 
16 passengers. This could mean a drag trade-off as the 3 seats 
abreast configurations will have longer fuselages and therefore 
lower skin friction coefficients. From a weight standpoint the 
4 seats abreast configuration is still preferable. 

These assumptions and conclusions are of course based solely 
on the wetted areas. Such factors as comfort, pressurization, 
ease of construction, etc., should also be considered before 
selecting a configuration. 

4.2 Approach to the Prediction of Zero-Lift Drag in a Propeller 

Slipstream 

When a rotating propeller is working in the presence of a body 
(wing, nacelle, fuselage), the flow through the propeller disc and 
about the body will change. The flow through the propeller disc 
will change because it will work in a perturbed flow field. These 
perturbations alter the local angle of attack and effective airspeed 
at each blade section, and thus the overall propeller characteristics. 
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The flow over the body will change due to the existence 
of a slipstream. The body can be either immersed in the slip- 
stream or be close to it; in either case the airflow about the 
body will change. As a result, lift and drag of the aircraft 
parts in the slipstream and propeller thrust and power-required 
will change. 


L+AL 



Figure 4.16 Change in Aerodynamic Forces 
due to the Slipstream 


The lift dependent drag is affected by propeller power in the 
following ways; 

1. The components of propeller thrust and normal force that 
are parallel to and have the same direction as the wing 
lift reduce the wing lift required, thereby reducing the 
wing drag due to lift. A method to compute the magnitude 
of this effect can be found in Reference 24. 
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2 . 


The propeller slipstream modifies the downwash and dynamic 
pressure over portions of the wing thus changing the wing 
drag due to lift. Two different methods are available to 
compute the magnitude of this change In drag: a) Reference 

25, Section 4.6.4 (essentially empirical) and b) theoretical 
methods, for example: Reference 26* 

Sophisticated methods exist tha,t predict lift dependent drag 
changes. Therefore, the topic of this investigation will be: changes 

in zero-lift drag and propeller performance. 


4.2.1 Zero-Lift Drag 

Methods to predict the increase in zero lift drag of bodies 
immersed in a propeller slipstream are presented in several publi- 
cations. According to all these publications, the drag increase 
is proportional to the average increase in dynamic pressure in the 
slipstream and an appropriate drag area. 

Reference 24: 

AC_ = ^ Aq„ dS (4.1) 

wher«d: Aq^ = T/irR^ 

Sj = area, Immersed in slipstream 


Reference 28: = ^/(Vq + cosa^)^ + sin^a^*(4.2) 

where: = resultant velocity about body 

''i - »/''" + f - 

= average velocity increase in 
slipstream far behind propel- 
ler according to momentiom 
theory. 


Reference 27: 


’^eff 

where; 


^ ^is 

EC S 

F = 1 - 1.558 


(4.3) 
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Except for an additional correction factor, this method is the same 
as the methods of References 24 and 28 because these can also be written 
as follows: 

with; Aq = T/irR^ „ „ 

F , , 4 (4.4) 

TT 

AD = AqCjjS^ 

4 

Instead of — a factor of 1.558 is used in the Reference 27 method 

IT 

(22% higher) . In all these methods the slipstream is represented 
as shown in Figure 4.17. 



Figure 4.17 Representation of Slipstream 

4. 2. 1.1 Theoretical Considerations 

Theoretically it can be proven that the spanwise axial velocity 
distribution behind a propeller is as shown in Figure 4.18. These 
results have been verified experimentally. 


tip boundary 



Figure 4.18 Spanwise Axial Velocity Distri- 
bution Behind a Propeller 
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Figure 4.19 Curves of H/q against x for NACA cowling 


Figure 4.19 (from Reference 29) shows distributions of total 
pressure increase (H) divided by q (ambient* dynamic pressure) . 

Considering these curves it seems unlikely that both slender 
and wide bodies are equally affected (as far as drag increases are 
concerned) by a slipstream. This suspicion is reinforced by experi-- • 
mental results presented in Figure 4.20 from Reference 30. 

Apparently, the drag increase of bodies with a small propeller 
(high d/D) is higher than the increase in zero lift drag of bodies 
with a larger propeller at the same conditions (i.e. same Aqg and 
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(a) 48- in. propeller; 12-in. 

(b) 48-in. propeller; 12-in. 

(c) 48- in. propeller; 16-in. 

(d) 36-in. propeller; 12-in. 

(e) 36-in. propeller; 16-in. 


nacelle (with spinner); d/D 
nacelle; d/D = .25 
nacelle; d/D = .33 
nacelle; d/D = .33 
nacelle; d/D = .44 


0.25 


Figure 4.20 The variation of slipstream-drag coefficients with 
apparent propeller thrust-loading coefficients. 
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AD. /pV^2d 2 

lus 0 p 


T = T/pV„2d 2 
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An obvious explanation for this phenomenon has not been found. 

It appears necessary to know the development of the slip- 
stream in streamwise direction accurately, in order to provide 
better predictions of zero-lift drag increases than those presented 
in Reference 24, 

Reference 29 presents dynamic pressure contours in the plane 
of the elevator hingeline of some propeller driven airplanes. Some 
of the results are shown in Figure 4.21. 



Figure 4.21 Dynamic pressure (q/q«) contours and inclination 
of the air stream in &e plane of the elevator 
hinge line. . Vectors show deviation of air flow 
from the free-stream direction. View looking 
forward. Circle shows projection of propeller 
center . 
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These figures show the velocity distribution in a cross-section 
of the slipstream some distance behind the propeller. 

It can be seen that in this case (where the slipstream is not 
influenced by a body) the slipstream diameter some distance behind 
the propeller is still the same as the propeller diameter. The dy- 
namic pressure distribution however is quite different from the one 
right behind the propeller, suggesting that mixing processes can be 
quite important. The results shown in Figure 4.22 are similar. 

Such a slipstream development rules out quick, but accurate, 
predictions of its influence on drag. 

Detailed theoretical or experimental investigations appear neces- 
sary to yield an accurate prediction method. A simplified approach is 
given in Section 4.2. 1.2. At the moment, however, it seems best to use 
the average Aq (cross sectional) in AC calculations. . 

The use of this ”easy'^ method seems even more attractive, after 
looking at Figures 4.23a and b, showing dynamic pressure contours of 
slipstreams that have been influenced by bodies. 
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An attempt has been made to show that dynamic pressure changes 
in streamwise direction due to mixing, etc. do exist and make accurate 
drag predictions virtually impossible (i.e. in a fast and simple way). 
But, even when a very simple theory is used (Momentum Theory), neglect- 
ing all effects like mixing and viscosity, it can be shown that the 
dynamic pressure (and static pressure) change in streamwise direction 
(Figure 4.24). 



Figure 4.24 Change in Pressure Behind a Propeller 

In the region of the fuselage nose, 4qg is apparently lower than 

far behind the propeller, while Ap is higher. The influence of both 

effects on fuselage drag however can be approximated by assuming that 

Ap = 0 and Aq = constant = A far behind propeller. This assumption 
a q 

might lead to erroneous results when short fuselages with a blunt nose 
are considered. 

4.2.1.2 a Simplified Approach to Zero-Lift Drag Prediction 

An approach to the zero-lift drag prediction problem might be: 

1, Calculate velocity distribution in slipstream (using continuity 
expression) . 

2. Determine dV/ds at nacelle boundary (s = distance to nacelle 
boundary) (Inviscid) . 

Determine influence of dV/ds on boundary layer development and 
wall skin friction. 


3 . 


This approach has not been developed completely. The following expla- 
nations might help in such a development. 


If local speed at a particular blade section is V-+V , then there 

U cL 4 
X 

must be a speed V-+2V ; in the fully developed slipstream. These speeds 

U 3 • 

are along one slipstream line. 
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Figure 4.25 Velocity Profile in a Slipstream 


With the continuity expression it follows: 

pir(r^. - r^. ,) x V_^ = pir(r^ -r^ ) x V ^ 

^ pi pi“l Pav. s^ ®i-l 

^pi pi-1 


(4.6) 


r r 
s s . , 
1 1-1 


where: r , , and r . are ends (inboard and outboard) of a blade section. 

pi-1 pi 


r .-r . , 

pi pi-1 


average speed at this blade section = 


r and r are boundaries of streamline tube (boundaries of 
®i-l ®i 

which at propeller disc are: ^pi^ slipstream. 

V = average speed through this streamline tube = 
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Figure 4.26 Velocity Profile in a Streamline Tube 


Therefore ; 


(r2.-r .2 ) X (2V.+V +V )2= 
px px-l" 0 a^_j^ 



r 

5 


2 ) X (2V_+2V^ 
i-1 ° ^ 


+2V )2 

i ^i-1 


(^. 7 ) 


With the kno\<rn axial velocity distribution at the propeller disc 
(theory or experiments), the velocity distribution in the slipstream 
about a nacelle can be obtained by numerical integration. 



Figure 4,27 Velocity Distribution in the Slipstream Around 
a Nacelle. 
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At nacelle boundary: ^ ~ ^0 • 

According to boundary layer theory: 


,6,3u 3v 
f (-y- + TT-) dy = 0 
0 ^3x 3y' 


^ V{6) - V(0) = - 1^ dy 

0 9x 


N»B- 6 is distance to wall, where viscous speed u is 

equal to inviscid local speed U = speed according 
to 1 of Section 4. 2*1. 2. 


(4.8) 


S f 3 ^ / N ^ -* 




3y= 


} dy = 0 


, .5 3^u j f/3u\ 1 *^0 /3ds 

where: - v^/ — dy - - — - V (^)j 


Rewritten: 


9y^ 


- { /^ u(U-u)dy} - /^(U-u)dy - v (-|r)f 

dx 0 9y 0 oy ° 


(4.9) 


(4.10) 


Assuming that U hardly changes over the boundary layer, and that 
r 0 and can be determined by the method explained in Section 
3. 1.2.1, the increment in skin friction draf^ can be determined 
using boundary layer theory. 

A first order approximation (assuming velocity profile in boun- 
dary layer doesn't change) is: 


A /Sn,, _ /3u^ 

= I"’ '^>6 ■ »Aii 

according 
to 1 


(4.11) 


Reference 30 presents a method to compute changes in lift due to a 
propeller slipstream. In the first part of the method velocities at 
the propeller disc and in the slipstream (including swirl) are calcu- 
lated. This method is more complicated than the one explained above, 
but could also be used Instead of the method in Section 3. 1.2.1 to 
yield similar results. 

4.2.2 Propeller Blockage. 

Only Reference 27 gives a method to account for propeller block- 
age effects: 
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(4.12) 


J 


ef f 


where 


(1-h) 


J. 

xs 


J. 

IS nD 


1-h = l-.32g 


By using •'eff instead of in propeller performance charts (foi^ 

isolated propellers), thrust and power required of a propeller 
(operating in a flowfield that is perturbed by a body) can be found. 


4.2. 2,1 Theoretical Considerations. 

In the presence of a body, the axial velocity relative to a blade 
section is not the velocity of advance Vq but a velocity u caused by 
the blocking effect of the body. As a result, the geometric pitch 
will be modified: 

instead of: P/D - V/nD (4-13) 


P/D V/u = V/nD 

It means that the blade sections will work at a different angle of 
attack. 

According to Reference 29, the true propeller efficiency may be 
computed according to the following relation: 


n 


= n 

is 




V dx 


dx 


1 

/ 

0 



(Approximate;) 


(4.14) 


In order to calculate n, both the velocity distribution through the 
propeller disc and the spanwise propeller load distribution have to 
be kno^>m. Both can be found either experimentally or theoretically . 

A Blade-Element theory could be used to predict theoretical load 
distribution. Source-sink distributions can be convenient when cal- 
culating velocity distributions in the propeller disc. Figure 4.28 
shows the body shapes that can be simulated by a combination of 1 sink 
and 1 source. 

Figure 4.29 shows some u/v distributions in the propeller disc 
according to this potential flow theory. 
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4,2,2. 2 Experimental Results . 

Some of the results presented in Reference 29 are shown in Figure 4,30, 
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Figure 4,30 Comparison of apparent propeller efficiency 
envelopes for four body shapes. 



Figure 4,31 Velocity distribution of propeller (pro- 
peller removed), 93 miles per hour, 

Depending on the accuracy required, either the method of Reference 26 
can be used, or a blade element theory combined with a method that predicts 
velocity distributions in the propeller disc- 

Further studies should be conducted to develop this approa^^h to the 
estimation of slipstream drag (along with the approach of Reference 26) 
into a programmable method. Eventually the intent is to include this 
method into the drag estimation portion of the GASP. 
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CHAPTER 5 THE GENERAL AVIATION SYNTHESIS PROGRAM, GASP 


As has been stated at several points throughout this report, 
the means by which all of the design optimization methods were 
to be implemented was the General Aviation Synthesis Program (GASP) 
developed for the NASA-Ames research center by the University of 
Minnesota, The KU-FIij research staff has been working continuously 
now for approximately one year to put an up-to-date version 
of the GASP into ah operational status on the University of 
Kansas Honeywell 66/60 computer. At this time the program is 
still not implemented on the University of Kansas facilities, 
although it is believed that for reasons that will be explained 
in the next section, this will be accomplished in the near future. 

The difficulties encountered in putting the GASP into 
operational status have also prevented the KU-FRL research 
staff from making revisions to the GASP to reflect the research 
that has been completed. This chapter will therefore discuss, 
first, the steps taken in attempting to transliterate the GASP, 
in Section 5.1, and second, how the design constraints were 
to be implemented in the GASP, in Section 5.2. 

5.1 Transliteration Process 

A copy of the GASP was first received at the KU-FRL in 
June 1976. The program was forwarded to the KU-FRL by Tom Galloway 
at NASA-Ames on magnetic tape. Also provided with the magnetic tape 
were a complete listing, a copy of Reference 31 describing the 
general flow of the program, a description of input parameters, and 
a sample output. Unfortunately, the magnetic tape parity was 
incompatible with the University of Kansas computer system. The 
Tape was returned to NASA-Ames for a copy with the proper parity. 

The new magnetic tape was in the following format: 

IBM compatible 
9-track 
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• 1600 BPI 

• EBCDIC -80 BYTE records 

' No Label 

Even Parity 

Transliteration procedures were started for this tape; and 
although all compilation errors were finally eliminated, the program 
could not be coerced to execute. At that time no explanation was 
available. 

A card deck version of the GASP was acquired from the Cessna 
Aircraft Company with the assistance of Don Halverstadt. This 
version was designed to be run on an IBM 360 computer system. 

After completing preliminary transliterations it was determined 
that the free-field input routines of this version and the previous 
NASA version were incompatible with the Honeywell 66/60 system. 

The Cessna input routines proved easier to revise. In early 
March two runs were made of the Cessna version with only minor 
difficulties. 

The Cessna version of the GASP was an older version and 
consequently it became necessary to update the program to reflect 
NASA revisions. Revisions were made to the program using the 
Honeywell 66/60 timesharing system editor and by storing the 
GASP on disc subroutine by subroutine. This process was proving 
to be extremely time consuming with no guarantee of success. 

It was decided to acquire a new magnetic tape from NASA-Ames 
with the intent of splicing the revised Cessna input routines 
onto a transliterated NASA version. This was the status at 
the end of June 1977. Fiscal year end administration difficulties 
prevented the completion of this process. 

Negotiations have been carried out with the University of 
Kansas Computation Center to gain what is referred to as 
"internal project status" for the GASP. This status has now 
been assigned by the computation center, and they will implement 
the program using its more experienced staff a^ charge . 
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Previously it was understood that such status would not be possible 
for the GASP, To have the GASP implemented on an "external project 
status’* would have resulted in prohibitive costs. It is expected 
that the GASP will be finally operational by the end of September. 

5.2 Implementation of Design Constraints 

In Chapter 1 it was assumed that the aircraft configuration 
could be designed to meet four specific constraints: 

a) Utility Constraints 

b) Stability and Control Constraints 

c) Mission Constraints 

d) Performance Constraints 

The GASP already will account for the last two constraints 
adequately and therefore those constraints will not be discussed 
here. The utility and stability and control constraints were 
in need of some revision, however. The methods that were to be 
used to revise the GASP accordingly will be discussed in the 
following paragraphs. 

5.2.1 Cabin-Utility Constraints 

In Chapter 2 a number of methods for sizing for the utility 
constraints were discussed. These methods were brought together 
in the form of the design mode of the program FUSE- It was 
intended to replace the fuselage sizing portion of the subroutine 
SIZE in the GASP with a call to a subroutine FUSIZE. FUSIZE 
would have consisted of the design mode and wetted area calcu- 
lation portion of FUSE. This could have been readily implemented 
with one minor exception. The Honeywell 66/60 allows only 59 
arguments for SIZE. This would mean that to add the needed data 
for FUSIZE, a namelist input method would be necessary* This, 
however, should present no great difficulties. 


5.3 


5-2,2 Stability and Control Constraints 


I 














Two forms of stability and control constraints were considered 
for implementation in the GASP — static constraints and dynamic 
constraints. The methods that were considered for each will be 
discussed in turn. 

The static stability and control constraints that were 

considered were: longitudinal static stability in the stick- 

fixed case, implemented through C ; and directional static 

a 

stability, implemented through . Both constraints were to 

have been accomplished in the empennage sizing process in the 
following manner. The first time the SIZE is called, the weight 
and c.g- data are usually lacking. As the c,g- location is 
required for both constraints, it was intended to use V methods 
to size the empennage. The next time SIZE is called, when a 
c,g, location is known, the subroutines STABAREA and VERTAREA 
(Section 2.3) would be used to size the horizontal and vertical 
tails for the stability constraints. In this manner the prelim- 
inary use of a V method would provide the required "seeds*^ for 
the iterative processes of STABAREA and VERTAREA. 

The dynamic stability and control characteristics presented 
more of a problem. To determine the dynamic characteristics of 
the aircraft all of the non-dimensional stability and control 
derivatives, and the airplane inertias must be known. It was 
intended that rather than to attempt to constrain the design 
for dynamic characteristics, the configuration would simply 
be analyzed for dynamic stability and control characteristics. 

To accomplish this the non-dimensional derivatives would be 
determined using the methods of Reference 14, and the analysis 
would be implemented by appending the appropriate portion of 
the program described in Reference 32 to the GASP, An inertia 
determining routine would also be required, but this was not 
felt to present any great difficulty. 

These characteristics, plus a number of others, will be the 
subject of the next chapter. 
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CHAPTER 6 PROPOSED FUTURE RESEARCH 


With the exception of integrating the design methods into 
the GASPj most of the design constraints mentioned in Chapter 1 
were accounted for* Only the stability and control constraints 
have not been adequately enforced. Considering this, a proposal 
for the continuation of the project was submitted to NASA-Ames 
(Reference 33) , placing sole emphasis on stability and control 
methods* This chapter will present what is essentially a reprint 
of the stater \t of work of that proposal. 

Care will be exercised not to duplicate any work already 
done in this area by other GASP investigators. 

6.1 Objectives 


It is proposed that work on the design optimization of short 
haul and commuter airplanes be continued with the following objectives 

a) To determine those stability and control characteristics 
which are critical to the preliminary design process. 

b) To evaluate stability and control analysis methods 
currently available to determine those methods most 
appropriate for the preliminary design function which 
GASP performs. 

c) To determine how the methods of b) may be used to 
provide the proper constraints and/or analysis functions 
for GASP. 

d) To develop the appropriate subroutines for the methods 
of c) and how they may be appended into GASP. 

In line with these objectives, emphasis will be placed on 
the stability and control characteristics of both jet and propeller 
driven airplanes. Also specific attention will be given to the 
determination of the following stability and control characteristics 
in the preliminary design process: 
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a) Static longitudinal stability; 

b) Static directional stability; 

c) Engine-out control; 

d) Calculation of rotation velocity, V_; 

e) Longitudinal dynamic stability; 

f) Lateral-directional dynamic stability; 

g) Trim at low speed and forward C.G. 

The following section will outline briefly the approach that 
will be taken to meet the objectives with respect to the above 
stability and control characteristics, 

6.2 Stability and Control Analysis Methods 


The primary references that will be used to determine and 
evaluate the stability and control analysis methods will be 
References 9, 14, 24, and 34, In the following paragraphs each 
of the stability and control characteristics listed in the previous 
section and its proposed analysis method will be discussed briefly. 

6.2.1 Static Lonaitudinal Stabilit^^ 


Using the methods of References 14 and 34, the configuration 

will be analyzed for static margin (dC /dC ) , static longitudinal 

m L 

stability (C ) , and the neutral point for both and stick-fixed 
m 
a 

and stick-free cases. To facilitate the calculation of these 
values a simple Multhopp integration procedure for the body (and 
engine nacelles in the case of wing^mounted engines; see Section 
2.3) will be used. Correlations with on-hand tunnel data actual 
configurations will be made. 

Using Equation 4.35 on page 4.23 of Reference 34, the stick- 
fixed neutral point will be defined as in Equation 6.1. The stick- 
free neutral point will be defined by Equation 5-99 on page 5.47 
of Reference 34 as shown in Equation 6.2. In both cases the variables 
are as defined in Reference 34. 



Stick-Fixed Case: 
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Stick-Free Case: 
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Propeller effects will be accounted for using the method of 
Reference 34, Chapter 4. 


6.2.2 Static Directional Stability 

Values for static directional stability, C , will be computed 

using the methods of References 9 and lA. Also propell' r effects 
on directional stability will be considered using the methods of 
References 9, 24, 34, and 35, 


6.2.3 Engine-Out Control 

The methods of References 9 and 34 will be used to analyze the 
configuration for the minimum engine-out control speed, Engine- 

out control will be considered from both the single-degree-of-freedom 
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and three-degree-of-freedom points of view. Drag due to stopped 
engines and/or propeller will be accounted for. 

6.2.4 Calculation of Rotation Velocity, V_ 


The speed required to rotate on take-off, V , will be calculated 

K 

using the method of Reference 9, 

6*2.5 Dynamic Longitudinal Stability 

Dynamic longitudinal stability characteristics will be considered 
using the methods of both References 9 and 34 for the stick-fixed case. 
The dynamic longitudinal stability characteristics upon which primary 
emphasis will be placed will be the short period damping and undamped 
natural frequency. Phugoid damping and frequency will also be 
considered* The methods considered will include both the two- 
degree-of-freedom short period and phugoid approximations and 
the complete three-degree-of-freedom solutions of Reference 34. 

Also the relatively simple dynamic stability relationships of 
Torenbeek in Chapter 9 of Reference 9 will be investigated as 
to their validity. 

One possible method for analyzing a configuration for dynamic 
stability in the GASP that will be considered will be a revised 
version of the dynamic stability and control analysis program 
document in Reference 32. 

6.2.6 Lateral-Directional Dynamic Stability 

In analyzing the dynamic lateral-directional stability the 
characteristics of the spiral, roll, and dutch roll modes will 
be considered. Using the methods of Reference 34, these methods 
will be analyzed from the approximate and complete three-degree- 
of-freedom points of view. 
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APPENDIX A USE OF POLYNOMIALS WITH 


T 

^ FRACTIONAL EXPONENTS TO APPROXIMATE 

EXTERNAL AIRCRAFT LINES 

Appendix B of Reference 9 describes a method for 
approximating the external lines of aircraft for the 
estimation of wetted areas, volumes, cross-sections, 
etc. This appendix will discuss that method and how 
it has been applied to fuselage shape simulation. . ■ 

A1 General Background 

In Reference 9, Torenbeek uses polynomials with 
fractional exponents to approximate the e.xternal lines 
of aircraft fuselages. These polynomials are of the 
form expressed in equation A-1. 

(if (ef = i; n > 1 (A-i) 

Figure A-1 presents a generalised curve of this 
type. In Reference 9, Tc~enbeek assumes that the 
plan view curves of aircraft may be generalized by 
equation A-1. As shown in Figure A-1, a shape para- 
n meter (|) may be determined from the plan view of an 

aircraft's external lines. This shape parameter may 
then be used to determine certain factors, k^ , k^, k^ 
and k^, according to Figure A-2, These k-f actors may 
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m 



Figure A-1 Definition of Shape Parameter, for the 
Generalized Curve 



_ Area of Section OBSA 
^a" a b 

. Circumference BSA 
"c a + b 


k\/= Volume of Bod y of Revolution 

Tia b2 


If = Wetted Area of Body of Revolution 
^ 277ab 


0 


OP 

OA 


0 


0 

Figure A-2 Nondiinensional Geometric Constants Defined 
by <J) 


^fmax 



Figure A-3 Fuselage Geometry Used to Define Characteristic 
Areas and Volume 


be used to calculate the characteristic areas and 


volume. Using the geometry of Figure A-3: 
Cross-sectional area*: 

Ap = k. bo ho (A-2) 

■^max max 

Circumferential length of the cross-section: 

Co = 2kp (bo + ho ) (A-3) 

max max 

Fuselage volume; 

Vj “ ^n ‘‘‘ (A-4) 

n t 

Fuselage wetted area: 


wet 


- 


■‘w In 
n 


+ 




(A-5) 


It may be noted that Torenbeek's method precludes 
the necessity of knowing the values of m and n in 
equation A-i. To use these curves to simulate fuselage 
shapes, however the values for m and n become a pre- 
requisite . 


*At the fuselage station where the width and height 
are maximum 
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As the intent was to develop a program to aid in 
the preliminary design of fuselage shapes, it was 
desired to find a method to specify cone shapes for 
the nose and tail that was relatively easy to use. 

The use of m and n in this regard proved to be rather 
cumbersome, as it was difficult to develop a 'feel* 
for their values. As an alternative, the shape para- 
meter, (J), was chosen to specify the shapes. In addition 
a second shape parameter was defined to enable the 
determination of the values of m and n within the pro- 
gram. Figure A-4 defines the geometry of the two 
shape parameters, <j>^ and 



Figure A-4 Definition of Elliptical Cone 
Shape Parameters 
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From Figure A-4 it may be seen that may be con- 
sidered, to indicate the taper characteristics of the 
cone, while (ji^ indicates the bluntness characteristics 
of the cone. 

These shape parameters may be used to determine 
the values of m and n in the following manner. From 
Figure A-4 it may be seen that the line ^ may be 
represented by the equation A-6. 

y = |x (A-6) 

Similarly the line OT may be represented by 
equation (A-7) . 



Substituting for y in equation A-1 gives the 
following: 


(If * (if - 1 

(-a)“ - (i)- ^ 

At the points S and T, x equals <^^a and <^ 2 ^ 
respectively. Substituting these values for x into 
equations A-S a and b respectively renders: 


(A-Sa) 

(A-Sb) 
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1 


(A-9a) 


(<f>2) + (t) " ^ (A-9b) 

By solving these two equations for m: 

mil - 

m " 

In 

Inll - 


Equations A-10 may be used, to solve for m and n 
by iteration. Such a method was used in the subroutine 
CONSHP, for which the listing is presented in Section 
A3. It should be noted here that the exponents m 
and n are purely a function of the shape parameters 
and a completely independent of the cone fineness 
ratio, a/2b. 

A2 Example Cone Shapes 

To provide a 'feel' for the cone shapes resulting 
from various shape parameters a program was written 
for a Hewlett-Packard 9100B desktop calculator that 
would compute the exponents for and plot various cone 
shapes. Figures A-5 and A-6 provide cone shapes for 


(A-lOa) 


(A-lOb) 
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cones of fineness ratios (a/2b) of 1.5 and 2.0 
respectively where is held constant and (j >2 is 
varied. Table A-1 provides the values of m and n for 
these figures, Figures A-7 and A-8 provide cone shapes 
for cone fineness of 1.5 and 2, where (jig has been 
held constant and (})^ is varied. Table A-2 presents 
the values of m and n for these figures. 

It is interesting to note that for each value of 
<j)^, there is a practical range for cpg* If <p 2 is too . 
small, the value of m drops below 1.0 and an inflection 
occurs near the nose of the cone. If 4>2 is too large, 
the value of n drops below 1.0 and an inflection occurs 
near the base of the cone. It has been found that the 
windshield shape can sometimes best be simulated by 
the former of these two cases. The inflection that 
occurs in cone N2 will normally be buried in cone N1 . 
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Table A-1 Values of m and n for 


Various Shape Parameters 


^1 

'^2 

m 

n 

.55 

.65 

.2350 

3.3990 


.70 

.7523 

1.6987 


.725 

1.1320 

1.1874 


.75 

1.5982 

.8118 

.60 

.70 

.2375 

4.2468 


.75 

.7808 

2.1769 


.80 

1.6865 

1.0749 


tn 

CO 

3.0995 

.4498 

.70 

o 

00 

.2970 

6.4413 


.85 

1.0216 

3.3258 


.90 . 

2.3345 

1.6001 


.95 

4.7857 

.5613 

.75 

.85 

.3810 

7.8740 


.875 

.7802 

5.5762 


.90 

1.3559 

3.9285 


.95 

3. 3668 

1.6596 



Figure A-7 


Various Cone Shapes 
for a/2b = 1.5 


ijid 
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a 

d) ((>2 = 0.85 


Figure A-8 Various Cone Shapes 
for a/2b - 2,0 
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Table A-2 Values 

of m and n for 



Various Shape 

Parameters 


^^2 


m 

n 

.70 

.55 

.7523 

1.6987 


.575 

.4482 

2.7038 


.60 

.2375 

4.2468 


.625 

.0902 

6 . 7696 

.75 

.55 

1.5988 

.8118 


.575 

1.1494 

1.3620 


.60 

.7808 

2.1769 


.625 

.4841 

3.3876 

.80 

.55 

2.9056 

.3239 


.60 

1.6885 

1.0734 


.65 

.8597 

2 . 7225 


.70 

.2970 

6.4413 

.85 

.60 

3.0982 

.4502 


.65 

1.8974 

1.3526 


.70 

1.0216 

3.3258 


.75 

.3834 

7.8530 

.90 

. o6 

3.5879 

.5566 


.70 

2.3345 

1.6001 


.75 

1.3559 

3.9285 


.80 

.5768 

9.4732 

.95 

.70 

4.7857 

.5613 


.75 

3.3668 

1.6596 


.80 

2.1867 

4.2647 


.85 

1.1343 

10.9631 
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A3 Listing for the Subroutine C0N5HP 


This section contains a listing of the iterative 
subroutine used to determine the values of m and n. 



1FCJ.EQ.100) GO TO 80 
*' DM1 a H1-M2 “ 

IFCDMI.LE.O.) GO TO 30 
ifc'omi.gY.dhS) go TO so 
0 N = N+C0.05/I) 
bM2 =''0Ht 
GO TO 10 

Q iFCAasrbHlTTLrrreOGO fo~30 


Ml = ALOGM.Q -CPHI1 **N)) /ALOGCPHII ) 

K'2 a'AVoG'cr.Q - fpiii2**N)rrAi.0G(PHi2/i7cn 
J a J+1 

r F“( I. e'qtto ar — iroTcra'a 

DH1 = M2-H1 

i F roMf, LT.O.) 60 TO 20 

If C A3S C DM1 ) .LT.TOL) GO TO 60 

i 6 “T 0~? 0 

so WRITEC6/1 ) 

H“s“170 

N a 1 .0 

— ro~Ti5~7n 

60 M a CHI + M2)*0.s 

FO^RETUPN 

80 URITEC6/A) 

GO' T0' 60 

1 FORMATC lOX^'**** ITERATION FOR H AND N D 1 VE RG ES * * * ' / 1 5 X ^ ' S E T 

X WZ'Z i , *ira'1T0'*7 / > 

A FORHATC 10Xi'****100 STEPS COMPLETE — DIO NOT CONVERGE***'//) 
R E T ur, N“ 
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^PENDIX B THE ROSIiAM/FimiAN PEDGRAM 


This appendix documents the Roskam/Pillman program’s use 
and provides a copy of the listing. Appendix B1 will define all 
input and output acronyms and explain how to prepare the ii^>ut 
(te.ta. Appendix B2 will provide a complete program listing. 
Appendix B3 provides an exanple output. 

B1 Preparation of Input Data 

Table B-1 defines all of the input and output ccaipxrter 
acronyms. 

Table B-2 presents the input card formats. All input cards 
must be in the order specified by Table B-2. Data for more than 
one baseline aircraft may be input by stacking data decks. For 
multiple data decks the last card of each set of data corresponds 
with the first card of the next set of data. The last card of 
the total data deck must be blank. This tells the program that 
the next set of data is a null set and terminates execution. 

Of particular interest are KL, K2, K3, K4, KEMPl and KEMP2. 
These cosrrespond with the weight correction factors of Section 
2. 3. 1.3. In this case an input value of zero will result in a 
correction factor of unity. An input value of one will result 
in the corresponding non-unity correction factor. 


B.l 


AEFF 

TIT’ 

Effective aspect ratio of 
the vertical tail (Reference 
14) 


AIECRAFT 

- 

Airci^t or configuration 
name 


AH 


Horizontal tail a^ect ratio 


ALT 

h 

Altitude 

ft 

AV 


Vertical tail as^ct ratio 


AW 

A 

Wing aspect ratio 


BHD 


Baseline horizontal tail 
span 

ft 

BVO 


Baseline vertical tail span 

ft 

BW 

b 

Wing span 

ft 

CBAR 

c 

Wing irsan aerodynamic chord 

ft 

CHO 


Baseline horj^ntal tail MAC 

ft 

QIA 

m 

a 

Static longitudinal 
stability 

1 — 1 
u 

a© 

\ 

Static directional 
stability 

rad ^ 

cvo 

c 

V 

Baseline vertical tail MAC 

ft 

DENSIT 

p 

Atmospheric density 

slugs/ ft 

DFUS 

D 

Fuselage equivalent 
outside diameter 

ft 





Table B-1 Roskamy 1 lllrt'iaii Program 


Acronym 

DIVE 

EP 

ETA 

FUSHGfT 

FDSWID 

HI 

H2 

HH 

IF 

K 

KL,K2,K3 

KAPAH 

KAPAV 

KAPAW 


Input/Output Acronyms (Cont'd) 


Variable 

Description 

IMits 


Dive speed 

Knots 

- 

Tail sizing span tolerance 
for iteration process 

ft 

n 

Dynamic pressure ratio at 



the Qipennage q-g/q^ 


H 

Fuselage outside height 

ft 

W 

Fuselage outside width 

ft 

^1 

Fuselage height at .25&g 

ft 


Fuselage height at 

ft 


Vertical distance from 
aircraft center line to 
horizontal tail root chord 
(positive down) 

ft 

- 

Number of different tail 
cones to be considered 


k 

Bipirical vertical tail 
factor fran Figure 7.3 of 
Reference 14 


” 

Weight estimation constant 
controls for fuselage (See 
text) 



Ratio of horizontal tail 
section lift slope to 2ir 


•'v 

Ratio of vertical tail 
section lift slope to 2ir 



Ratio of wing section lift 
slope to 2ir 



B.3 



Table B-1 Roskam/Filltnan Program 


Input /Output Acronyms (Cont'd) 


Acronym 

Variable 

Description 

Units 

KEIMPI 

% 

Ekipennage weight estimation 
consrtant control (See text) 


KEMP2 


Lk^nnage weight estimation 
constant control (See text) 


LBOID 


Baseline ibselage Length 

ft 

LCl 

- 

lliniraum tail cone length to 
be considered 


liCF 

- 

J/Iaxifflum tail cone length to 
be considered 


LOOLD 


Baseline tail cone length 

ft 

ISOLD 


Baseline hori 2 K>ntal tail 
moment arm 

ft 

LN 

^N1 

Nose length 

ft 

LU 

£ 

U 

Cabin length 

ft 

LVOLD 

V 

Baseline vertical tail 
moment arm 

ft 

MACH 

M 

Mach number 


PHICl 

1 — i 

o 

Tail cone shape parameters 


PHIC2 




PHINl 

^N1 

Nose cone shape parameters 


EHIN2 

'*’N2 



SBSCl 

- 

Side view projected area 
of the baseline tail cone 

ft2 
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Table B-1 Roskam/Flllmaii Program 


Input /Output Acronyms (Cbnt’d) 


Acronym 

Variable 

Description 

Units 

SBSN 


Side view projected area of 
the nose 

ft^ 

SBSU 

- 

Side view projected area of 
the cabin 

ft^ 

SBO 


Baseline horizontal tail 
area 

ft^ 

svo 


Baseline vertical tail area 

ft2 

sw 

S 

Reference wing area 

ft^ 

mm 

H 

Baseline horizontal tail 
leading edge aveep 

deg 

SWPVO 

- 

Baseline vertical tail 
leading edge sweep 

deg 

TAPRH 


Horizontal tail taper 
ratio 


TAPRV 

V 

Vertical tail taper ratio 


TAPRW 

X 

Wing taper ratio 


TAS 


True airspeed 

ft/sec 

TEST 

_ 

Test run number 


TEnCKH 


Horizontal tail thickness 
ratio 


THICK7 

(t/c)y 

Vertical tail thickness 
ratio 


TOL 

- 

Tolerance for the iteration 
in OGSISEIP (See Appendix A3) 


ULTDOAD 

'^ULT 

Ultimate load factor 

g's 
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Table B-1 Roskam/Fillmaa Program 



Input /Output Acronyms (Cont'd) 


Acronym 

Variable 

Description 

Units 

VISOOS 

u 

Kinematic viscosity 

ft^/sec 

VIGROSS 


Gross weight 

lbs 

\0«PLE 

^LE 

Wing Leading Edge Sweep 

deg 

XBARHO 

- 

Baseline horizontal tail 
aerod 3 ?namic center location 




as a fraction of c 


XBARW 


Wing aerodjmamic center 
location as a fraction 




of c 


XCBAR 

”” 

Location of the leading 

edge of c, relative to the 
nose 

ft 

XDG 

’'eg 

Baseline c.g. location 
relative to the nose 

ft 

XOffl 


Baseline horizontal tail 
c.g. relative to the nose 

ft 

XCGV 


Baseline vertical tail c.g. 
location relative to the 
nose 

ft 
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Table B-2 Roskam/Flllman Program Input Card Formats 


INPUT CARD NO. 1 - (13A6,A2) 

TEST AIRCRAFT 

I T - - " ' ' -| 

1 7 80 

INPUT CARD NO. 2 - (7F10.0) 

TAS MACH DENSIT VISCOS ALT WGROSS 

r — i — I I 'I r — I 

1 11 21 31 41 51 

INPUT CARD NO. 3 - (7F10. 

SW BW AW CBAR WSWPLE TAPRW XCBAR 

1 — r I I I I I I 

1 11 21 31 41 51 61 

INPUT CARD NO. 4 - (7F10.0) 

LBOLD DFUS LN LV LCOLD FUSWID FUSHGT 

I I I I I '• r I I 

1 11 21 31 41 51 61 



Table B-2 (Continued) 


INPUT CARD NO. 5-(7F10.0) 




HI 

H2 PH INI 

PH1N2 

PHICI 

PH1C2 TOL 

r ■■ 

1 1 _ 

1 _ 

1 _ 

1- 1 1 

1 

11 21 

31 

41 

51 61 

INPUT CARD NO. 6 - (412) 




K1 

K2 K3 

K4 



□ 

□ 

□ 

□ 



1 

3 5 

7 



INPUT CARD NO. 7-(7F10.0) 




SHO 

BHO AH 

SWPHO 

TAPRM 

THICKH 

L 

1 I . 

1 

1 _ 

l_ _ 1 

1 

11 21 

31 

41 

51 


INPUT CARD NO. 8-(7F10.0) 
LHOLD HH CHO 


31 


1 


11 


21 



1 


11 


Table B-2 (Continued) 



{7F10.0) 

SBSCl DIVE VLTLOAD 

I I 

21 31 41 
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Table B-2 (Continued) 


INPUT CARD NO. 13-(7F10.0) 
XC6 XC6H XC6V 


1 11 21 
INPUT CARD NO. 14-(7F10.0) 


XBARW 

XBARHO CMA 

CNB 


1 

1 "T 

1 

□ 

1 

11 21 

31 


INPUT CARD NO. 15 - (7F10.0) 



KAPAW 

KAPAH KAPAV ETA 


1 

1 


□ 

1 

11 21 

31 


INPUT CARD NO. 16-(3F10.0, 

IX. 13) 


LCl 

LCP IFl 

IF 


I _ 

1 . 

1 1 - -J 

1 

11 21 

3132 




' Table B-2 (Continued) 


INPUT CARD NO. 17 - SHOULD BE BLANK AS THE LAST CARD OF THE TOTAL DATA 

DECK. OTHERWISE. THIS CARD CORRESPONDS WITH 
THE INPUT CARD NO. 1 OF THE NEXT DATA SET. 


B2 Program Listing 

Hiis section provides a con^ilete listing of the Roskam/Fillman 
Program. 
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LABEL 
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1 DF HEW SI ON LCOD(500)-^LC 5 00) ^SVLCSOO) ,SHL (300 ) ,lCODO Cl )^L0(U#CD0(1 

I *)/C00L( 5U0WSV0 (1 (1 )/SHOt 1 ) / dHOn )^SMPllOCl ) / SWPVL *; 1 ) ^ S WE T F 0 ( 1 

3 »),SW£TEOn)#SljeTWa<1),SU£TCOn)r5U£TUOM)/WAFTfO<1)/UAlTEO(l)^BVL 

k * (5 00) ^OHLCSOO) / SWPHU5 00)^S WPVL ( 5 00 ) <•$ U F ( 500) # S WE ( 5 00 ) / S WC C 500 ) i^SW 

5 *N( 500 > #SWU (5 00 ) /WAITF C 500) ^ WA IT£ ( 500) # L 6L C 500) r XC G L < 5 0 0 ) OH C 5 00) 

o ^yCBVC 5t)0)^SWHC5 00)^SU\/£ 5 00) / SW T ( 500 ) #L ABELl < 1 4 ) /LABEL2 ( 1 4) # 

7 *LA0HL3n4)^ARCRFT(13) 

8 DF WEN SION XBvJ ( 500 ) / X 0D ( 5 00 ) , X8 WO ( 50 0) x X OH C 5 00) # XBA ( 500 ) 

9 COMMON TAS^AMACH/VI SCOS#SW^0W^ AWi-CBAR#WSWPLEj»T APRW#XCBAf)^L0OLD/ 

10 10FUS/LN^LU/LCOLD/FUSWlD^FU5HGT^Hl/H2^PHiNl#PHlN2rPHICWPHlC2/TOL# 

11 2Kl ^K2 ^K3/K4/SHO^BHO/AHi^5WPHO/TAPRH# THI CKH^LHOLDi^HH/SVO/BVO^AV# 

12 3SWP\/0 4»TAPRV^THICKV^LVOLOj^AEFF^K^K£HPWKEHP2#StiSN^S35Ui'SOSC1^0I VE# 

1 J AOLTLOAD^ XC G^ XC GH/ XCGV^ XDAR WB/ X BARHO^^ CM A / CMB^KAPAW/ KAPAH^ KAPAV^ ^ 

14 SETArDELTLC^EP 

15 COMMON/STAB/XBARW^XBAR ACVOXACBrXOARH 

16 COMHON/RHO/DENSI T 

1? COMMON/WEIGHT/WN^WU^WC^WF^ WHrWV/WEMP/CGERR/UGROSS/WGft ' 

18 REAL LCOD/L/LCOOO/LO^LBL/LN/LU/LCOLD/LHOLD# KAPAW^KAPAHi^ 

^^--[9 *LV0LD^LB01D/KAPAV#LC,L/LC2^ LQNEW^LCNEW ' ‘ 

20 REAL LC1#LCF/IF1 

21 DATA (LABEL! C I ) ^ 1 = WS) /6HLCOO V^6HS, CD0/6H -'SCA^aHLEO VA^ 

22 lAHLUESy 

23 DATA (LABEL2(1)/I = W5)/6HLC VS,^6H SH - ^6HSCALED^6K VALUE^IHS/ 

24 data (LA0EU3( I)/ I-l /S )/6HLC VS-^6H $V - ^6H$CALED^6H VALUE^lHS/ 

25 ■ " DATA TERM/6H / ‘ 

26 10 READ(5#1 )TE$T#(ARCftFT( 3) 

27 IF CTEST,EQ*TERM) STOP ' ^ 

2 8 RE ADt 5/2>T AS/AMACH^OENS IT/ VlSCOS.rALT#WGROSS 

29 READC5/2)SW/BW/AW/CaAR/USWPLE/TAPRW/XC0AR 

30 RE ADC 5/2)L0Ol,D/DFUS/LN/LU/LCOLD/FUSWlD/ FU5HGT 

READC5/2)H1/H2/PHIN1/PFUN2/PHIC1/PHIC2/T0L ' 

32 READ(5/3)K1/K2/K3/K4 

33 READ(5/2)SM0/0M0/AH/SUPH0/TAPRH/THICKH 

34 READ(5/2)LHOLD/HH/CHO 

35 REA0C5/2)SV0/3V0/AV/SUPV0/TAPRV/THICKV 

36 REAP(5/2)LV0L0/AEFF/K/EP/CV0 

_ . READ ( 5V3)KEMP1 /KEHP2 ‘ 

3 8 RE ADC 5/2)S0SN/S0SU/SeS Cl /D I VE/ UL TLOAO 

39 REA9C5/2)XCG/XCGH/XCGV “ ’ 

40 READC 5/2)XBARW /XB ARHO/ CHA/ CN0 

41 READ(5/2)KAPAIJ/KAPAH/KAPAV/ETA ' ' 

42 READC54r4)LCl/LCF/IFl/IF 

OLC^CLCF-LCl) /in ‘ 

44 DELTC = 0, 

45 call SUETCPHIN1/PHIC1 / SHO/ S VO/ LCOLD /SWETNO# SWETUO/SWET £0/ 

46 *SWETH0/SW£TV0) 

4 7 SUET F 0 = SWE TNO + SWETUO + S WETC 0 

48 SWETEO=SUCTHO + SWE TVt) 

A9 “ ■ CALL CDO CL0OLO/ SUPV/SWPH>3WETH0 /SWETVO/SWETFO/CHO/ rVO# CDOBHV) 

50 CALL FLfSWGTC S METFO/SWE TN 0/ S WE T U 0 / S WETC 0 /Wf / WN/ UlU/ W C) 

51 CALL EttPWGTC SWETHO/SWETVO/ S WET EO/ SHO/GHO/ SV 0/8VO/ S WPH# S WPV# WV/ WH/ 

52 *WEMP) 
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LAOEL 


PAGE 


2 


C- 




55" 

54 

55 

56 

57 

58 
59" 
60 
61 
62 
63 
6A 

^6r 

66 

67 

68 

69 

70 

‘?r 

72 

73 

74 

75 


‘CALL~"XBARCG(0> j^BCG) 

SWPW = WSWPLE * 57.3 
WftITE(6^5l TEST^rtRCRTT 

Wn ITE(6#6) ALT^TAS^ AMACH/DEHblT^VISCOS 

WRITE (6/7) SW^BW^ AW/CBAR/XC 8AR/ S UPW/TAPR W 

WR ITE (6/8) LBOLD/LN/LU/LCOLO/FUSWIO/ FUSHGT/DFUS 

"WRlTEC6/9>SH0/BH0/AH/CH0/TAPRM/rHlCKH 

WRITE (6/11 )SVO/BVO/AV/CVO/TAPHV/TH1CK\/ 

WR JTE C6^ 1 2 ) TEST/ARCRf T/WGROSS 

WR ITE (6/l'3)tJN/kU/WC/WF/WH/WV/WEnP/XaCG 

WR ITE (6/14)TEST/ARCRFr 

WRITE (6 / 15) SWE_TNO,SWETU^S WEJCO/SW^fO/SWETHO/SWETyO/S W£T£0 
"SWPHa'SWPHO 
SWPV=SWPVO 
SW£T F=SWETFO' 

SWETN=SWETNO 

‘SMETEsSWETEO 

SWETC=?SWETCO 

SWETU = SweTUO ■" ' ■■ 

SH = SHO 
BH = BHO 

SV = SVO _ _ 

BV a DVO ' ■ ^ 


76 

CH=CHO 


77 

CV=CUO“ “ 


78 

LCODO n } 

- LCOLO/DfUS 

79 

"COOd) = 

COOBHV 

80 

SVQM ) = 

SVO 

61 

■ SHUC1 ) = 

SHO 

82 

BVOd) = 

0VO 

83 

BHOd V = ■ 

BHO' 

84 

SMPHQ (1 ) 

SWPH * 57. 

85 

■ ■ SVIPVO (1) 

= SWPV * 57. 

86 

SUETFOd ) 

= SWETFO 

87 

■■■ SWETEOd ) 

= SWETEO 

88 

SUETNQd ) 

— /* . fe n y ^ 

= SWETNO 



WA IT F0< 1 ) = WF 
WAITEO( 1 ) = WEMP 
LOCI) = LCOLD 
DO 100 1=^1 /IF 
LC2=LC1+(CI-1>*DLC) 

LBNEW = LMfLU + LC2 

OELTLC=LC2-LCOLD 

CALL SWETCPHINI jrpHlC1/SH#SV/LC2/3WeTN/S WETB^SWETC/SWETH/SWETV) 
CALL FUSWGT (SWETF/SWETN/SMETU/ S WETC/WF/ WN/WU/WC ) 

CALL EMPUGK SWETH^SWETtf^SMETE/ SH/BH/SV/GV/S WPH/SWPV/WV/WH/WEMP) 
CALL X8ARCGCI/XBCG) 

I 'MAY NOT BE REDEFINED IN CALL OR ABNORMAL FUNCTION 
CALL ST ABASE A( XBCG/HSWPL E/SH/8H/CH) 

CALL VERTAH£A( VSWPLE/S V#GV#CV> 

CALL SWET(PHINl /PHIC 1/SH/SV/LC2/SWETN/SWETU/SWETC/ SWETH#SWETV> 
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PAGE 


104 SWEtr’a" SWET«‘'f'’SWETU SW^TC 

105 SWETE a SWETH + SWETV 

1G6 X.ALL C00(‘L8NEW/VSWPLE/HSWPLE#SWETH#SWETV#SW£Tf^CH#CV#C00aHV} 

107 SWV(I) a SWETV 

lOS ' ‘ SWHU ) a SWETH 

109 SMTUJ = SWETf + SWETE 

110 LCODCn = LC2/DfUS' " 

111 UI) = LC2 

112 X0WCI) = X0ARW 

115 XB8(I> a DXACO 

IK ^ "X0WBU) = XBARWB . .. . 

115 xeiun c xoAftH 

116 X8An> is'XBAJiAC ^ 

11? SVL(n a sv 

“118 ■“ “ '^’ SHL(I) a SH " 


td 

Ol 


119 COOt(l) = COO0HV 

■ 120 8 VL< I ) = ev ■ 

121 0HLC1) = BH 

122 SUPHLCIJ = HSWPLE * 57. 3 

123 SUPVL(I) = USWPLE » 57.3 

124 ■ SWf <1 > = SWETF 

125 SME(I) = SMETE 

-126 ^ SMCdr* SWETC ■ 

127 SWNCI) = SWETN 

12B SWU(0" = SWETU 

129 UAITF Cl) = WG« + WC 

130 WAITE (I) = UE«P“ ■ 

131 L0L(I> » LBNEM 

132 ~ XC6CCn = XBCG 

133 CBHCI) = CH 

134 CBVm = CV — 

135 1Q0 CONTtNUE 

■ 136 - . WBirE(6/16)TEST/AnCRFT,UGfi#(a,La),LBLCnrWAITfa)/WAITea)» 

137 *XCGL( I))^I=1^IF) 

" ISa WRITE (6#17)TEST»ARCRFT,CMA /AH#TAPRH 

139 WRITE (6^18) ( C I #L U )>SML f I) # BHI. ( I ) <■ CQH( I)#SMPHU I) ) ^ 1=1 / I f) 

Hi 5A(2/L#SHL/!F,L0ySH0^1,LA8£L2) “ 

141 WRITE (6» 19) TEST,ARCRFT/CMl}#AV»TAPfiV 

142 WHITE C6#18) ((I/LCn/SVL(IJ#BVL Cl )^^CBV( I )#S«PUL( D) »I=1 #IF) 

143 CALL H15AC2,L/SVL#IF,L0,SV0#1, LABELS) 

- -TAA WRITEC6,2S) T EST# AfiCRF I # CH A 

145 WR ITE C6#26) CCI#LC I)» X0WCI) /XBB Cl ) «XaWB C I) rXBHC I ) /X :AC I ) #XCGLC I ) ) f 

146 *I = 1^IF) - - . - 

147 WRITE C6, 21 )T£ST,ARCRFt 

148 ■ ■ WRITEC6<^22)CCI#LCI)»SWF(I),SWHtI)#SUVCI)#SWlCI)),Ia1,IF) 

149 WRITE (6«23)TEST«ARCRFTfCD0t1) 

150 - WRITEC6»24)<CI,LCl)/LC0DCI),C0nLU)),I=UIF) 

151 CALL H15AC2,LCOO/^COOL/ IF#LCOOO,CDO#1/LABEL1 ) 

152 GOTO 10 ■ — - - 

153 1 fORHATC13A6»A2) 

154 2 FORMATC7F10.0) 

155 3 FDHHATC4I2) 


I 


i 


d 


c 
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156 " 6 f0finATt3F1CI.0r1X,15)' '■ 

157 5 F0HMAT(1H1>10X, ‘TABLE ‘i.A6/'.1A ROS KA(1- F I L LMAN HE T M 0 D ’ / 37X , 1 2 A6# 

158 *A2/27x» ‘ORIGIMAL FLIGHT CONDITIONS AND • /27X , ‘ 6E0METR K DEFINITIONS 

159 *’//10X,65(' *' ) //) 

160 6 FOBMATt lOX^ ' FL IGHT C ON 0 1 T I ONS ' // 1 5 X , * A L T I TU D E ' , T50 y F 1 0 . 3 <■ * FT'/ 

161 *15X#‘ AIRSPEED*, TSO/^FIO. 3/' FPS ' / 1 5 X, ' H A CH • ,T 5 0 , F 1 0 . 3/ 1 5 X# ' OEM S I T Y 

162 i';T50VFl0.5,‘ SLUGS/CU.FT. • /15x;‘KINEHATlC V I S C OS I T Y ' , T 5 0, F 1 0. 5# 

163 •• SO.FT./SEC'//> 

■ 16A ■ 7 FOfiMATdOX/'GEOHETRIC* //15X,‘UING'/20X, 'REFERENCE AR EA ' , T5 0, F 10 .3, 

165 *' SU.FT.'/20X,*SPAN',T5D/F10.3/' FT * /2 OX, ' A SPE C T R AT 10 * , TS 0, Fl 0. 3/ 

166 *20X, 'H, A.C. ',T50, F10.3, ' F T ' / 2 2 X , ■ F . S. ' , T 50 , F 1 0. 3, ' FT ’ / 20X , ‘ L. E. 

167 *SWEEP',T50,f10.3,' 0 E G ‘ / 20 X , ' T APER ' / T5 0 , F 1 0 . 3// ) 

168 6 F6RMAf(i5X,'FUSELAG£'/20X,'0V£RALL “LENGTH'/tSO/FlU.S,' d'/20X, 

169 *'NOSE LENGTH' ,T5Q,F1C. 3, • F T ' / 2 OX , • CAB I N L E NGT H ’ , T 50 , F 1 0. 3 , ' FT*/ 

170 *20X,*TAIL CONE LENGTH* ,T50, f 10 . 3, ' FT'/20X,'MAX M I OT H ', T 50 , FI 0. 3, 

171 *' FT*/20X,*HAX HEIGHT', TS0,F1O. 3,' FT ' / 20X , ' EQ U I V. 0 1 AM ETE R * , T50, 

■ ■ '172 ■ ■ *F10,3,* FT*//) 

173 9 F0RHAin5X, 'HORIZONTAL T AI L ' /20X , ■ R EFEH ENCE Afi E A ' , T 50, F 1 0. 3, 

**SQ. FT. ' /20X,'*SPAN',TS0,F10.3, *‘ F T ' /20 X , ' A S PECT ' R AT 10 ' , T 50 , FI 0. 3/ ' 
175 •20X,*M.A.C.‘,T50,F10.3,' FT'/ 

■ 176 ■ *20X,*TAPER‘ ,T50,F10.3/20X, * THI CKNESS‘.-T50,F10.3//) 

177 11 FORHATdSX,' VERTICAL T A I L ' / 2DX , ' RE F ERE N CE A REA ' , T5 0, FI 0 . 3, 

■ “ 178 ■ ♦•SQ.FT.*/20X,'SPAN',TSO,F10.3,' FT • /2Q X , ' AS PEC T R A T 1 0 ’ , T 50, FI 0, 3/ 

179 *20X, ' H.A.C. ' ,T50,F10.3, ' FT'/ 

*20X, 'TAPER', T50,F10.3/2DX,‘ THICKNESS*, l5n,F10. 3//) 

181 _ 12 FORHATdHl,10X, 'TABLE ',A6,'.1B RO S KAH - F I i_l MAN MET HOD ' / 2 7 X , 1 2 A6, 

182 «A2/27X, 'ORIGINAL WEIGHT AND BALANCE ST A T £ Mt NT • / / 1 OX , 65 C ' * ' ) / / 

183 *T52, *{LBS.)'//10X, 'GROSS WEIGHT', TS5/,F10. 3//) 

W 184 “ ' ■ “"13 FORMAT! 10X, ' APPROXIMAT ED SHELL WEIGHTS’, " //15X, 

^ 185 ** FUSE LAGE* //20X, * NOS E S E CT I ON ' , T 5 0, F 1 0 , 3/20 X, ' C AB I N S E C T 1 0 N ' , T5 0, 

^ — 1B6T *F10.3/20X, ' TAIL CONE S E C Tl ON ' , T 5 0, F 10 . 3 // 20 X , ' T OT AL FUS E LAGE ' ,T 50, 

1 87 • no. 3 // 1 5X, 'EMPENNAGE* // 20 X,' HORIZONTAL TAIL',T50,F10.3/20X, 

188 " " “'VERTICAL TAIL', T5D,F10.3//20X, 'TOTAL EMPE NN AG E ' , 1 50, F 1 0 -3 / /1 0 X , 

189 •'CENTER OF GRAVITY — FRACTION OF H . A. C . ' , T 50, F 1 0 . 3 / / ) 

'190 ' 14 FORMATdHl,10X, 'TABLE ' , A6 , ' , 1 C ROSKAH- F I LL M AN ME T HOD • / 2 7X , 1 2 A6, 

191 • A2/2 7X, * OR IGINAL WETTED AREA □ fl E AK D OWN ' / / 1 0 X , 6 5 t * * ' ) / / ) 

* "192 ir FORMATdOX, • FUSELAGE ', T51, • (SQ. FT. )'// 1 5X, ' NOSE S E C T 1 0 N ' , T 5 0, F 10.3 

193 »/1SX, 'CABIN SECTION', T50,F10.3/15X, 'TAIL S £ C T I ON ' , T 50, FI 0 . 3 // 

194 »15X, 'TOTAL FU S EL A GE * , T 5 0 , F 1 0 ,3 / / 1 0 X , ' E M PENN AG E ' / / 1 5 X , ' H 0 R I Z ONT AL T 

195 *AIL',TSO,F10.3,/1 5X, 'VERTI CAL T A 1 L ' , TS 0 , F 1 0 . 3 / / 1 5 X , ’ T OT AL EHPENNAG 

196 •£' ,T50,F 10.3//) 

197 16 F0RMATdH1,10X, 'TABLE ',A6,'.2A ROSKAH- F I LL MAM MF T HOD ' / 27X , 1 2A6, 

*A2/27X, 'OPTIMIZATION WEIGHT AND Q AL ANC £ ' / 27 X , ‘ S T A F EHEN T * / / 1 OX , 6 5 ( * 

199 »•• )//10X, 'GROSS HEIGHT! W/0 TAIL OH EH P . ) ' , T 5 5 , F 1 0 . 3 , ' CLBS.)' 

200 • //Tl 2,' I ', 118 , 'lC, T31, 'LB' , T41, ' FUSELAGE *,T54, • EMPENNAGE', T69, 

201 «'C.G.'/T44,'Wt’,T5e,'WT'/T17,'(FT)',T30,'CFT)',T43,'(LB)',T57, 

202 •• !LQ) *,T67,' !ffl C BAR > ' / / ! 1 OX, I 3 , 5 C F 10 . 3, 3X ) / ) ) 

203 17 FORMAT! llll ,10X, • TABLE *,A6,'.2B RO S K AH - F 1 LL MAN HE T HOD ' / 2 7X , 1 2 A6, 

.A2/27X, ' OPTIMIZATION S TABl LITER ' /27X, ' SIZING' //27X, ' CHA = ',T55, ' 

205 »F10.3/27X, 'ASPECT RATIO = ' , T 5 5 , F 1 0 • 3/ 2 7 X , ' T AP E R RATIO = ',T55, 

206 •F10.3//10X,6S!'»')//T12,'I*,T1B,'LC',T31,'SH',T44,'BH',i57,'CBARH‘ 

207 *,T68, 'SWEEPH* /T17,'! FT) ' ,T 2 8, ' C S Q . FT. ) • , T4 3 , ' C F T) ' , T58, ' C F T ) ' ,T 69, 


^ - ; ; ; ; ; | ‘ . . • , . e . . c 
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20a *‘CDBG)'//) ■■■ 

209 18 F0fiMAT((10X/J3/5CF10.3^!x) /)) 

210 19 FOBMAK 1 HI #10X, ' TABLE *^A6,’.2C HOS RAM- F 1 LL MAM HE TIIOO • /27X , 1 2 A6# 

211 *A2/27X#‘0PTIHIZAT10N VERTICAL T A I L ' 727 X / • S 1 ZING * //2 7x, • CNB = %TSS 

212 •,F10.3/27X/'ASPECT RATIO = ’ ,T 5 5 » F 1 0 . 3/2 7X^ • TAPE R RATIO =*,T55» 

213 *f10.3//10x.65('*')//T12#>I',I18,'LC%r31»*SV’,TA4,'ev*,T57,*C0AfiV 

214 *,T6B,^ 'SUEEPV' /TIZ^'CFT) ■ ,T2 8, ■ C S Q . FT .) * 3 y ' I f T1 • ,T58r ' C FT ) • ,T69, 

215 ‘‘(DEG)'//) 

• “216 21 FORHATHHI ylQXy 'TABLE ' , A6y ’ . 3 A ROS KAM- FI LLMAN METHOD ' /27X y 1 2A6y 

217 *A2/27X^ 'SENSITI VI TV OF WETTED AR E A • /27 X , • T 0 TAIL CONE L£N6TH'//1 X 

218 “ ■ *,65( ■»') //Tl 2y ’ I ' yT20y ' LC yT29y ■ FUSELAGE' yT42y* HOR IZ',T55y 'VERT* 

219 *yT67y 'TOTAL '/T19y*(FI)'yT29y‘fSQ.FT.)'yT41y'(SQ.FT.)'yT53y'(SQ.Fr. 

220 -)',T65y'CSQ.FT.)'//) ^ 

221 22 FORNAT t ( lOXy 13 y2Xy 5( FI 0. 3/2X)/) ) 

■ 222 23 FORHATdHlylOXy'TABLE 'yA6y*.3B fiOS KAH- F I LLMAN HE T HOO * /2 7X y 1 2 A6 y 

223 *A2 /27Xy ' SENS ITI VI TY OF Z £R 0 -L I F T • y / 2 7X y ' OR A 6 TO TAIL CONE LENGTH* 

224 ■■ »//27Xy'0filGINAL C DOBH V ' y TS Oy F 1 0 . S 

225 *//10Xy65C *') //T1 2y' I' , T 23 y ' L C ' y T 35 y 'L C/0 ' y T4 7y ' CD08HW ' /T22y •( FT) • 

226 *//) - -- - - 

227 24 FORMAT ( ( 1 1 X y I 3 y 2 X y 2 ( F 1 0 .3 y 3X ) y F 1 0. 5/) ) 

... 25 FORHATdHlylOXy’TABLE 'yA6y'.20 ROSKAM-FI LLMAN HE THOO » /27x y 1 2 A6 y 

229 *A2y/27Xy'STATIC MARGIN C OH P ONE NT S ' / /27 X y • CHA = ' yT5 Oy FI 0 . 3/ / 

- - 230 *10Xy65( '*•) //) 

231 26 FORMATd 2Xy * I * ySXy 'LC ' y6Xy ' XQW'y4Xy 'XBB 'y4Xy 'XBWB* y2Xy ' XBH'y4Xy 

- — 232 v»XBA'y4Xy'XQCG'//d0XyI3yF10.3y6F7.3/)) ' 

233 END 

1470 EOUALTTVOR NON-EQUALITY COMPARISON HAY NOT BE MEANINGFUL IN LOGICAL IF EXPRESSIONS 
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1 “ subroutine SWEt {PHtN.?PHlC/SM^SV#TALUGT^SWeTNrSweTU/SW£TC# 

2 ^ *SWETH/SWETV> 

3 ^ ^ - COMMON TAS/AMACH4rVISC0S#SW#flW#AW#'CDAft, WSwPLE/TAPRW#XCBARxLBOLD# 

4 IDFUS^ NOSLGT^CADLGT^LCOLOxFUSWlD^FUSHGr^Hl ^ H 2 # PH I Nl / PH I N 2^P H I C 1 # 

5 * PHIC2# TOL^ 

6 2K1 ^K2^K5yK4#SHO/BHO^AH^SWPHO^TAPRH/THKKH^LHOLO^HH ^SVO^BVO^AV/ 

'7 3SUPVO/TAPRV#THICKV^LVOLO/AEF F^K/KEMP1#KEMP2/SSSN^SBSU^SBSC 1 ^DI VE/ 

B 4ULTL0AD#XCG^ XC GH^ X CGV^ XB AR W X B A RHO ^CHA / CNS ^ K AP AW^ K AP AH *-KA PAV^ 

_ 9 ■' ■ SETA^DELTLC^EP 

10 REAL KW/KV/KC/KA/NOSLGT 

■ 11 DATA PW3. 14159/ - . 

12 CIRCUM = OFUS * PI 

13 SWETU = CtfiCUM CABL6T 

14 call C0NPAR(PHIN<KA^KV#KC#KW) 

.... SWETN = CIRCUM * NOSLGT * KW 

16 CALL CONPARCPHfC^KA^KV^KC^KW) 

17 ■ ■’ SUETC = CIRCUM * TALLGT * KW 

18 CKP = .52 * THICKH 1,987 

19 5METH"- CKP * SH 

20 CKP = .52 ♦ THICKV + 1.987 

21 SWETV = CKP * $V 

22 RETURN 

“ '23 END ^ 
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LAUtL COO 


PAtiE 


1 


cct 

* 

^3 



2 
■ 3 
A 

5 

6 

~ 7 * 

a 

9 

10 

11 

12 

‘1 3 “ 

14 

15 

16 

17 

18 
" 19 ' 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 - 

32 

33 

34 

35 

36 

37 

38 

39 

40 


CDOFUS 
R E NU M H 


650 

'660 ■” 


6 70 
680 — 


”'SlJflflOUnMe"COiJ'fAL0^SWPVT;SWPHT#$weTHr/SWETVTy^SweTf /CBAfiH/CDARV/ 

1 CDOBH V) 

■ COMMON TAS#A«ACH^ VISCOS^SW^DU^AWrCBAR# WSWPLE^T APRw^JCCBAfi rLBOLD^ 

1 OFOS ^LN^LiJ^LCOLD^f USWJ 0 / FUS HGT / H WH 2 ^PHIN 1 ^PHIN 2 /PHIC 1 #PH 1 C 2 ^T 0 L^ 
2 kT /K 2 ^K 5 /K 4 ^SH 0 -'BH 0 ^ AH/SUPHO/TAPRH^THlCKIi#LHOLD/HH/SV 0 ^ 090 ^AV^ 
3 SWPV 0 ^TAPRV^THiacV^LU 0 L 0 ^A£FF^KxKEHPUKEMP 2 /SaSN^SBSU^ SGSCI^DI V£^ 
" 4 ULTL DAD^XC G^XCGJUXCGV^ XfiAft tJQxXHARJtO^CMA^CNB^KAPAW^K APAH/*CAPAV^ 
SETA, D FLTLC,EP 

RENUM? = RENUH(DE«SIT,VISCOS^TAS,ALB> 

CFnO = 0.4 55 / CAL 0610 (RENUMF) ** 2.5 8 ) 

C FCM = 6 . 899 * AMACH- 1 6 . 226 * AH ACH* *2 + 1 5 . 74 1 *AMACH + * 3 - 5 • 4 89 4 * AH A CH * * 4 
SKIN - CFHO *CFCH 

“FINE =”ALB A" Of US 

SKIN*SVIETF/SU *<1 + 6 0 /FINE** 3 +. 0025 *FINE) 

RENUHCDENS I T, VI SCOS^TAS,C 0 ARH) 

CFHO = . 455 /<AL 0610 (REWUMH) ** 2 . 58 ) 

SKIN = CFMO *CFCH 

HSWPC 3 = ATAnC (SIN tSWPMT > /COS <SNPHT) ) -C 4 ,/AH* C . 3 3 3 * C 1 -T AP RH) / 

*(l . + TAPRH) ))) 

VSWPC 3 = ATAN((S 1 N(SWPVT)/C 0 SCSWPVT))-( 4 ./AV* <. 333 *tl.-TAPRV)/ 

*n . + TAPRV) )) > 

SWPH = COS ChS'^PC 3 ) 

'SWPV = C 0 S(VSVJPC 3 ) 

IF CAHACH.lt, 0 . 25 ) GOTO 650 

^Y = 13 . 36 - 79 . 39 *AMACH + 2 l 2 . 8 *AnACH** 2 - 232 . 9 *AMACH** 3 + 94 , 3 *AHACH**‘ 4 ' 
RLlFT = y- 16 , 12 *SWPH + 3 3 . 1 5 *S UPH* * 2 - 2 7 . 82 * SWPH * * 3 + 8 . 3 2 *SW PK ** 4 
GOTO 660 

RLIFT= 3.54 1 - 16 . 12 *SWPH + 3 3.1 5 * SWPH* * 2-2 7 , 82 * SWPH **3 + 8 . 32 *SWPII* + 4 
■ CDOHT = SKIN* SWETHT/SW*RLIFT * (1 + 1 . 2 *THIC KH + 1 00 *T H ICKH ** 4 ) 

RENUnV = REWUM(DENSIT,VISCOS^TAS,CBARV> 

- CFHO =' 0 . 455 / CAL 0 G 10 CRENUH V)** 2 . 58 ) 

SKIN = CFHO *CFCH 
IF(AHACH.LT.D. 25 ) GOTO 670 

RL IFT=Y- 16 . 1 2 *SWPV + 33 - 1 5 *S WPV * * 2 - 2 7. 82 * S WP V **3 + 8 , 3 2 *S WP V * * 4 
GOTO 680 

RLIFT= 3 - 541 - 16 . 12 +SWPV+ 33 . 1 5 *S WPV** 2 - 27 . 82 *SWPV** 3 +B. 32 *SWPV **4 

CDOVT = SK IN*RLI FT*$WETVT/SW* n + 1 . 2 *TH I C K V + 1 00 * T H ICK V * + 4 ) 

COO 0 HV = CDOFUS + CDOHT + COOVT 

RETURN 

END 


s 
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LABEL fUSWGT 


U# 


1 rUBRbUTlNE FjJSWGT(SWErF;SW£TN/SHETU;SljETC-^UF^UN/ljU/WC) 

2 COHHOM TAS,r AMACH^VISCOS/SW^BW/AW^CBAR/ WSMPL£^TAPRW#XCaAR,UQOLD# 

3 l&f US#LN/IU/LC0LD,FUSWID^FUSHGT^H1 ,h2^PHIM1^PH1N2/PHIC1 ^PHIC2^T0L^ 

4 2K1 #<2^K3^K4^SHO^BHO^AN^SUPHO#TAPRH^THlCKH/LHOLD^HH^SVO/BVOi'AV# 

5 3SWPVG^TAPRV^THICKV^LU0LD^AEfF^K/KertPWKEMP2^S0SN^S8SU/SQ5Cl #DIVe^ 

6 4UL TLOAD^XCG^XCGH^XCGV^ XB AR X 8 A RNO A ^ CNB # K AP AW^ K APAH/ KA PAV^ 

7 5£TA^DElTlC^EP“~ ' ^ ' 

8 REAL LHOLD 

9 C 

10 C THIS SUBROUTINE^ FUSWGTi^ ESTIMATES THE FUSELAGE WEIGHT WHEN HTE 

n C DIVE SPEEDS WETTED AREA/ FUSELAGE GEOriETRYi^ AND FOUR CORRECTION 

>2 C FACTORS TO ACCOUNT FOR PRESSURIZATION/ ENGINE LOCATION/ AND 

15^ ”C landing gear LOCATION ARE KNOWN 


14 C 

15 ■■■■■ C 

16 C DETERMINATION OF FUSELAGE CORRECTION FACTORS FROM OPTIONCODES 

17 C 

18 If (Kl-1 >130/120/130 

19 120 FUS<1=1*DB 

20 G0T0140 

21 130 FU SKI =1.00 

22 UQ IFCK2-D 160/150/1 60 

23 ■ " 150 FUSK2=1 .04 

24 G0T0170 

25 160 ‘FUSK2 = 1, ' ' 

26 170 IF (K3-D190/180/190 

“ 27 180 FUSK3nl,07 " " 

28 GOT0200 

- 29 19D ^ FUS<3=1 . 



30 

20G 

IF<k4-1 >220/210/220 

, 

31 

210 “ 

FUSK4n0.96 

o 

32 


G0T0250 


33 220 ’ FUSK4=1 . 

34 C 

55 C FROM THE INPUT DATA AND THE FUSELAGE CORRECTION FACTORS/ THE 

36 C ESTIMATED FUSELAGE WEIGHT CAN 0£ COMPUTED 

38 230 FUSK = f USK1*FUSK2*FUSK3*FUSK4 

39 ■ OCWHT = LHOLD + DBLTLC 

40 A = DIVE * OCWHT 

41 B = FUSWID fUSHGT 

42 C = A/e 

‘43 D“= SMETF**r;2 

44 £ = SGRTCC) 

45 WF = 0.021*FUSK*0*£ 

4i AFUSK = WF / SWETF 

47 WN = AFUSK • SWETW 

48 WU = AFUSK * SWETU 

WC = AFUSK *‘5WETC^ - ’ 

50 RETURN 

51 END 
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LABEL EMPMGT PAGE 


1 


* _ 

N5 


i 

3 

A 

5 

6 

7 

8 
9 

10 
11 
12 
^13" 
1 A 
15 
1 6 
17 

ia 

19‘ 

20 

21 

22 

23 

2A 

-25- 

26 

27- 

2B 

29 

30 
'31- 

32 

33 
3A 

35 

36 

37 

38 

39 
AO 
A1 
A2 
43^ 
AA 
AS 
A6 
A7 
AS 

49- 

50 


C 
C 
C ■ 

c 

~c 

500 

520 

530 

5A0 

“550- 

560 

C 

c 

c 

-570- 


” SUdftOUTlNE EMPWGT (SWETtUSWETV^ SWETE/HT AftEA/HTSPAN/VTAJlEA^yTSPAN# ^ 
*HS UP LE#VSWPL£/VTWGT^HTWGT/ EHPWGH) 

COMMON TAS, AHACH^ VI SCOS^ SIJ# OWr AU,CDARr MSWPLE^TAPRW^XCOAR^LBOLD^r 
1 OPUS# LNxLU^ LCOLO^f USWI D^fUSHGT / H2/PH INI rPHIW2^PHlC1/r'HIC2/T0L^ 

2Kl ^K2 #K3 /KA #SHO^[3HO^AM#SWPHO^T APfiH, THI CKH^LHOLD#HH ^SVO^DVO^AV# 
3SWPV0^TAPRV-*TinCKV^LU0LD^A£FP^K>K£MPl,KEMP2,SBSW^SGSU^SBSC1 ^DIVEjp 
AUL TLOAD^KC G/XCGH/XCGV, XB AR W8 B A RHOi' CM A # CNB ^ K A P A W/ K AP A H / K AP A V # 
SETA^DELTLC^EP 

THIS SUBROUTINE/ EHPWGT/ CALCULATES THE ESTIMATED WEIGHT Of AN 
EMPENNAGE WHEN THE DIVE SPEED/ WETTED AREA/ AND EMPENNAGE 
GEOMETR Y ARE KNOWN 




EMPENNAGE WEIGHT FOR DIVE SPEED LESS THAN 250 KTS 

IF (D I VE,GT,25 0) GOTO'SOD 
EMPWGH=.rA*(ULIL0AD*SWETE*^2)** ,75 
HTWGT =: i,MPWGH*SWE7H/SWETE 

VTWGT =■ EHPWGHiSWETV/S WETE ' ' ' 

GOTO 630 

EMPENNAGE WEIGHT WHEN DIVE SPEED GREATER THAN 2V0 KTS 
DETERMINATION OP EMPENNAGE CORRECTION FACTORS 


IF (KEMPl -1 >530/520/530 

EMPKHat.l ^ - - - .... - 

G0TO540 

EMPKH=1 . . - - ^ - 

IF CKEMP2-1 ) 560/550/560 

•EMPKV=1 .+,15* (HTAREA*HTSPAN/VTAREA7VT5PAN> 

G0T0570 

EMPKVsl, 

DETERMINATION OF HORIZONTAL AND VERTICAT TAIL CORRECTION FACTORS 


580 

-50- 

600 

6TQ 

63J 


— HTCF-HTAREA**.2*DIVE/1 000,/ (saRT(COS(HSWPLE)>)' 
VTCF=VTAREA**.2*DIVE/1 000- / (SQ RT ( COS< V S WPL E > ) ) 
If (HTCf.GT,! .40)6010580 
KTWCT=:( 3 -9*HTC F-,419)* HTAREA^EMPKH 
G0T0590 

HTwGT“t 5.?4-2-2 9*ALOG CHTCF ) ) *HTAREA*£rtPKH 
-IFCVTCF-GT.r.40)GOT060D 


VTW6T=(3,9*VTCF-.419>*VTAREA*EMPKV 
VTW6T=( 5.74-2,29 *AL06<VTCF) ) * VT A R EA*EM P KV 
G0T0610 

EHPWGH = HTWGT + VTWGT 
CONTI NUE 

■ RETURN 

END 


a 

# 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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LABEL X8ARCG PAGE 1 


tJi 

* 

]SJ 


1 STJenOUflNE' X&ABC6(l#)(BCer" ' ■ 

2 COMMON TASf AMACH,VISCOS/SW^dU«'AU^CBAK/USUPLE«TAPRU<XCeAR#LBQLO/ 

3 10FUS#LH/LU/LCOLO,FUSUIO/fUS«fiT,Hl/H2/PHIN1#PHIN2-PHIC1#PHIC2rTOL/ 

4 2K1,K2 J.K3/K4 ,SHO/BHO,AH»SUPHO»rAPRH#THI CKH/L HOLD»!IH/SVO/ BVO/AV/ 

" ' 5 T 3SWPVO/TAP(JV,rHICKV,LVOLO/AEFF<^K/KEMP1/KEMP2,SBSN^SBSU,SBSC1#OIVE» 

6 4UL TL0AD#XC6»XCGH»^XCGV#XBARW0»XBAnH0,CHA/CNB#KAPAW/KAPAH»KAPAU^ 

y 5ETA#DELTLC,EP 

a COMMON / UEIGHT/WN«UU/UC«UF/UK#UV<^U£MP^CGERR#UGROSSaUGR 

■ 9 COHMON/CEEGEE/XCGI 

10 COMMON/SHAPE/HN^NN^HC/NC 

11 C ■ " 

12 C this subroutine ESTIMATES THE LOCATION OF THE C.G. OF THE 

13 C AIRCRAFT'ALOMG THE LONG I TUB I N AL ' Axl S . ' THE ME TM 00 ' USED ISHOt' 

U C INTENDED TO ACCURATELY LOCATE THE C.G.» BUT RATHER TO AID 

“15 ■ - c ■■ IN LOCATING THE C.G. SHIFT DUE TO INCREMENTALLY LENGTHENING THE 

16 C AIRCRAFT TAIL CONE. INPUTS REOUIREO FOR THIS SUBROUTINE 

17 C ARE LISTED BELOUi, AND MUST BE ENTERED INTO THE SUBROUTINE 

18 C THROUGH EITHER READ STATEMENTS OR COMMON. 

19 C LN NOSE LENGTH (FTJ" ■" 

20 C LU CABIN LENGTH (FT) 

• 21 ‘ C LC TAIL CONE LENGTH (FT) 

22 C WFUS FUSELAGE WEIGHT (LB) 

... — gj ^ empennage WEIGHT' (LB) 

24 C PHINU NOSE CONE SHAPE PARAMETERS 

25 C ~“PHIN2 

26 C PHICl» TAIL CONE SHAPE PARAMETERS 

■ 27 C PHIC2 

28 C XCGH HORIZ. TAIL C.G. LOCATION (FT) 

. ^ VERT. TAIL C.G. LOCATION" (FT) 

30 C I TAIL CONE INCREMENT NUMBER ^ 

32 REAL M*NfLN/LU#LC fMN^NN^MC /NC#LCCLD 

■ 33 '■ X(A/H,N)- =(.5-(1.0/(M* (N+2.0)))+(1.0-H)/(4.0*(n**2)*(2,0* 

34 •ta^2.0))-((1.0-3.0*H+2.0*H**2)/ (6.0* (M* *3)* ( 3.0*N+2 .0) ) ) ) *A / Cl .0-(1 

— 35 <..0/CM*(Nt1 .0 ))) + ((! .0-M)/(4 .0*(M»*2)*(2 .0*N + 1.0)}) - (1.0-3.0*H + 2.0* 

36 *(H**2))/(6.0*(M**3)*(3.0*N*1.0))) 

38 C COMPUTE FUSELAGE COMPONENT C.G. LOCATIONS 

- j, .. - (. 

40 LC = LCOLO ♦ DELTLC 

- 41 CALL C0NSHP(PHIN1 /PKIN2#T0L*MN/NN) 

42 I F (I .GT.O>GOTO SO 

yQP _ U5R0SS (WC+WV+WH) - 

44 50 CALL CO NSH P (PHI Cl *PHIC 2 # TOL /HC / NO 

45 XC = X(LC/HC»NC) + LN ♦ LU 

46 ARM1 = WGR * XCG 'F WC * XC 

47 C 

48 C COMPUTE EMPENNAGE COMPONENT C.G. LOCATIONS 

50 XV = XCGV + DELTLC 

51 XH = XCGH + DELTLC 

52 ARH2 = UH*XH ♦ WV*XV 




09.9Q5 


0 I t 
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si C 

S'; C eSTUIATE AIRCBAfl C.G. LOCATION 

55 C 

56 UEST s U6R 'T UC MH * WV 

57 XCGEST = (ARMI * AHH2}/NEST 

58 _ _ _ IfCI-GT.O) GO TO 100 

59 C 

60 C COMPUTE CORRECTION FACTOR — CGERR 

61 C 

62 CGERR » (XCG - XC6EST) 

63 XC61 s XCG 

64 GO TO 200 

65 10b XCGr« XCGEST + CGERR 

66 200 XBCG =(XCG1 -XCBAR)/CBAR 

67 RETURN 

68 END 






SUQROUM^rE STABAf?EA(XUARCG/HSUPLE4rSKrBH^CH) 

COMMON/S rAQ/XQARU#XBARAC /D XACB^XBARH 

COM'IOM TAS ^AMACH#VlSCOS^StJ/UW/A W^CBAR^USMPLE^TAPRW^XCBAR^LOO 
lDFUS^LN^LU#LCOLD/FUSWlO/FUSHGT^Ht^H2,PHIWl,PHIW^^PHICT#PNlC2 
2K1 #K2 ^K3/KA^SHO#aHO/AiUSWPHO^TAPRH^ THI C KH ^ LH OLD ^ HH ^ S VO^ 8 VO A 
_ SSWPVO^T APftV^THICKV^LVOLD^AEfF^K#KEMPWKEMP2^S0SN#SaSU^SBSCl ^ 
AULTLOAD/XCG^XCGH/XCGV/ XBARWO^XBARHO^CMA^CNB/KAPAW^KAPAH/KAPA 
5ETA,OELTLC#EP 

COMMDN/CLOATA/ASPECTi^HATA^ S WPC5/KAPPA 

REAL L»0LD^LHNEW#HH^KU8/KAPPA,KAPAH#KAPAUI#K A APR# KH 
SH=SHO 

BH-BHO 

HSWPL E=SWPHd 

THI S PART Of THE PROGRAM OPTOM12ES THE SWEEP ANGLE OF THE 
HORIZONTAL TAIL# FIRST WE CONVERT THE INPUT VALUE OF THE LE 
TO THE HALF-CHORD SWEEP ANGLE 

■ HSWPC 2 = ATAM( { S I N ( H S WPL £ > /C OS (H 5 WPLE ) >* C 4 . / AH * ( . 5 * < I - -T A P RH ) / 
*<1 .fTAPRH)))) 

FIRST WE CONVERT THE INPUT VALUES OF SURFACE SWEEP ANGLES 
TO useable FORM FOR CALCULATION PURPOSES 

WSWPC2 = ATAN( CS INCWSWPL E) /COS(WSWPLE))- C4*/AW* C , 5* C1\-TAPRW) 
*/{1-+TAPRW)>) ) 

WSWPCA=ATAW(<SlN<WSWPLE)/COS£WSWPLE))-(4./AW»C.25* (I.-TAPRW) 
•/n.i-TAPRW J))) 

BETA^SQftni .-AMACH**2.) 

1 CLLH1=0, _ _ 

HSWP1 =0, ■■■ ■ ■' 

.DO 10 I = l#181 
HSWP2SH SWP1 + .00872 7 
ASPECT=AH 
0ATA=BETA 
SWPC2 -HSWP2 

““ KAPPAcKr.PAH ^ 

CALL LSLOPE (SLOPE) 

CLAH=SLOPE 

LHNEW = LHOLD+DELTLC + < ( ( BH/6- )* t1 . +2* *TAPftH) / ( K +TAPRH) ) * 

*( (SIN (HSWP2) /cose HSWP2 ))-( S I N ( H S WP C2) / COS C HSWPC2 > ) ) ) 

CLLH2=CL AH*LHNEW 

" IF (CLLH1.G£,CLLH2 ") G0~ TO 20 

ClLHl =CLLH2 
10 HSWPli=HSWP2 
20 CONTINUE 

NOW WE COMPUTE THE PARAMETER:; TO TH ^ E QUA T I ON FOR SH 

FlRSt WE COMPUTE clam ~ 


SWPC2=HSWP1 

CALL LSLOPE (SLOPE > 


m 



n 


f> 




39967 01 10-1 8-76 09.905 


LABEL STAUAR PAGE 


Cd 

Ui 


5J XLAH^SLOPe “■ '■ " 

5A t 

55 C NOU UE COMPUTE XBARACH 

56 C 

57 xeiARH=X0ARHO4^DELTLC/CBAR + C((BH/6.)*n, +7. *TAPRH) / < 1 . + r APRH) ) * 

58 *(CSIN<HSWP2J/C0S (HSWP2 JJ - C SI M ( HS WPC2 ) / CO S ( HS WPC 2 ) ) ) ) / C Bfl R 

59 ■■ - c 

60 C NOW WE CALCULATE D EPSILON / D ALPHA 

■ 61 C 

62 ASPECT =AW 

63 SV)PC2=WSWPC2 

6A 0ATA=8ETA 

KAPPA=KAPAW ~ 

66 call LSLOPE (SLOPE) 

67 CLAWH=SLOPE 

68 0ATA=1. 

69 CALL LSLOPE (SLOPE) 

70 CLAWO=SLOPE 

KA-CI ./AUj-d./d.tAU**!.?!) ' 

72 KTAPn=C1Q,-(3.*TAPHW))/7. 

73 KH^n CHH/QW))/( (2.*LHNEW)/0W)**,3333 

7^ OEPDAL=<CLAWn/CLAWO>*( A.^iA^tKA*KTAPR*KH*SQRT(C0 5CWSWPCA>)><f*1.19) 

' 75 C 

74 C NEXT COMPUTE CL ALPHA W0 

78 KWa 1 ,-(.25*(OFUS/BW)*»2) + .O25*DfU$/0W 

79 CL AWB=KW0*CLAWH 

80 C 

. .. „ COMPUTE THE TOTAL AIRPLANE CL‘ ALPHA 

82 C 

- ’ ■ gj CL AA=CL AWQ+CLAH*ETA*SH/SW* (1.-DEPDAL) 

8^; c 

85 C . COMPUTE X0ARAC 

86 C 

87 - X0ARAC=XOARC6-CMA/CLAA 

88 call MULTOP(LHNEW/CLAWM^DEPOAL/DXACO) 

X0ARWB =X9ARW '^OXACU ‘ " 

90 C 

91 C ■■ ■ NOW THAT ALL PARAMETERS ARE KNOWN UE CAN COMPUTE THE NEW AREA 

92 C 

93 SHNEW“( XBARUB-XOARAC)/ ( tX0ARAC-XBARH)* C CLAH/CLAWG) * 

9A iKeTA/SW)*n.-DEPDALn 

95 0HNEW=SQRT C SHNEW* AH) ' 

96 BHDI F =ABS (BHNEU-BH) 

97 IF (QHD 1 F.LC. EP) GO TO 3G 

98 SJr=SHNEW 

99 OH=aHNEW 

1UD GO TO 1 

■ ~'101 • “■ 30 CONTINUE ' 

1 02 SH=SHN£U 

103 0H=3HNEW 

10^ CH=(A.*SH)/C3.*BH)-<n.+TAPRH*TAPRH**2.>/C1.+TAPRH)**2.) 


V. 
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LABEL VERTAfi .PAGE 


e 


[S?_ 


1 

2 

3 

5 

6 

”7 

B 

9 

10 

11 

12 

1 4 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 
‘ 25 

26 

27 

28 

29 

30 
^31" 

52 

33 

34 
■35" 

36 

37' 

38 

39 

40 

41 

42 
4 3 

44 

45 

46 

47 

48 

49 

50 

51 

52 


SUOnOUUNE VEHTAREACVSWPLE / SV/OV/CV) 

COMMON TAS/AMACH^ VISCOS/SW^BW^AW^CBAR, WSWPLEi'TAPRN^XCDAR,LOOLD# 
10FUS^LN/LU^LCOLO^FU5WIO/FUSHGT^Hl#H2#PHlWl#PHlN2/PHiCl#PHIC2^TOL# 
2K1 ^K2 ^K3i^K4^SHO#BHO/AH#SWPHO#TAPRH^THICKH/LHOLD^HHr SVO^aVO^AV# 
3SWPV0,TAPftV^THlCKU^LV0LD^AEFf/^K^ KEnPUK£MP2^SBSNySBSU^ SBSCWDIVE^ 
4ULTL0A0/XCGirXCGIt/XCGV# XQ AR V'B^ X 13A RHO ^CM A# CNB / K APAU# K AP A H ^ KAPAV^ 
SETAi-OELTLC^EP 

cc::mon/cee6EE/xcgi 

REAL K/LVOLD#LVWEW^KN^KRL^KA^LCOLO,LCNEU^LBOLO/L0NEW#LB2SSS 
real hi i-H2vHlH2/M1 / HWRAT/M2^KAPPA/KAPAV#KVr KC/KW/LN 
C OMMON/CLD ATA /ASPECT ^13 AT A/ SWPC2/ KAPPA 
trV^BVO 

VSWPLE'=SWPVO 

THE FIRST PART OF THIS ROUTINE 0PT0M12ES THE SWEEP ANGLE Of 
THE VERTICAL TAIL FIRST WE CONVERT THE INPUT VALUE OF LE SWEEP 
TO HALF CHORD SWEEP 

VSWPC2=ATAN( CS INCVSWPLE) /C05<VSWPLE))-(4./AV* ( , 5*n .-TAPRV) / 

*(1 .+TAPRV) ))) 

B£TA=SQRTC 1 .-AMACH**2. ) 

1 CLLVl =0. 

VSMP1 =0. 

DO 10 I=m8l 

VSWP2=VSWP1+.00B727 ’ ' ' ' 

_ ASPECT=AEFF 
BATA:^BETA 

SWPC2=VSWP2 _ 

KAPPA = KAPAV ^ 

CALL LSLOPE (SLOPE) 

CLAV=SL0PE ■ " ^ ' 

XVOLD = LVOLD - CXCGl-XCG) 

* LVNEW=XV0L0+DELTLC+<((BV/3.)*(2.+2.*TAPRV) /{l.+TAPRV ))* 

_*C(SIN CVSWP2)/C0SC VSWP2 ) )-( SINC VSWPC2)/C0S(VSWPC2>) )) 

CLLV2 = CL AV^LVNEW • . . 

IF (CLLV1.GE.CLLV2) GO TO 20 

CLLVl “CLLV2 

10 VSWP1=VSWP2 
20 CONTI NUE 


THIS PART Of THE ROUTINE COMPUTES THE NEW VERTICAL TAIL AREA 
AND VERTICAL TAIL SPAN 

>IRST' WE CALCULATE KRL ■' 

lbnew^lgolo+deltlc 

RENUM F=RENUM(OENSIT^VI SCOSi^ TAS/LBNEW) 

KRL= -1 ,830754+. 20494*AL0GCRENUMF) 

^NbW WE calculate THE DOO Y S I DE ' A R £A OF 'THE NEW AIRPLANE*" 

lcnew=lcold+deltlc 

CALL CONPAR ( P H I C 1 /K A , K V /K C /K W ) 


© 


9 


1 


3 
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LABEL VERTAR PAGE 


SU5ECl-KA*LC0Lb*DFUS 

SBSRAT = 1,00 

SQSEC2=KA*LCNEW*DFUS 

CALL C0NPAR(PHIN1 /K A^JC / KW ) 

SBSN = KA»LN*0FUS 
Sl3S=SBSNfSBSU + SBSHAT*Sd$EC2 

MOW WE CALCULATE KN BY A THREE STEP 

METHOD IT WAS ARRIVED AT BY APPLYING CURVFITTING TECHNIQUES 
TO the chase around chart in DATCOM 


LQ2SUS = LBNEW**2,/SBS 

'i F' (L02SBS.GE.8.) GO T 0 5Q 

SH IF 11 = 6,094 2-1 .3 516*LB2S0 S+,1 5525i»l02SBS** 2,-,O062*LO2SBS **3,+ 

■ ■ *,0001 *LB2SBS 
GO TO 50 

30 IF CLQ2SBS,GG,12,) GO TO 40 
SHIFT 1 = -, 12* LB 2 SB 5*1, 91 

GO TO 50 ^ 

40 SHIFT1=-,05875*LB2SUS+1,175 
50 CONTINUE 

YVALUE=2,8353*XCG/LBNEW-,4 1667+SHIFTl 
fn H2 = SQRTCH1 /H2) 

HI =3.64 9 7-3,5796*h1HH-,3 9*H1H2 **2, + 2,31 49*H1H2**3.-.6946*H1H2**4, 

2VALU E=YVALUE7M1 ' ^ 

HWHAT=FUSHGT/FUSWID 

M2 = C-1.014 7 + 4,4649*HWR AT-3,362 6*HWRAT**2. + 1 .0794*HWRAT**3.-.1217- 
*KWRAT**4,)*.001 

KN=W2*2VALUE-,0005 “ ' ^ 

SWPC2=VSWP1 

CALL LSLOPE (SLOPE) 

CLAV=SLOPE 


COMPUTATION OF THE NEW PLANFORM PARAMETERS 

SVNEM=;(SW*BW* CWB + 57,3* KN *K R L *S B S * LB NEW) /(K*CLAV*LVNEW) 

CALCULATION OF THE NEW SPAN AND COMPARISON WITH THE TOLERANCE 
OF THE ITERATION 


OVNEW=SORT CAV*SVNEW) 

BVDI F =ABSCBVNEW-BV> 

IF CDVDIF.LE.EPI 'GO'TO no ' ' 

SV=SVNEW 
OV=BVNEW 
GO TO 1 
1 10 CONTI HUE 
SV-SVNEW 

BV^OVNEUI - . - . 

CV = (4,*SV)/C3.*BV)*C(1 ,+TAPRV+TAPRV**2, )/C1 • -FT APRV ) * *2 . ) 
VSWPLE = ATAN( CSINCVSWPl ) /COS fVSUPl ) )“(4 , /AV* <-,5* Cl ,-TAPRV)/ 
*C1 , + TAPRV)))) 


10 
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LABEL HULTOP .PAGE 




W 

w — 
o 


1 ■ SL/BftOijTlNE 'r'lULTOP(LKnEW/CLAWM^bEPDAL/DXACB) 

2 COMr^ON/RHO/DENSlT 

3 COWMOW TAS/ AHACH/.V1SCOS#SW#EJ,J/AW^CBAK^WSWPLE/TAPRW/XCBAR/LBOLO^ 

A 1 DFUS/LN/LU^LC0LD/FUSWI0^FUSH6T/H1>H2/PHINWPH1N2/PHK1 / PHIC2/T0L/ 

5 2K1 /K2^K3/K^^SHO/eHOi.AM^SWPHOj-TAPRH,TBlCKlULHOLD/HH^SVO/nVO^AVy 

6 _ 3SUlPVO,TAPRV,THICKV,LVOLD^AEFF,Ki‘KEMPUKEMP2/S0SN^SBSU^SBSCWOlVE^ 

? ■■"AULTLOAD^XCG^XCGH^XCGV^XBARWe^XBARHO^CMA^CNB^KAPAW^KAPAH/KAPAV^ 

8 SETA^OEL LC^EP 

9 COHnON/-HAPE/MM^NN^MC>^NC 

10 DATA CQ^Cl ^C2/C3^C/i / I . 90^ - 1 . 6 95 8^ 1 . 57 5 9^ -0 . 729 2^ 0 . 1 30 2 / 

11 DATA OOaDI ^02/ 03.- da /6. 2 503 #-23*0908/60. 35 53 #-78.45 A0^38. 2895/ 

12 REAL L«/NN#MN#LT#LBOLO#LCOLO#LCNEW#NC#WC#LH #LHNEW _ 

iy POLYfX/C#Cl.pC2VC3/C4) = C + Cl *X + Cy*X**2.*C5*X**3p+C4*X'**4. 

U CLAWM := CL A WM/ 5 ?• 29 5 8 

15 QBAR » 0.5*DENS IT*TAS*TAS 

16 XLE=XC0AR-<<(BH/6*>*(1. + 2.* TAPRH) / Cl . f T APRlO )- . 5 *0 F US > IM CWSWPLE) 
M 412 BH IS MOT DEFINED 

17 _ */C0S (USWPLE) 

l6 ■ CR = y . 5 *CBAR* f1 . +TAPRU ) /( 1 . VT APRU + TAPHW-^*2 J 

19 CHF = DFUS *CR* (TAPRW-T . ) /BW +C ft 

20 DX I XLE/50. 

21 SUM = 0. 

22 ■ ■ " !?0 100 1 = 1 #50 

23 XI - XLE - <DXi*Q.5>- 0XX*CI-1) 

24' I F CXI ,LE , CKLE-LN) ) GO TO 10 " ' " ' ' " ' 

_ 25 _ WFI = OFUS * C (1 .-( ( (LNf XI-XL E) /LM) **NN) ) ** (1 , /HN)) 

'26 ■ ■ GO T0'20 

27 10 WFl = DFUS 

28 20 XCf = Xl/CRF 

29 DEOA = POLV(XCF#CO/C1/C2^C3#C4> 

IF(1,GE.45) OEDA = POLYtXC F#D0#D1#02#D3#D4 J 

31 DEDAI = DEDA*CLAWH/A.5S 

32 ' 100 SUM = SUM + ( 0 X I * W F I * W F I *D E DAI) 

33 LT = LBOLO + DELTLC - XLE - CRF 

"34 " "OXI^LT/SO* 

35 LCMEW -LCOLO ^DELTLC 

(JQ 1=1 #50 ■ ■ ■ ■ 

37 XI = DX[*Cl-l>+(0.5*DXI) 

"38 1 f <X I .LE, CLT-LCNEW) ) GO TO 30 

39 WFI = DFUS * ( n •-< ( CLCMEU+X 1-LT ) /LCNEW) **MCn A* CT./MC) ) 

40 ■ ““ GO TO 4 0 

41 30 WFI = DFUS 

42 40 LH = LHMEW-CRF +. 25 *CBAR + C CBW/6 , > * Cl ,+2 •* TAPRM) / (1 . + TAPRW) -DFUS /2. ) 

43 *^SIN(WSWPL£)/C OS C WSWPL E ) 

44 DEDAI = (1 ,-D£PDAL)*XI/LH 

45 200 SUM =i SUM + C DX I *WF I *W F I *D E DAI ) 

46 DnOA=CQBAR/36-5)*SUn 

47 DXACB = -1 .*DMDA/(QBAR*CBAR*SW*CLAWH> 

™ return 

49 END 




i..= : 


IZ 
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LABEL CONSHP f^AGE 


u> 


209 


^ 

2 

3 

A 

5 

6 

7 

8 

9 

10 

11 

12 

13 

U 

15 

16 

17 

18 

19 

20 
21 
22 
23 

2A _ 

_ 26 

■ 27 

28 
_ 29 

30 

31 

32 

- 33 ■“ 

3A 

35 ■ ^ 

36 

37 

38 

39 

STATEMEWT 

40 


10 


20 


30 


40 


50 


.60 
■ 70 
80 


SUBROUTINE C ONSHP C Pil 1 1 / P H! 2 / tOL ; 
REAL n#W#MWM2 
1 = 1 
J = 1 
W = 1,Q 

DM2 = 10.0 

Ml = ALOG<t.0 “(PHn*»Nl)/ALOG(PHI 
m 2 = ALOGd.O -CPHI2**N>)/AL06(PHI 
J = J+1 

IF(J.EQ.IOQ) GO TO 80 
DM1 - Hi -M2 

IFCOHl ,LE.O.) GO TO 30 
F“i DM i : G t . D M2 1 60 TO 5 0 
tJ = N + (0.05/l) 

DM2 = OHl 
GO TO ID 

IFCAQS(D«1) .LT.T0U60 TO 60 
I = 1*5 

N = N - (0,057l) 

Ml = ALOGd.O -(PHI1 **N))/AL06(PH1 
«2 = ALOGCI.Q -(PHI2**N))/AL0G<PHI 
J = J+1 

IF CJ.EQ.100) GO TO 80 

DH1 = M2-H1 

"IF COMl.LT.O. )“ GO' TO 20 ' 

If CABSCDMD.LT.TOL) GO TO 60 

GO TO 40 

WRITE(6#n 

M « 1 , 0 

U = 1.0 
■G'O“t0^7 0 
n = CMl + M2)*0,5 
RETURN 
WRITE (6^4) 

GO TO 60 


N) 


n 

2 / 2 , 0 ) 


1 > 

2 / 2 . 0 ) 


1 F0ttMAT{10x,****!TERATlON FOR H AND N DIVER 

x*722x^‘w=1 .0*//) ■■ ■ ‘ * ■ 

4 fORMATC 1QX#'***100 STEPS COMPLETE— DID NOT 
.RETURN * ' 

cannot be reached 

END 


6ES***V1 

CONVERGE 


5X / < SET- M=1^q 
A***//) 


=. -7 


<5 


n; 

U 


a ■ . ,. J 


fr 


ii ■ : - J 


. y 


1L ,^,.,.^1 if..i,:.,^,.:i 

t \ J U : - ■ T 




rl- 


^ ■'5 
t/ 


&r-^3 0^] 
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LAOEL CONPAR .PACE 1 


2 

3 

5 

6 

8 

9 

10 

11 

12 

“B- 

u 

15 

16 

■ 1 ? ■ 
18 
--19 


SUUROUT tNE CONPARCPHI/ KA^kV^KC^KWj ' “ ' 

REAL KA/KV#KC^K« 

C THIS SUBROUTINE COMPUTES THE AREA CORRECTION FACTORS 

C KA^ KV# KC/ AND KW WHEN THE SHAPE PARAMETER PHI IS INPUTTED 

C THIS SUBROUTINE WAS DERIVED FROM TORENOEEK PC. A47 

C 

DATA AK^AKl , A K2 / AK5# AK A /-O * 59^ 3 , 8 T 0 9#-5 . 721 #6*A1 68/ -2*9167/ 

DATA VK/VKI / VK 2/ V K3/V k 4 /-. 9 095 / 5 . 8 0 3/- 1 3*092 7 / 1 6. 5 927/ - 7 * A 074/ 
DATA CK/CKl/CK2/CK3#CK4/2.96/-1 2.0488/ 22. 87 52/- 1 8. 3335/5 .5 556/ 
data WK/WKI /WK2/UK3/UJ<4 /-1 72.7 21 /1008. 8388/ -21 6 2 , 2 5 / 202 3 .9 B77/-69 
*7. 9097/ 

P0LY(X/C/Cl/C2/C3/C4)sC + Cl *X + C 2 * X * *2.+ C 3* X * *3 . + C4 * X **4 , 

KAcpdLYCPHi;AK;Ak1/AK2/AI<3/AK4l " 

KV=POLY CPH 1/VK/VK1#UK2#VK3/VK4 ) 

KC=P0tY(PHI^CK/CKl/CK2/CK3/CK4) 

KW=POLY<PHI^WK/WK 1 /WK2 /UIK3 /WK4 > 

CONTINUE 

RETUR^N _ __ 

""‘EN D ■" • - ■ ^ - - - ^ 


w 


U) 

(O’ 


I 

I 



B.34 


i t I t ; - t t,.,- i i 

LABEL RENUt^ PAGE 1 

FUHCtlbW RENUn^DENSlT/VIS‘cbsVTAy#AL»G^ 

THIS FUWCTIOH COMPUTES THE REYHOLOS NUMBER OF A QOOV 

RENUn=(TAS*ALNGTH)/VlSCOS 

RETURN 

END 





C 




o 




B3 Exanple Output 


This section provides Rosicam/Filltnan Program output for a 
conceptual aircraft configuration. This conceptual aircraft is 
not to be confused with any of those mentioned in Chapter 4. 
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B.36 




< I 




I 


g. 


TABLE 2.15.TA ROSKAM- F I LL MAN METHOD 

CONFIGURATION J (31 PAX/ 2/2) 

'original* flight conditions and' 

GEOMETRIC definitions 


flight conditions 


ALTITUDE 25000.000 FT 

AIRSPEED ” ■ ‘ ‘ ■ 366, 750 F PS ~ 

MACH Q-561 

DEMSITY 0.00107 S L UGS /CU , FT , ' 

KINEMATIC VISCOSITY 0.00030 SQ.FT, /SEC 


GEOMETRIC _ 

UING 

REFERENCE AREA 440,000 SQ.FT. 

SPAN 66.330 FT 

ASPECT-RATrO TO- 000' 

M.A.C. 6,8ti0 FT 

f ,S. ^ " ■ 19,840 FT ’ 

L.E, SWEEP 0, DEG 

TAPER " ’ 0. 500 


TUSEtAGE — 

OVERALL LENGTH 56,350 FT 

— NOSE LENGTH U;620 FT 

CABIN LENGTH 20. QQQ FT 

TAIL CONE LENGTH ""■21.760 FT 

MAX WIDTH 8.080 FT 

— MAX HE I GHT 8: 080 -FT 

EOUIV. DIAMETER 8,080 FT 


HORIZONTAL TAIL 

REFERENCE AREA 
SPAN 

ASPECr‘RATlO 

M.A.C. 

TAPER ■ 

THI CKNESS 


55.600SQ.FT, 
14.910 FT 

-4,000 

3. 870 FT 

■ 0. 500 ■■■" ” ■ 

0.090 


VERTICAL TAIL 

REFERENCE^-AREA- 

span 

ASPECT RATIO 
M.A.C. 

TAPER 

THICKNESS 


4a.380SQ-FTV 
6,960 FT 
1-000 
7.210 FT 
0. 500 
0,0 90 


TAOLE 2.15,1a ROSKAM-f ILLM.^W METHOD 
_ COMF I GURATION J (31 PAXj- 2/2) 

original weight awd balance statement 


(LBS. ) 

GROSS WEIGHT 

22000.000 

APPROXIMATED SHELL WEIGHTS 


FUSELAGE 


N(DSE SECTION 
CABIN SECTION 
TAIL CONE SECTION 

75^.560 

12A6.399 

126A.082 

TOTAL FUSELAGE 

3265. OAO 

Empennage 


HORIZONTAL tAIL 
VERTICAL TAIL 

115.226 

_ 3 30. A 31 

TOTAL EMPENNAGE 

4A5.657 


CENTER Of GfiAVITy 


FRACTION Of M.A, C 
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TABLE 




2-15-IC 


ROSKAM-f ILLHAM rtETHOO 
COWFIGURATION J <31 PAX^ 2/2) 
ORIGTNAr WETTEb 'Al^EA 'BREAKDOWN 






FUSELAGE (SQ.FT.) 

NOSE SECTION 3C7.346' 

CAOIN SECriON 507.601 

tail section 514. 803 ' 

TOTAL FUSEL AGE ■ ' 1 329.910 " ' 

EMPENNAGE 

HORIZONTAL TAIL 113.079 

VERTICAL TAIL 98.395 

TOTAL EMPENNAGE 2 1 1 . 4 75 








1 

15.000 

49.620 

2 

1 5.250 

_A9,B70 

3 

15.500 

50 .1 20 _ 

ti 

1 5.750 

50.370 

5 

16.000 

50.620 

6 

16.250 

50,870 

7 

1 6.500 

51 .1 20 

0 

1 6.750 

51 .370 

9 

17.000 

51 .620 

10 

1 7.250 

51 .870 

-- cxJ — 

. 11 

17-500 

52 .120 

\ 12 

1 7.750 

52.370 

1^ 

1 B.QOO 

52.620 _ 

1 4 

10.250 

52.070 

1 5 

18.500 

53 .120 

16 

1 8.750 

53.370 

1 7 

1 9.000 

53.620 

10 

19.250 

53-870 

1 9 

19.500 

5^. .120 

20 

1 9.750 

5^.370 

21 

20.000 

5^i.620 

22 

20.250 

5^.870 

23 

20.500 

55 .120 



20.750 


55.370 


N METHOD 


i J (31 PAX^ 

2/2) 


WEIGHT and balance 

. ) 

20290.261 

(LBS. > 

fuselage 

EMPENNAGE 

C .G. 

WT 

WT 


(LB) 

(LB) 

(F0 CBAR) 

21066,239 

439.115 

0.148 

21062.91 2 

793.682 , 

0.205 

21080.119 

832.91 5 

0.216 

21097,460 

797.579 

0.214 

21114.933 

791,ufi2_ 

0.217 

21132.538 

^ 788.696 

0.221 

211 50.275 

786.583 

0.224 

21163. 142 

784.563 

0.228 __ 

21186,141 

782.504 

0.232 

21 204,269 

782.219 

0-236 

21222.527 

780.21 3 

0-240 _ 

21240.91 5 

777.981 

0.243 

21 25 9, 432 

__ 7 7 5. 6 2^6 

0.247 

21278,077 

773.334 

0,251 

21296.850 

771,047 

0.255 

21 31 5,751 

768.685 

_ 0.259 

21334,780 

766.367 

0.263 

21353.936 

764.168 

0,267 

21 373.218 

761.934 

0.271 

21392.627 

759.563 

0.275 

21412.161 

7 5 7,1 95 

0-279 

21431.821 

754.873 

0.283 

21451 .606 

752.560 

0.2B7 
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0,295 

26 

21.250 

55 ,870 

2151 1 ,709 

745.562 

0,299 

27 

21.500 

56,120 

21531 .991 

74 3-2 49 

0.303 

28 

21 .750 

56,370 

21552.597 

740.889 

0 . 307 

29 

22.000 

56,620 

21 572.926 

730 . S 66 

0,312 

30 

22.250 

56.870 

21593.578 

736.295 

0.316 

3 T 

22.500 
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21 656.266 

729.313 

0.329 

3 ^ 
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0.333 

35 
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21698,667 

724.656 

0-338 
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23.750 
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21720.049 

722,344 

0.342 

37 

24.000 
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2 1 741 . 5 51 

720,008 

0.347 

38 

24,250 

58,870 

21 763.1 74 

71 7.700 

0.351 

39 

24.500 

59.120 

21 784.916 

71 5,396 

0.356 

4 J 

24.750 

59.370 

21006.779 

71 3.066 

0.360 

4 1 

25,000 

59.620 

21 82 8.761 

710,760 

0.365 

42 

2 5.250 

59,870 

21850.862 

708.514 

0,369 

43 

25.500 

60,120 

21873.081 

706.203 

0,374 

44 

25.750 
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21895.420 

703,907 
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45 

26.000 
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2191 7.876 
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61 ,1 20 

21763.1 43 
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21985.953 

694,747 
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27.000 
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22008.880 

692,450 
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22031 .925 
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687,907 
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62 .870 

22125.265 

681 .076 
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2 a - 5 00 

28.750 
29.000 
29.250 
29-500 

29.750 


63.370 

63.620 

63.870 
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64-370 


221 ^ ^ a . 815 
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22244.542 669-856 
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2.15.20 fiOSKAM-F ILLHAN METHOD 

CONFIGURATION J <31 PAX^ 2/2) 

— 0PtlMl'UTidN'“StABtLlZER' 

SIZING 

CMA = -1.000 

ASPECT RATIO == A. 000 

TAPER RATIO = _ _ ^-500 

I LC SH BH CBARH SUEEPH 

CFT) CSQ*FT.) CFT) (FT) (DEG) 

1 1 5701)0 3r.33T 71 772 0 STTB S ^T‘574 7 2‘ 

2 1 5.250 4 9;3 3B 14.048 ' 3.6 4 2 19.854 

3 15.500 51 .61 5 1 4.369 3.725 19.854 - - 

4 1 5 ;7 5 a 50 ;T5^2 2 3 0 3^ 6 89 T9 ^ B 34 

5 1 6.000 ’ 50.847 14.261 “ 5.697 19-854 

6 ‘ 16.250““ 51 .238 14.316 3.712 19. 854 

7 16 :500 5T-66 7 1 4‘376 3 7727 19:3 94 

8 1 6.75 0 52.081 “1 4.4 33 3,742 19.3 94 

- 9 - "17,000 5274 92 ““ 14,4 90 3.7 57 19.394 

10 17.250 52:978 J-4';557 3.774 19.394 
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■ • 19.394 ■ 

1 2 ■ 
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’ '■ 19,3 94 ^ ■’ ■■ 
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19; 39 V ' 
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- 19.394 
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- -18.933' - - 

20 

1 9,750 
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3.926 
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'■ ”18.933 ■ ■ ' ■^ “■ 
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2,1S,2C ROSKAM-F ILLHAi^ METHOD 

COMFIGURATIOW J (31 PAX/ 2/2) 
OPTlHIiATlON VERTICAL^TAIL 
SIZING 

CNB ~ 0.120 

ASPECT RATIO = 1.000 

TAPER RATIO = 0-500 


******** ***** *********** ************** **** *********************** 


I LC SV 13V CBARV SWEEPV 

(FT) * “(SQ.FT.) (FT) (FT) ' (OEG) 
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15.500 ■ 

TT7".11 9 

" 1 17.577 

■ 110, A2a - ‘ 
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■ 1 1.245 
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■^69Vi:Br 

6 9.1 83 ‘ ■ 
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— /p 
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- 6 
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10^7786^ 

' 109.238 

10.478 
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10.866 

10.839 ‘ 

68.2/2 

■67,971 ■ ■ 
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10.426 

10.812“ 

67.673 - - - 
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6". 5 00 — 

TOB^I 5"6 

1 0.'4 00^ 

1 0,TB5“ 

67-376 
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1 6.750 

“”107.614 "■ 

"10.374 

10,758 

67-OaO ^ 

’ 9 

■ 17.0Da 

““107 ,380^ 

‘■■““10.36 2 

10.746' 

66.786 ■ 

^10’ 

rr;-?so ” 

TO 6781 

1 0.335 

10;718 

66 • y*t 

11 

17.500 

106. 2«« 

10.307 

■10.t89 

66.203 

12 

“17.750 

~ 105.658 “ ■ 

10.279 

10,660 “ 

66.203 ■ 

■13 

-“~i 8;ooo 

1057085 

10.251- 

“ 10.631 

"65.914' - ■ - 

1 A 

18.250 

104.51 1 

10.223 

10.602 

65.626 

15 

18.500 

103 .94 0 

10.195 

- - 10.5 73 

65-340 ■ ‘ - 

“16 

^ 1 S . 7 5 0 

1 03.355 

1 0.1 66' 

1 0.5 43 

65.340 

1 7 

1 9 .000 

102.780 

10.138 

10.514 

65.054 

1 6 

19.250 

102.207 

10.110 

10*484 

64.770 

- 19 

19.500 

‘ — 101 .633' ' — 

”10.081 

' 10.455 - - 

64; 4 88 "■ 

20 

19.750 

101 .050 

10.052 

10.425 

64.488 

21 

20.000 

100.475 

10.024 

10.395 

64.206 

-22 

?0.250~- — 

99. '901 

— - '9.995 

10.365 

63. 926 

23 

2 0.500 

99.51 8 

9.966 

10.335 

63.926 
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25' ■ ' 

21.000 

■ 98.1 70 
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10.275 

63,369 ”■ 
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21.250 

97. 5^96 

9.879 

10.245 

63.092 
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21.500 

V7.D1 5 

9.850 

10.214 
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21 .750 

96.441 

’9,820 

10.184 

62.816 
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22.000 

95.868 

’ 9.791 

10.154 

62.541 

- - ■ 

.. . ... . - 
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22.250 

■■ 95.289 

9.762 

10,123 

62.5^1 




~‘2 2;TOO'" 

94,71 6 

9.T32 

10.093 
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32 

22.750 

94.144 

’ 9.703 

■" 10.062 

61.994 



33 

23.000 

93.567 

'9.673 

10-031 

61.994 ’ 




”3 '2 3 ; 2 5 0 9 2V9 9 5 9 3 1 Q , Q 0 1 6 1,721 

35* 2 3.500 9 2 .4 24 9.614 9.9 70 61.4 50 “ 

36 23.750 91 .B49 '9.584 9.939 61,450 ~ 


3 7 2l\ 000 9 r. ?7 8 9" 5 5 4 979'Q S 51 “I 80 
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““ 9.524 

9.877 
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39 
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60.910 
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4s, 
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4U 

■ '24'.'7SD" 

89,5 66 

9.564 

9.814 

so. 641 


41 

25.000 

” ■ ‘ 88,998 

9,434 

9,783 ‘ 

60.375 


42 

25.250 

88.426 

9,404 

■ 9.752 

60.373 


’ 43 ’ 

■“'2 57500' 

S7.d59 

9,T73~ 

9:720 
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44 
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87.292 ' 

9.343 

' 9.689 

59.838 


45 

26.000 

85.722 
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9,657" 
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"8571Ti" 


^ 59 ". 5 72 -' 


47 26.500 85.592 9.252 9.594 59,307 

48 26.750 ■ 85.024 9.221 9.562 59.307 

49 '27.000 r4.'45'9 ”^9.190 9;S31'~ 59.042 

50 27.250 83.896 9,160 9.499 58.777 

51 27.500 83.330 9,129 9.467 58.777 

52 2 7,750 82V76B" 9.098 '9:435 58: 513 

53 28.000 82.203 9.067 9.402 58.513 


38.250 


81 .642 


9.G36 


9.3?0 


58.250 


54 

55 28.500 “ 81 . 081 ” " 

56 28.750 80.519 

57 29.000 79.959 

— 58 29.250 79V400 

59 29.500 78.839 

60 29.750 78.281 

6 7 3 0700Q 77*;72'1 


9.005 “ 9.338 57,987 

8.973 9.306 57.987 

8.942 9.273 57,724 

5.917 9.247 57.462 

8.879 9.208 57.462 ^ " 

a. 848 9,1 75 ■ 57.200 

8;'8 1‘S' 9TT4 2 57 '.’Z 0 0 
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TABLE 2*15.20 ROS K A M- F I LLHA N HCTHOD 

CONFISURATION J (31 PAX^ 2/2) 
S T A T I C "M A fi 6 1 N C O M f» 0 N £ N T S 

CMA “ -1.G00 


"n,C 


'XBM~ 


XBA 


A 

' 5 
7 

- 8 
— 9 — 
10 
11 

13 

U 


“XST X0UB XBR‘ 

15.000 0.230-0*023 0.207 3.771 0.333 

15.250 0.230 -0.023 0.2Q7 3.821 0.388 

**15:500— 0V230'^D’23^^ 0.207 — r.H58“07398 
1 5,750 ■ 0.230 -0.023 0.207 -3.894 " 0.396 

16.000 0.230-0.023 □,207‘3.930 0.399 ’ 

16 725 0 — 0t 2 3 tT^OrO 25 — 07207 5 .^66 — 074^0 3“ 


— XBCG' 
0.148 
0.205 

'"0721 r* 

0.214 

D.217 

177221 


16.500" 0.230 -0.023 0.207 3 .999 Q .4 06 0 . 224 
16.750 * 0,230 -0.023 0.207 4 .035 0.4 10 -’0.228 
"1 7,000^-aT23a™Q7lT23^’'0.'207 — 4 .071 ’"07 t;i 3 — tT7Z37- 
1 7.250 0.230''-D.024" 0.206 ~4.108 0,417 - 0.236 

17.500 0.230 -0.024"' 0 . 206^ 4 . 1 4 4 " 0 .4 2 1 — D.240 

1T .75 0 — 0 723tl**-=ini324 — 0 ,’20 6 — 471 8 1 — 074 2'5 D72 4"3"' 
18.000 0.230-0.024 0.206 4 .21 7 0 .4 28 ^ 0. 247 

18.250 0.230 -0.024 "0.206 4.253 0.432" 0.251 



- 1 5 — 

“18; 500“ 

"07230"-0:024-’ 

"0;206" 

4 -.290'*"’ 

0.4 36- 

"0V255"" 


16 

1 B.750 

0.230 -0.024 

0.206 

4.326 

0.440 

0.259 


17 

19,000 

Do230‘-0.025 

0.205 

4.35B 

0.443 

0.263 

... — .... 

■-18-- 

19.250 

072SU--0.02S-- 

0.205— 

4.395 - 

■0.447 ’ 

0.267 - 


1 9 

19,500 

0.230 -0.025 

0.205 

4.431 

0,451 

0.271 


20 

19.750 

0.230 -0,025 

0.205 

4,467 

0.455 

0.275 
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■ 21" 

‘20,000” 

-■0.230-^0:0 25 

0.205-^^ 

4.504 

0.459" 

■ 0,279 '■ 

^ ' 

22 
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0.230 -0.025 ■ 
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4.540 

0.463 

0.283 
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0. 230 -0. 025 
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-2ti 

-2 0.-75O-" 
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-O'. 20 4— 

■ 4 . 6 I 3 * 

0.471- 

-0.291"- 


25 

21.000 
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4.649 

0.4 75 

0,295 
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“2 7 
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'0.230~-0.026“ 
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~BT2 3 0~^ro'2 T” 
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4 -972 


3 5“ 
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■fl.2£l8“ 

's:o4 5 ■ 
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0.202" 

5.117 ■ 


■■ 39 

“24.500^ 

‘“0'.230"'-070 28“ 

'O'. 202~ 

5;i5'4“ 

V 

4 0 

24.750 

0.230-0.028 

0.202' ■ 

5-190" 


41 25,000' 0.230 -0,028 0.202 “5,222 


■■ ■ ~ 42 2 5 72‘5T3~ 0.23 0“- 0':D"2 8 ” D . 2 0 2~5". 258^' 

«3 25.500 0.230 -0.028 0,202 5-294 

44 ‘25,750' 0.230 -0.028 0-202 5.331 

‘ 4 5 ?6.0a0'~D. 230'-07D29'~D,20r"5‘.367“ 
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0,452 

-- 6 0 - 

29.750 

■ Q.230'-0.031 ■ 

0.199 

■5,907 

— 0.6 3 5 

-'0,458 

61 

30,000 

0.230 -0,031 

0,199 

5.94<i 

0.6(0 

0.463 


P, r • -;-rj t 






Tg’a 






— li' 


r. 


6 

21.25D 

131 7.842 

120.511 

198.492 

1636,845 

1 

21.500 

1323.758 

121 .294 

197*309 

1642.361 

8 

?1 .750 

1 329.673 

122.070 

196.142 

1647,886 

9 

22.000 

1335.589 

122.885 

194.976 

1653.450 

0 ■ ' 

22,250"- 

■~T34r,5D4 

123.665 

1937800 

1658.969""' 

1 

22.500 

1347.4 20 

124-440 

192.634 

1664.495 

2 

22.750 

1 35 3.3 35 

125,217 

791.471 

1670.023 

5 

■ 23-000" 

“" 1359.251 

125.994- 

1907296" 

1675.540 “ 

4 

23.250 

1365.1 66 

126-767 

189.134 

1661 *068 

5 

23.500 

137T.082 

127.543 

■ 187.972 

1686.598 

6 

23.7S0-^ 

F3767997 — 

; 128,317-^ 

186,80 2- 

169 2;-T1 8 

7 - 

24.000 

1382.913 

1 29,093 

185.642 

1697.647 ■ 

8 

24.250 

1385.828 ' 

129.866 

18 4 .4 83 

1705*178 '■ 

9 — 

24.5 00 ^ 

T39 4;r4 4 — 

1 307639 

18 3.3r6T 

1 70 8:699 

0 ‘ 

24.750" 

1400.659 

131.411 

182,159 

— 1714.228 " 

1 

-" 25.000 " 

• — 1406.575 “ 

132.220 

IBl .004 

1719.799 

2 “ 

25.250" 

141 2.4 90 ^ 

1 32*998 

17 97'84r 

1725:329 

3 

25.500 

■ ■ 141 8.406 

133.770 

178.687 

1730.862 

4 

25.750 

1424. 321 

134.543 

177.555 

1736.399 

5 

26 . 000 

143072 37 

■ 155.317 

17 6,376 

1741,930 

6 

26.250 

1436.1 52 

136,088 

175.225 

- 1747,465 

? 

26.500 

- 1442.068 

136-862 

174,076 

1753.006 

8 ' 

26* 750 ” 

144 7. 983 

1 37.637 

172i921 

1 758.541 

9 

27.000 

1453.899 

138.410 

171.774 

1764.082 

0 

27.250 

1459,814 

159.184 

170.629 

1769.627 

1 — 

27.500 

1465.729 - 

139.960 

169.477 

1 775.1 66 ‘ 

2 

27*750 

1471.645 

1 40,732 

168-333 

1780.710 

3 

28,000 

1 477, 5 60 

141.506 

167.185 

1786,252 

; 

28^.2 50“ 

^ 4 8 3 * 4" 7 6 

I 42.280 

” 166.043 

1791 ,799"" 

5 

28,500 

1 489.391 

143,055 

164 .903 

1797,350 





4 — , 

1 t ^ 

, . t 
~t ■ t 

a :. ~ i 

e j ^ 

1 » j 

JS -1 • ■* 

f- • 4 . ■ J 



36 

78,750 

U95.307 

U3.S66 

163.759 

1802.932 


■"57^ 

29,0Q0'“’ 

-150K222 " 

“’144.645‘ 

162.621 

180 8.4B8 ^ 

— - 

50 

29.250 

1507.138 

U 5*4?3 

161 .984 

1814.044 


59 

29.500 

1513.053 

1^6.199 

160.343 

181 9.595 


■ 60 

29.750 

“1S1 8^969 

‘U6.975 

159.208 

’1825.152 


61 

30.000 

152ii.889 

U7,753 

158.070 

1830,707 




TABLE 


2,15.33 ROSKAM-FILLMAM METHOD 

COIJf IGURAIIOH J (31 PAX/ 2/2) 

SENSITIVITY OF ZERO-UIFT 

DRAG TO TAIL CONE LENGTH 

OR IGI NAL C DOOHV 0,00957 


1 LC LC/D COOBHV 

(FT) 

1 15,000 1 .aS6 D,D089A 

7 15. "2 50 r.'8B7 D7009fl 

■3 1 5.500^ 1-918 Q, 00909 

A ■ 15.750 ‘ 1.949 ■ 0,00909 “ 

5 T 6.000 T;980 ““0.00910 

6 16,250 ' 2.011 0,00911 

' 7 ' 76-500 ” 2 .04 2 ’ 0.00912 " " 

8 16.750—^ 2.073 0.00913 

9 ' 17.000 ■ — 2.104 0-DD915 

10 - 17,250 2.135 0,00917 

"”11 I7v500 2.166 0,00918 

12 17.750 2.197 0.00919 

13 18,000 " 2.228 0.00921' 

- U 18.250 -2.259 0;00923 

1 5 1 8.500 2 .290 0,00924 

16 18,750 - ■ 2.321 0.00926 

“ 17 - 19.000 2.351 -0.00928 

18 19,250 2 .382 0.00930 

19 19,500 2,413 0-00932 

20 19,750 2'.444 0, 00933 

21 20,000 2.475 0,00935 

22 20,250 2.506 0,00937 

23 20.-500 2 .53? 0,00939 

24 


20, 750 


2,568 


0, 00941 


25 

21 *000 

2*599 

0.00943 

26 

21.250 

2.630 

0*00945 

27 

21.500 

2.661 

0.00947 

28 

21 *750 

2.692 

0.00949 

29 

22*000 

2,723 

0*00951 

50 

22*250 

2.754 

0.00953 - - 

51 

22.500 

2.785 

O.Q09SS 


22.750 

■2T816 

0.0d95 8 

33 

23.000 

2.847 

0,00959 

'■34 

■ 2 3.250 “ 

2 *877 “ “ 

■■■■Q. 00962 


35 23.500 2 .'90B 0.00964 


36 " 

23*750” 

"2,939 

0.00966 




37 “ ' 

24*000 

2 *970 

0.00968 

' " - - 




8 24 '.^SO 3VO0 1 'O; 0097 1 


39 

24.500 ^ ■ 

3.032 

0700972 

40 

24.750 

3.065 

0,00975 

4 r 

~^5*0D0 

3.094 

0700977 

42 

25.250 " 

■ ■'3.125 

(1,00979 

43 

25.'50l) 

■ '3.156 

0,00981 

44 ■ 

— 2 57750 

3':i87 

a;oo984'' 

45 

2 6,000 

■'3.218 

0*00986 

46 

26.250 

3.249 

0-Q098B 

4 7 ‘ 

26*500 

3.280 

0*00991 

48 

26.750 

3.311 

0.00992 

49 

27.000 

3.342 

0.00995 

50'* 

27.250 

3.373 

0.00997 

51 

27, 500 

3,403 

0*00999 

52 

27.750 

3.434 

0,01002 

53 

‘28.000 

3*465 

0701004 


54 


28.250 


3.496 


0.01006 



B.56 


t - .- ■ -i. ,' 51 .. ’ 


55 

28.500 

3.527 

0 . D 1 D 09 

' 56 

28.750 

3.558 

0,01011 

57 

29. 000 

3.539 

0-01013 

58 

29.250 

3.620 

0.01015 

59 

2 9 . ' 500 ' 

3.651 

0.01017 

60 

29.750 

3.682 

0.01020 

61 

30.000 

3.713 

0.01022 


- - 

• 

-- 


APPENDIX C FUSELAGEE SHAPE SLMGLATION PEGGRAM, 

FUSE 


Tliis appendix dociiinents the program FUSE and how to use it. 
Appendix Cl will describe the preparation of input data for both 
the simulation and design modes of the program. Default values 
for the design mode are also documented in this section. Appendix 
C2 will explain how the program was operated on the Honeywell 
66/60 timesharing systen and about the different operation options 
that were available. Appendix C3 provides a copy of the program 
listing. 

Cl Preparation of Input Data 

FUSE may be operated in either of two modes — the simulation 
mode or the design mode. In the simulation mode input data are 
read iron a separate disc file set up like a card data declc 
(Logical unit no. 03). In the design mode data are input by 
means of an interactive question and answer sequence. Each of 
these methods will be dealt with in tiim. To aid in understanding 
the ini'ut computer acronyms. Table C-1 has been prepared. 


Table C-1 FOSE Input Acronyms 


Acronym 

Variable 

Description 

Units 

H 

H 

Cabin height 

ft 

LED 

VD 

Fuselage fineness ratio 


LC 


Tail cone length 

ft 

LCD 


Tail cone fineness ratio 


HIT 

% l /^2 

Ratio of nose cone lengths 


im 

^N1 

Cone N1 length 

ft 

LN2 

%2 

Cone N2 length 

ft 

LU 

£ 

u 

Cabin (utility) length 

ft 

NEOKX . 

- 

Array containing number 
of lengthwise divisions 
per fuselage section 


NFOEX(l) 

- 

Number of lengthwise 
divisions for the nose 


NIDRX(2) 

“ 

Number of lengthwise 
divisions for the cabin 


NP0RX(3) 


Number of lengthwise 
divisions for the tail cone 


UJFUS 


Number of fuselage sections 
to be simulated; £ 3 


NRADX 

- 

Number of radial divisions 


PHll 

'*’11 

Cbne N1 taper shape param- 
eter (Side view) 


PH21 

'*’21 

Cone N1 bluntness shape 
parameter (Side view) 
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Table C-1 FUSE Input Acronyms (Cont'd) 


Acronym 

Variable 

Description 

Pm2 

*^12 

Cone N2 taper shape param- 
eter (Side view) 

PH22 

'*’22 

Cone N2 bluntness shape 
parameter (Side view) 

mn 

hi 

Cone N1 taper shape param- 
eter (Top view) 

PHY2 

^Y2 

Cone N1 blvintness shape 
parameter (Top view) 

PHCl 

hi 

Tail cone t£^r shape 
parameter 

PHC2 

h2 

Tail cone bluntness shape 
parameter 

Q21 

hi 

Coefficient of (x') for 
the cone offset equation 

002 

h2 

Coefficient of (x')^ for 
the tail cone offset 
equation 

(m 

%L 

Coefficient of x for the 
cone N1 offset equation 

QN2 

*^N2 

Coefficient of x^ for the 
cone N1 offset equation 

RBI 

^bl 

Cone N1 round-off radius, 
bottom 

RB2 

^b2 

Cone PE round-off radius, 
bottcxn 

RBC 

^bc 

Cabin round-off radius, 
bottom 

RBT 

^bt 

Tail round-off radius, 
bottom 


Units 



Table C-1 FUSE Input Acronyms (Oont'd) 


Acronym 

Variable 

Description 

Units 

RTl 

^tl 

Cone N1 round-off radius, 
top 


RT2 

^t2 

Cone N2 TOmd-off radius, 
top 


ETC 

^tc 

Cabin round-off radius, top 


ETT 

^tt 

Tail round-off radius, top 


W 

w 

Cabin width 

ft 

mL 

%l 

Half width of Cone N1 at 
^ ^ Sll 


2N0 

^0 

Vertical offset of Cone N1 
at X = 0 

ft 

ZNl 


Half hei^t of Cone N1 at 

ft 

zzc 

H 

Tail cone offset 


ZZN 

2^/H 

Eatio of nose cone heights 




Cl,l Simulation Made Input Data 


As was stated, in the simulation mode data are read from a 
disc file set up like a data card deck. In other words, each line 
of the file represents a different data card. The file must be 
sequential. Also only one simulation data deck may be run at 
a time. Table C-2 provides a guide for the "card" formats. 

Most of the input data required for FUSE are relatively 
sinple to derive from a drawing of the fuselage desired. The 
nose, however, and the offset cones in general can sonetimes 
present difficulties. In checking out FUSE several methods 
for deteimining the shape parameters for the nose were tried. 

The method that will be described in the following paragraphs 
consistently produced the best results. 

Figures C-1 and C-2 provide the top, side, and front views 
for the nose of the Gates Lear jet 35/36. These views will be 
used to determine the nose shape input parameters for FUSE. 

The length of the nose, is chosen to extend to the 

point vibere the fuselage cabin cross-section can best be said 
to become constant. Usually this will occur just aft of the 
crew coopartment . Hie elliptical cones are modelled such that 
the planform curves are perpendicular to the cross-sectional 
plane at the base of nose section (x = This should also 

be considered when the nose section length is chosen. The 
length of Cone N2, is determined by fairing the windshield 

curve down to the fuselage centerline as shown in Figure C-1. 
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Table C-2 (Continued) 


INPUT CARD NO. 5 - (7F10.0) 
PHYl PHY2 



INPUT CARD NO. 6-I7F10/0) 

YNl ZNl ■ ZN0 QNl QN2 QCl QC2 

I--- L - I I -| I -T 

1 II 21 31 41 51 


INPUT CARD NO. 7 -(7F10.0) 



INPUT CARD NO. 8 - (7F10.0) 




S.63- 


I 




NOTE : SCALE: 1/40 

NOSE SHAPE PARAMETERS (SIDE VIEW) 


N1 

= 10.58 

ft. 

^0 = 

-1.83 + 

N2 

= 5.63 

ft. 

'^11 = 

0.651 


= 5.48 

ft. 

11 

0.825 

'NO 

= - 1.83 

ft . 

>12 = 

0.598 

N1 

11 

CD 

O 

ft. 

^22 " 

0.749 


..00834X* 


CALC 


CHECK 



REVISED I date [Figure C-1 Gates Learjet 35/36 
Nose Shape Input Parameters 
(Side View). 


UNIVERSITY OF KANSAS 
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NOTE: SCALE: 1/40 


NOSE SHAPE 

PARAMETERS ‘ 

(TOP 

VIEW) 


II 

10.58 

ft. 

W 

5.48 

ft. 

II 

1 — 1 

2.61 

ft. 

♦yi ‘ 

0.644 

ft . 

■ +y2 “ 

0.805 

ft . 





IS REVISED [date I Figure C-2 Gates Learjet 35/36 

Nose Shape Input Parameters 
(Top View). 


UNIVERSITY OF KANSAS 
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The values for the Cone N2 shape parameters, (j»^ and 
are easily calculated by constructing the box and diagonals of 
Figure A-4 around the windshield cone. This has been acconplished 
in Figure C-1. The values for these shape parameters are als 3 
presented in Figure C-1. 

As Cone N1 is deformed and offset, the procedure becomes 
sli^tly more difficult. The first step is to fair the curve 
of the upper surface back to the base of the nose section. 

Again, in accomplishing this it should be ranaribered that the 
planform curves must be perpendicular to the cross-sectional 
plane at the base of the nose section. The centerline of the 
deformed cone is then constructed. By applying a parabolic 
curve fit routine the equation for is established. In this 
case a parabolic regression analysis was performed using a curve 
fit routine in the BASIC language. A listing of this pragram 
is provided in Appendix D. All offsets are referenced to the 
centerline of the fuselage. 

The half-height of Cone Nl, z^, measured as shown in 
Figure C-1. This will always be a positive value. 

The shape parameters, and be properly 

determined from the deformed cone. Instead the cone is projected 
onto a straightened centerline of length, as shown. The 

shape parameters are then calculated by the rectangle and diag- 
onals. The method for determining the gecraetry of a deformed 
elliptical cone may also be applied to the fuselage tail cone. 



Figure C-2 pictures the top and front views of the nose 
section. The front view is used to determine what the maximum 
width of Cone N1 at the nose section base must be. To avoid 
discontinuities at the base of the nose section, should be 
no greater than the y-value of the cahin a distance above 
the bottom of the cabin. This has been danonstrated in Figure 
C-2. The shape parameters and <j)y 2 are then calculated in 
the usual manner. 

The cabin input pai-araeters are relatively straightforward 
and should be easily determined. The length of the cabin is 
Jly - or that part of the fuselage which may be simulated by a 
cylinder. The width of the cabin is W aitd the height is H. 

The tail cone shape parameters should be conputed in a manner 
similar to that of the nose cones. The only differences are 
that the taper shape parameter, is calculated as the average 
of the top and side view taper shape parameters, and that the 
bluntness shape parameter, <|» 2 q. is the average of the top and 
side viev/ bluntness shape parameters. 

Table C-3 below provides a listing of the input data file 
for the Gates-Learjet 35/36. 

Table C-3 Lasting of the Data File for the 
Gates-Learjet 35/36 

3 7 7 15 10 

GATES-LFARJFT 35/36 — 530 PA?IFL£ 


IG. 55 

5.63 

10. 3A 

25.5 3 

5.4c 5.45 

0.651 

0.SS5 

0.598 

0.749 

0.60S 0*777 

0. 6A4 

0.005 



2.6t 

1 .50 

-1.90 

0.195 

-0.008 3 4 0.00217 0.00238 

1 . 

1 . 

1 . 

1 . 

I . 

1 . 

1. 

I . 
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Cl. 2 Design Mode Input Data 


Design Mode data are input by an interactive question and 


answer sequence. The roost effective way to describe this procedure 


is to follow the procedure for an exairple design. Such an exanple 


will be provided here. In this exanple underlined statsnents are 


those typed by the opeiator. The symbol '* + ” will be used to 


indicate a carriage return. Note that all of the following 


questions and answers documented here occur after the program 


has been coipiled and execution begun. 


DESIGN DR SINULfiTION MODE? 

ENTER D OR S 
=D % 

PRINT DIMENSIONS? Y FDR YES N FOR NO 
PRINT FUSELAGE COORDINATES? Y DR N 
PRINT TOTAL WETTED AREA? Y OR N 

^ ft 

^RINT PANEL AREAS AND RATIOS? Y OR N 

PLOT FUSELAGE? Y DR N 
=N 4 

BUSINESS JET DR PISTON COMMUTER? 
ENTER J OR P 

=£♦ 


The questions about printing data and plotting the fuselage 
are asked for either the design or simulation mode. They are 
included here for the sake of coninuity but will be discussed in 
detail in Appendix C2. The question as to wiiether a jet or piston 
aircraft is to be designed determines which default values are 
to be used when needed. These values will be discussed later. 
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"T '^~L 

= rC^*^ITn?ATIQw — gj FT C-CT^Fg 

r"7'!3r?^ OF PAKELS?- ~ ^ ’ V 

= ^rn t 

no. "of HAEIAL civ* each !t1EE? 

= irx ^ 

Elliptical qp. cihculap cpoee-efctichs? ''■' op ' 

■ COLE Ml : HQl-Pir-OFF PAtl I : TOPjBOTTGrr 

= - Ej iV -t 

COME N2 POUflC-OFF‘PAri ! : TOP^BOTTOM 
— 

CAPIP5 ROUfJE-OFF 'PAEII : -^OPjEaT-pM 
— • ^ j .3 

“^.IL POTtir-OF«‘ PAtl I : "CP^EO^TOM 
= . t ’ •0 ^ 


If the question, "ELLIPTICAL OR CIRCULAR CRDSS-SESCTION?" , 
had been answered with yes, FUSE would have assumed that all 
round-off radii were equal to 1.0. 


r^"'^EP:"'?ASEFT’IGEHE^ SEATS ^SREAETjPITCF 
=25,3^30 

SEAT WIDTKj AISLE WIDTH 
=20^ 1R ^ 

DESIRED INSIDE CABIN HFIOHT? CIM.) , 

^70 i 

FUSELAGE WIDTH = 6.9 S FT 

SHOULD H=W ? Y 0PM 

IWPUT.HTYO-M 
=N T 


l^SELAGE HEIGHT = 6.^^ 

INSIDE CA3IM HEIGHT = 70,00 IN. 


If the operator had decided that ;1 should equal W, then the 
desired inside cabin height is over-ridden. In this case the inside 
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cabin hei^t is caiputed according to tbe McDonnell Douglas method 
documented in Chapter 2, and is output. Due to a progranming 
oversight, if H is chosen to be input, the desired inside cabin 
height is output althou^ this may no longer be valid. 

•3A'7r:*GE. COMPA'^’TME^IT'? Y OP. 

=Y t 

=>> t- 

F'=> PH ! ! ^ PH? I w 1 S.., PV.22j 

— j j j j j ^ 

LED^LCD^ZZC' 

=•» t 

For this exanple a baggage canpartment was chosen to be 
included. The baggage conpartment is sized according to Section 
2.4. The shape parameters were chosen to be default values in the 
exanple above. The use of cannas alone inplies values of zero 
for each parameter. A value of zero tells the program to select 
the appropriate default value. The parameters above are all 
defined in Table C-1. 

With the data as input above the program will determine the 
aircraft size and conpute the node coordinates. 

The default values that have been mentioned are tabulated 
in Table C-4. These values have no statistical value. They were 
chosen during the programning stage rather arbitrarily. At that 
time it was felt that they did however generate designs which 
were representative of each of these two aiaplane tspes. 
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Table C-4 FUSE Default Values 

I 


T 

Variable 

Jet 

Piston 

11 

liLN 

2.5 

2.5 

1 


,75 

.75 


ZN0 



nn 


.65 

.65 


^21 

.825 

.85 

ts tj 

‘*’12 

.60 

.60 


‘*’22 

.75 

.725 



.64 

.65 


‘*’y2 

.825 

.85 

1 ) ■ 

IH) 

6.0 

6.0 


LCD 

4,0 

3.5 

" -i 

zzc 

.75 

.75 


'*’lc 

.60 

.60 

1 ■ 

‘*’2c 

.775 

.80 
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C3 Program Operation 


This program was written for the Honeywell 66/60 timesharing 
systen. Consequently, many of the coninands that will be described 
here will not be the same or perhaps even applicable on other 
systems. An attenpt will be made to e^^lain the reasons for each 
so that it will be possible for the reader to evaluate his own 
application. 

C3.1 Program Initialization 

As FUSE is written to be used with the Tektronix PLOT 10 
graphics software package, and because it is sometimes necessary 
to read from, or write to, other disc files, it is necessary to 
link the program file with several others before ccmpiling. 

This is accoiplished as described in the following paragraphs. 

Three files were required to operate the PIDT 10 package. 

The names for these on the University of Kansas system were 
ADEODT, ADEIN, and OBJECT. To be linked with FUSE these files 
had to be placed on the Available File Table (APT) for the KU-FRL 
project conputer account. Each project account is identified by 
what is referred to as a Project Identifier (PI). To place a 
file in the AFT three things must be given: the PI under which 

the file is stored; the catalog name (if any) under \diich the 
file is stored; and the filename. Then the conmand is: 

GET PI /CATALOG/FILENAME t 

If a logical unit number (L) is to be assigned to the file 
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when it is stored in the MT, then the conmand hecomes: 
GET PI/CATALOG/FILENAME’’L" + 


For the case of FUSE this might be acconplished as; 

SYSTEM? GET 76 84 WALLACE /PLOT 10 / ADEOUT f 
SYSTEM? GET 7684WALLACE/PL0T10/ADEIN + 

SYSTEM? GET 7684WALLACE/PLOT10/QBJECT + 

SYSTEM? GET 7102DESIGN/LEARJET’*03” + 

SYSTEM? GET 7102DESIGN/LEARCQORD"Q7" t 
This would gather the PLOTIO files, and LEARJET and lEARCOORD 
onto the AFT. LEARJET wo\ald be assigned the logical unit number ”03'' 
(to be read frcm) and IEAH(X)0ED, the logical unit number ”07” (to be 
written to). 

To load FUSE onto the AFT for use as a FORTRAN program the 
conmand was: 


SYSTEM? FORTRAN + 

OLD FILE OR NEW - OLD FUSE + 

And to conpile and execute: 

* RUMH *; ADEOUT; ADEIN; OBJECT + 

FUSE, in this case, is referred to as the star-level file; 
hence the * after the command RUNH. LEARJET and LEARCXX3RD are input/ 
output files to be referred to by the program and therefore are not 
needed in the run ccranand. Having carpleted these comnands, the 
program is ccsipiled and executed. 
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C2.2 Output Options 


There are six output options available for both the design 
and siMilation inodes: 

1) Print Dimensions 

2) Print Coordinates 

3) Print Total Wetted Area 

4) Print Panel Areas and Ratios 

5) Plot Fuselage 

6) File Data 

If dimensions are output, the result will appear either as 
in Table C-5a or as in Table C-5b, depending on \itiether the design 
or simulation mode is selected. Note that at the end of each 
printout, the operator is asked whether the data were satisfactory 
If not, he is given the choice of resubmitting his input or 
terminating execution. IVhen the dimensions are satisfactory, he 
may choose to file data. 

If coordinates are output, the result will appear as in 
Table C-6. The I and J value locate the node longitudinally 
and radially, respectively; and the X, Y, and Z values provide 
the coordinates in feet. 

Total wetted area of the fuselage appears as shown below; 

TOTflL UIETTEP fiPEfi = 769.52 
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Panel areas and ratios are output as shown in Table C-7. 

In this case, the I and J values locate the panels longitudinally 
and radially. The area of the panel is in ft^. RATIO J is the 
ratio of the area of that panel to the (I,J-1) panel, and RATIO I 
is the ratio of the area of that panel to the (I-1,J) panel. 

If a plot is requested, the fuselage is plotted on the graphics 
terminal similarly to the plots presented in Figures 2.43 and 2.44. 
The major difference is that the actual aircraft fuselage lines 
are not superimposed. 

If it is requested that data be filed, the fuselage coordinates 
are filed in a manner carroatible with the input data formats of the 
NCSU BODY program of Reference 20. 


C.19 



f 




cl 




f 


Table C-5a Dimension Printout for Design Ifode 


■■ V.- 

o 




T 


COHGn“P;ATI 0?: -- Si ’'A-!EE ^niTT-C'”" 




FU=F PP-OGRAJ: iJOCE — D 


fl-fV 

fai 


n VIS IONS ' , 




PAE-IAL' (HALF FrSFLAG-F)- 

!0 



LENGTHWISE 


« 

> * 

nose'; ’ 

7 



‘ CABIN- 

12 

*• ^ 

■j ■ 


tail; 

9- 

ii 

£ .• 

'• ' 

. • CABI W GONFI G UP ATI ON . 


T ' 

i; 

1 

no. or PArSFNGFFS 

?5 



seats pep -gw 

3 

. . 


■ ■ no. OF POPS . 

9 



F^'SELAGF GFOMET’^V 


i ^ 


nosE co?-TF :-:i , 




LENGTH 

10.96 



HALF 

P .AG 

* 


HALF HI ETT-: 

3.3S 

i 


*^KA?E PA'’A''!FT“’''^ 


L - 


PHI i 

c. 



PH2I 

K.S5C 

] 


PliV I 

G .‘6 SO 

1 .. 


PKY G 

C.G5F 

i ^ 


NOSE COME fia 

. 



LRIGTK 

6.26 



SH.APE PApAMET£’5C 


•, 


PHIS 

0 .69C 

i ' 

- 

PHC2 • . 

9.725 



CABIN 




LENGTH 

29.17 

I.: 


HEIGHT 

6 .4A 

* w 


■ -y, HI ETH 

6.92 

4 

11,' V 


TAIL CONE 




LENGTH 

23. 3S 

^ T" 


. SHAPE PAPAMETErS 




^KCl 

9 .Ae-G 



• PKC5 

G.S99 

,v 

TATA 

E ATI £ F/i cfc^Y ? Y 0'^ ?J 


dxr 

—V 



1 

7ILE 

=^T 

FUSELAGE EATA? Qt N 


> — 1 




11 
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Table C-5b Dimension Printout for Simulation Jfode 


.GA'^TS-LFAPjr*^ 35/36 — 5‘ 
‘fuse PrOGrAt! UOrE S 
iriuisious 

F A E I AL • AL F FT ' 5 ELAG E ) 

' LEAIGTK^TIS F 
. • fJOSE 
v:..., CASIFi;- 

t' /- TAIE.^ • , 

FUSELAGE GEOUETFY 


MOSF COME U1 
LEJiGTH 
HALF HEIGHT 
HALF HI r-ri-- 

cpp_p^ r-T>irpc 

PFl I ' 

PHGl 
PITY I 


iG-SS 

_I-:9G 

S.61 


? OSE COHE 


cuaPT pp'oaM 


5.63 


CAEI N 

'lehgtk 

- HEIGHT 
T' I ETH 


i0. 3A 
5.48 
5.43 


TAIL COME 
L ft; GTK 

SHAPE ?A PAM FT EPS 
PH c 1 


25*G3 


TATA S GP 


FILE ^’’‘''^ELA^E EA'^A? M 
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Table C-6 Example Coordinate Printout 


I 


I 

I 

J 

r $ 

Y 

2 

1 

1 

1 

0 . 

0 . 

-0.964 


1 

£ 

0 . 

0 . 

“0.964 


1 


0. 

0 . 

”0 .964 

k 

1 

4 

0. 

0. 

-0.964 


1 

5 

0. 

0. 

-0.964 


1 

6 

0. 

0. 

-0.964 

1 

1 

7 

0 . 

0. 

-0.964 

d 

1 

8 

0. 

0. 

-0.964 


1 

9 

0. 

0. 

-0.964 

d 

1 

10 

0. 

0. 

-0.964 

=1 

1 

11 

0. 

0. 

-0.964 



1 

1 .312 

0 . 

—2 .261 

■^1 

O 

o 

CL 

1 .312 

0.453 

-2,212 

2' 


1 .312 

0 .903 

-2 . 054 

M if 

O 

4 

1 .312 

1 .320 

-1 .762 


£ 


1 .312 

1 .60S 

“1 .326 

^ n 

iZ. 

6 

1 .312 

. 1 .662 

-0.303 


c! 

t* 

1 .312 

1 .460 

- 0 .329 



8 

1 .312 

1.113 

0.005 

7 - 


9 

1 .312 

0.731 

0.203 


£ 

lb 

1 .812 

0.359 

0,303 

*■' 

o 

11 

1 .312 

0 . 

0.333 


3 

1 

3 . 625 

0 . 

”2.300 


3 

o 

c 

3 . 625 

0 .676 

—2 .724 


3 

3 

3.625 

1 .336 

-2.431 


o 

4 

3.625 

1 .922 

-2.039 


3 

5 

3 .625 

2.305 

-1 .392 



6 

3.625 

2.334 

-0.643 


3 

7 

3 .625 

2 . 0 09 

0.010 


3 

'3 

3.625 

1 .503 

0 .453 


3 

9 

3 .625 

0.931 

0 .707 


3 

10 

3 .625 

0 .430 

0 .834 


3 

11 

3 ,625 

0 . 

0 .372 


4 

1 

5.437 

0. 

“3 .191 

1 ' 

4 

2 

5.437 

0.343 

-3 . 092 


4 

3 

5.437 

1 .663 

-2.778 


4 

4 

5.437 

2.331 

-2.212 


4 

5 

5.437 

2.319 

-1 .398 

II -fir 

4 

6 

5 .437 

2 .306 

-0.432 


4 

»■' 

5.437 

2.375 

0.290 

n rk 

4 

s 

5.437 

1 .760 

0 .796 


4 

9 

5.437 

1 .136 

1.081 

wt» 

4 

10 

5.437 

0.553 

1 .221 

ns-i 

4 

11 

5.437 

0 . 

1.262 


5 

1 

7 . 249 

0. 

-3.491 

C-. 

5 

2 

7.249 

0.991 

—3 .372 


5 

3 

7.249 

1 .942 

-2 .995 

/r;> 

5 

4 

7.249 

2.753 

— c' .8 2 1 


nr? 
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Table C-7 Example of Panel Area and Ratio Output 



<i jj;> 
HREfi 
RRTIGJ 
RRTIOI 


t I !• J> 
flPER 
RBTIDJ 
RflTIOI 

I ?J> 
fiRER 
PRTIDJ 
RRTIDI 


a »j> 
RREB 
RRTIOJ 
RRTIOI 


<i »j> 
flPER 
RRTIOJ 
RRTIOI 


A 

I 5 1 


£ ? 1 


3 > 1 


4f 1 


5> 1 



0.512E 

00 

0 . 1 OSE 

01 

0.1 4£E 

01 

0.1 70E 

01 

0.193E 

01 

I 

0 , 1 OE 

01 

0 . 1 OE 

01 

0 . 1 OE 

01 

O.lOE 

01 

0 . 1 OE 

01 

A 

O.lOE 

01 

0.£1E 

01 

0.13E 

01 

0.1£E 

01 

oaiE 

01 

T 

1j £ 


£j 2 


3.- 2 


4» £ 


5.- £ 


i 

0 .537E 

00 

0.112E 

01 

0.147E 

01 

0.1 75E 

01 

0.198E 

01 


0 . 1 OE 

01 

0 . 1 OE 

01 

O.lOE 

01 

O.lOE 

01 

O.lOE 

01 

T 

0 . 1 OE 

01 

0.21E 

01 

0 .13E 

01 

0.12E 

01 

O.llE 

01 


1» 3 


£ 3 


3j 3 


4» 3 


5, 3 



0.580E 

00 

O.llSE 

01 

0.1 53E 

01 

O.ISIE 

01 

0.£15E 

01 


O.llE 

01 

O.llE 

01 

0 . 1 OE 

01 

O.lOE 

01 

O.llE 

01 

uu 

oaoE 

01 

0.£0E 

01 

0.13E 

01 

0.12E 

01 

0.1£E 

01 


1? 4 


£j 4 


5 !• 4 


4- 4 


5j 4 



0 -6££E 

0 0 

0.1 £3E 

01 

0 . 1 57E 

01 

0.1 S4E 

01 

0.1 93E 

01 


O.llE 

01 

0 . 1 OE 

01 

O.lOE 

01 

0 . 1 OE 

01 

0.90E 

00 


O.lOE 

01 

0.20E 

01 

0.1 3E 

01 

0.1 £E 

01 

O.llE 

01 


1 j 5 


£5 5 


3j 5 


4> 5 


5- 5 



0 .64 IE 

00 

0.1 £3E 

01 

0.156E 

01 

0.1 S£E 

01 

0.£04E 

01 

^7 

O.lOE 

01 

O.lOE 

01 

O.lOE 

01 

0.99E 

00 

O.llE 

01 


O.lOE 

01 

0.1 9E 

01 

0.13E 

01 

0.12E 

01 

O.llE 

01 


1 j 6< 


£ ? 6 


3 > 6 


4i> 6 


5 j 6 



0.6££E 

00 

0.1£OE 

01 

0.15 IE 

01 

0.174E 

01 

0.£07E 

01 

«*■> 

0.97E 

0 0 

0.97E 

00 

0 .97E 

00 

0.96E 

00 

O.lOE 

01 


0 . 1 OE 

01 

0 . 19E 

01 

0 . 13E 

01 

0.12E 

01 

0.1 £E 

01 


1 j 7 


£ j 7 


3> 7 


4» 7 


5» 7 



0 .56£E 

0 0 

0 . 1 1 OE 

01 

0.1 37E 

01 

0.1 53E 

01 

0.214E 

01 


0.90E 

00 

0.92E 

00 

0 .90E 

0 0 

O.SSE 

0 0 

0 . 1 OE 

01 

kJ.lJ 

0 . 1 OE 

01 

0.20E 

01 

0.1£E 

01 

O.llE 

01 

0.14E 

01 


1 j 8 


£ !> 3 


3? 9 


4; 0 


5y S 


rt 

0.489E 

00 

0 .964E 

00 

O.llSE 

01 

0.1 44E 

01 

0.££0E 

01 


0 .S7E 

00 

0 . <58E 

00 

0 .S6E 

00 

0.94E 

00 

O.lOE 

01 


0 . 1 OE 

01 

0 .£0E 

01 

0.12E 

01 

0.1£E 

01 

0.15E 

01 

>;■ 

1 > 9 


£j 9 


3 > 9 


4j 9 


5» 9 



0.43 IE 

00 

0.S54E 

00 

0 . 1 04E 

01 

0.1 47E 

01 

0.££3E 

01 


0 .SiSE 

00 

0.39E 

00 

O.SSE 

00 

0 . 1 OE 

01 

0 . 1 OE 

01 

T 

0 . 1 OE 

01 

O.£0E 

01 

0.1£E 

01 

0.14E 

01 

0.15E 

01 

l-Lr 

1 1 0 


£ ) 1 0 


3.-10 


4.-10 


5 1 0 



0.40 IE 

00 

0 .796E 

00 

0 . 96 OE 

00 

0.150E 

01 

0 .££5E 

01 

t 

0.93E 

00 

0.93E 

00 

0.93E 

00 

O.lOE 

01 

O.lOE 

01 

Jir 

0 . 1 OE 

01 

,0.£0E 

01 

0.12E 

01 

U • 18E 

01 

0.15E 

01 
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C3 FUSE Listing 

This section provides a caiplete listing of the fuselage 
shape simulation program, FUSE. 
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ft |P€ ws ton i:f OfiX < ^ os 1 < zn* il ^ SNOS 2 C >rj 4 » i3'/3 ) IT 5 ros e 60# isn> 5 
KftH* I C 2 > # m I Lrn SX^DE F J^IILT S > #X X <1500 » # vrc i 5 00 > # 22 <1 500 > # 
itiftcsov5r>>Aoc5r>#3ni #«a ifi le < so > # rS 1 1 o«< 50 # rat i om 56 # lO) " 

fttAtt W1 #LN?^tU#^4N1 #«ftf #niV2#0rJ7 #;<C #rJOL OLL^Jr LOft 

ACiC HTWtt;Vi#LCD " ' “ ' " ’ 

CHARACTER KOftf r^ODClC^TYf ErTYPCKrlOOfrOf 


OAIA < &EFAOLrtWll irial #1»)yi,75*.75#.65#.H25 
i , <V4 # , 8 ? S # ^ , C # . 70 # > 6 7 7S # 6 . C / 
ft At A < ft€ f AOt T < 7 ^ I > # I r 1 , M ) / 1 , 75 # . 75 # . 65 ^ . R STTftTTmT" 
A aftS »■ AS # 6 # #■# 70 #,»6 r » 3 • S / ^ 

f All <Kt^Sl?HX> /cost 

rbftCx#A#>^#»r#o) a 0^ tn ,-<x/Ar**p)» • < i . / o > > 


‘ ITYPfl 

Ipt^t 


||PR»1 

rco^j 


fAPtA«{f 

rpA*o 


I FILE •'I 

PRINT : ••Of sisti OR SIHOLATION HOOE?** 


PHIHT:'*E»1iER D OR S*" 
WEADrPOftCK 


PRIMll-PRINT OINfNSIfiNS? Y fOR YES^ » FOR NO^* 
Rf Aft; or : 


|F;(OT»EQ.IOUT>IPR» 

PRINT;-PR|!PI fOSftAGE COORftIMATES? Y OR 


P€:AO:Or 

If tor.eo. iiOuT>ico<i 


PRINIs'*PRi|<JT TOTAL UETfEO AREA? Y OR N 
REAOiOt . .. • 


IFtOT.EO.IDUri lAREAitl 

PRINT :**rRIft!T PANEL ARC AS AND RATIOS? Y OR «-• 
REAOiOT ““ 

t f COT , Eg, iOU T X | PA«1 


PRlKf :-PLOT FUSELAGE? Y OR N 
REAP; Of ^ . 


IFCOI ,Eftp|OUTy|PL«r 

F t ■^OftCK^ EO^ NOPE > GO TO 90 

NFllS^i 

PR t tir;** ^llSINESS JET OR PISTCM COt RQ TER?** 

PRINT {'•FMTER J OR P** 

REA P; Tt PCX . 


I Ft T YPCK •EO.T YPE> I T YPY? 
PRINT:;*T | tiE ?** 


REA DC S# 1 1 ) C f T I TLE C 1) , t «1 « 1 A ) 

FR| NT Of PANELS?** 

READ: P?fi 

PRINT ;*•»(), OF RAOiAL Oltf. EACH SIDE?" 


READrNRAOX 
IPNi * E N L 


PR|KTs"FLL|PTICAL CH circular CROSS-SECTIONfr Y OR K" 
READ; or 


IMCr.NE.fGUnGO TO ’5 

RTt sT ^ ’ . . 


R01 -I • 

RT2M, _ 

R07-1, 




CUOtJOlMO 

fc OgLii OZO 

COjSG'USir 

Ltjyacft4g 

ciijooosg 

Qugoooftg 


cyggoc9j 
EQ30QT03 
003001 T'J 
„(? 33001 ?g_ 
£u3Q&133 
CO0QG14O 


CQ3001 53 
C3300160 


003 OG 173 
OOQOCIdJ 


C030G193 

00300203 


00D3Q21G 

C0caj220 


00330233 

00300240 


CaCG32S3 

10300263 


TOntltiiPm 


000032 DO 


OOOOG29J 
00300500 


3CDGj5ia 
t 03 00 123 


D03Q0550 

^003 341 

G3 j'uG 3 53 
10330363 


00C03373 

C033Q3S0 


00300390 

e3J}004G3_ 

COOGuAIO 

CQOOJ423 

00300453 

000:00440 


0030C45C 

CnCD:46L 


00300473 

000 03433 

00301490 

OGtOOSOJ 


00000513 
00300523 


G3003553 
00330540 


0C30C55J 

00000560 


0033G573 
G0DDG5SJ 
t 0^00593 










|»rKnmM»f4 V ■xnv|(»»r4f..i L^-ssuv^ l^r^ytitjir!4 








i»cti»oir*fsB»,tf, fF«*nnotfSi8bj»"l' ■ - fftnUTTyJ 

Ltl»CiBOw«J»l/tg.»tL 0aU012L>.l 




ttF«OtlLN#/tN#tMr 

— tnrrarni 

Oti}QI220 





|FCttM.F».r.HLK*(lteff ULttllV(>*t > 
irCt2».F(l.r.t7II.'<0tF««IL1 CI1fl>*?> 

U0JJt241 





F«iNT:"EH;TF(» P)ill#f>H?t#Pril2#PH??#PIIVl#KHV?*‘ 
BCBPrPHl tVpH»l»«»>i;1 ?#PH??»PM¥t »PMT2 

0030 126:: 





TfTP«TiTf fl.C , » PMl t ■OF. F BULK 1 f tP» J» 
I M PU ?1 . Ea.c . »P!*>1 •6E r »0L t I I I VP #4 > 

smryn — 

3033 12 S3 





■ IFIPHii . E 1 , 0 , 1 PM 1 2 «»C F •Util I T » #51 

tM PM2?. FA. 1 , > PH??>0E F 4UL II 1 1 »P#M 

■ cJ3tii2^3 —— 

10301303 



- 


■ f I fHl VI . E ). C . Vph Vl««Ef Alii, till V P * 7 1 
|F|PilV?.E4.b.)p|iv>BSErAU(.M|fVPf5) 

tojirtJin 

10301323 





Pill nr t •»€ fj IE H "laii»L”cb# t if’*^ 
RE«»tL.lP#iCO#IFC 

'otijytlii 
00301 34!) 





1 FCL^«.Eb.r . ILaP'OEFAUtTI II VP#9> 

1 FI LC0.E n.r . )IC0*DE FAULT 1 1 1 VP# 1 1 1 

&O301350 
C 030 1360 





IME*C.Ea.r.»IEC«l)!EFAULTI|IVP#1'’J 
PR|Nlt**EtllEII P»CI#PMt?" 

r!j?f'*T2B 

0O3OV5SJ 





ReAPiPMCt#FHC? 

IMP|ICl.ER.C.>PMC1>0EF*ULIIITVP#t1 ) 

00301 3V0 ^ '■“ 

C0001400 





IFCPIIC2.E<).0.}PIIC2«I)E~FAULIIITVP#12I 
I F 1 is b . E a * f . 1 2 •I'* » 12 i « -1 . > * H / ? . 

QU3G141J 

C033I423 



- 


CO TO 95 

S3 coritiriuE 

ooiiFfcTS “ 

00301440 





C0U»IIUS*5r.>*li.*RT2)*l.9t*U*l2.1 1/12. 
IFfCRU.CI.(.>GO 10 85 

003014513 . " 

00001460 





CO TO SA 

CC301473 

D03Q1433 




a 

S5 Cit2aC»IHt>9.47*RT?*R*saRTt1.>{Ct1M/'{.R4*RT2*U>>«{CaU/l. 
«>0.A7*M*I1.>RI2> 

94#|tf2*^)>) CQ3C1490 

COGOnOQ 



1 

to 

VI 

'll viic • A H5 1 C MU 7- C B 7 1 

Lfi?*CRLl 7ln.-ll7,iVRC/l.9««Hn *MK2))*Mt.7NN2)> 

t533TJiJ5 ' ' 

00301510 





iUI*LL«*L«? 

{n1«22N*M/2. 

ooJMTSo ^ 

C 0301 540 

'em 


- 


2»>i«2Bl-H75, 
QN2s|{Nb>irv>7{tKl*LM 1 

iSlffUTT 

00301560 





aNl«>?,*UR?*LR1 

VHla('>/2.t«SQRTI1 .•l2.*ItlK7H)*(7.*2DX/Hn 

i»03ai57» ' 

GOwOlSSU 





tC»ti»0*<R*M>>?.-LU-LNl 
Oat '«»»>/?. 

«b331595 

00301600 



' 


I F 1 IL C /M .li.iC BH Cal C D* D 

oci«:. 

00001617 ” 

00301623 






I'O^uTbii 
0030 1643 





KNtatIL 

*a|.T.«LBl F?.*li|*|^5*lC)7n. *KRL> 

,03301653 ■ ■■ 

C03C1663 





»FbR«tl>«;LM/V 

NFO<l«C2)atU/{1.5*V> 

C0301673 
C 030 16 S3 





mf 6*<< I ilaRL'BF O-'A (Vl'RFiiftXI 2) 
50 10 95 

uojbuVo 

OC'iC173J 





or REROi l#i F.'lFUS#KF0Sv#RRAtl * 
READII#11 MlirHFt|)#|at,tn> 

cc::i\7ij 





«E* Dl 3#? lint #litT*:lu#L C#U#H 
RE«DI5#?)fHl1#FM2l #FMI?#PH22#PHCt#FHC2 

C03Ll1?t3 

nU30124i) 



y 


«E*Di i#>>pivt#pitvK ~ 

RE * 01 3 #2 > VM # 20 1 # 2 1’,? # OMI # CM ? # OC 1 # t| F ? 

C0331750 
00331 760 




(ft > Sf t t 

jo3:i7?: ' 




V 
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01.1 COWSHftf Bl t Lff**** «WHl I 

"?ii 1 C 6«IS HP (ptt1>«PH?> » fi^oni »H lili tHHii 

c*it . coH|«ipjjpMCi eP»c2*f*.aqi #fic.sici 
' oil co«iiip|piTti#i*iiV?Vr.o(Jl »M»I *K* 

|f«.‘<ODC<.riE.HOP':> C O fO ^0 

lo'v? 

tf»ITt < 6 r 1 ?> ( l 

3 il fF Cft#UrxQ«CK 

wqirM 6 «tS)N««Dx,(Kro«xfj>*i«t* 5 > 

if('WftO^jEb.*40i>E>^'c^ to lOOl 

UR I fc <4*1 IS) 1P«K*I SF#I ROU_ 

uRlff <A*r>utif *r»t #vMi| #piii i »pm? « *piiri*i»iiY 2 

URI fMA*t ‘t)L*t2>PHl2*PU?^J^#H«U«LC *PHC1»PHC2 

~PRl»ll"oiti “siiTisf •ctaiY? y' or «" 

RfROtOf 


tFCOf.E 9 .IOur>SO to 97 

PRI‘tti'*RF Sm>.'iIT DRTO? Y OR U" 
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|f|«b«<0A1 •.LT.fOl) 00 10 A3 
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00305753 
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0 * UC : ^ ■ 
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w * Ul : ■ 

m to 7^ 

00305573 

90395433 




RP ^ » Cw| ♦ F?'»*C.5 
7C AT IMIJf 
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10335350 
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50 13 V': : 

0U3055?0 
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APPENDIX D PARABOLIC LEAST-SQUARES 


CURVE-FIT ROUTINE 


This appendix presents a listing for a time- 
sharing program in BASIC to compute the coefficients 
for a least-squares parabolic curve-fit of the general 
form: 

T7 L 1 2 

Y == + a^^x + agX . 


D.l 



0010 HIM D< 1 0 j£> ,i3> jBCS::’ jC'-S 

00^0 PPIHT "NUMBER DP DRTR POINT 
0030 INPUT N 

0 040 PRINT "INPUT DflTfi 5 X-.Y" 

0050 FDR 1=1 TO N 
0 06 0 INPUT D a j 1 > j n 'll !■ £ ) 

0070 NEXT I 

0080 FOR 1=1 TO N 

0090 PRINT D<I ?1> flKI j£> 

0100 NEXT I 
0110 X 1=0. 

0120 X£=0. 

0130 X3=0. 

0140 X4=0. 

0150 Y1=0. 

0160 Y£=0. 

0170 Y3=0. 

0180 Y4=0. 

0190 Y5=0. 

0200 S1=0. 

0210 S£=0. 

0220 FDR 1=1 TD N 
0230 X1=X1+D< I •.!> 

0240 x£=x£-i-ri I !. 1 :? ♦n a i > 

025 0 X3=x3+n c: 1 7 1 ;> +D I ) 1 > +d a ? i > 

026 0 X ■■'•=X4+Ii ( I j 1 •t4 
0270 Y1:=Y1+IKI j2> 

0280 Y£=Y2+B C I j 1 > ^Ii < I » £ > 

029 0 Y3 = Y3+B <1 j 1 > ^IK I > 1 > ♦D a , 2 > 

03 0 0 Y5=Y5+B < I ? £ > < I j 2 > 

0310 NEXT I 

0320 Y4=Y1-^N 
0330 E<1>=Y1 
0340 B<£>=Y£ 

0350 E';:3>=Y3 
0360 fl<l)l>=N 
0370 P':;2!>1>=X1 
0380 Pc:i7£>=Xl 
0390 8^.2 72.^— X2 

04 00 R';1j3>=X£ 

0410 R>;3sl>=X2 
0420 fit:£73>=X3 
043' 0 fi 3 5 2 =X)3 
0440 R ( 3 7 3 ) =X4 
0450 MRT R=INVCR> 

0460 NRT C=R+B 


I 

I 

I 

T 

147 0 S£= liL Y5-I- 1 > 1 -L ■; c > t iii > ♦’r‘3 •••• N 
i4S0 SS=:?£-^S£ 

1490 FOR 1 = 1 TO H 
15 Cl 0 Y= V+ ai< I !. 2 .. -Y4 > t£ 
i!51Q NEJa I 
0I520 S1=Y.-'N 
0530 R£=1.-S2/S1 
0540 PRIMT "1=10 ="jC<1> 

0550 PRINT "fll ="»i:<:£> 

0560 PRINT "R2 ="jC(:3> 

0570 PRINT "CDEFF. OF DET . jR£ =" ?R2 
0530 STOP 
0590 END 
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